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This  report  describes  the  results  obtained  during  the  two  years  of 
funding  extending  from  1/88  to  1/90  in  support  of  our  project  entitled  "The 
problem  of  Robust  Compensation  for  Systems  with  Unmodeled  Dynamics."  The 
project  has  consisted  of  several  lines  of  research  which  are  quite  distinct, 
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theory.  Our  goal  has  been  to  explore  the  problem  of  designing  feedback 
control  systems  that  are  insensitive  to  the  presence  of  high-frequency 
dynamics  not  accounted  for  explicitly  in  the  mathematical  model  of  the  plant. 
Some  of  our  previous  work  suggests  that  it  is  possible  to  design  controllers 
which  simultaneously  stabilize  a  given  nominal  system  as  well  as  a  large  class 
of  small  singular  perturbations  of  the  system.  Attached  are  six  papers 
summarizing  work  which  has  been  supported  all  or  in  part  by  the  present  grant 
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Abstract 


This  report  describes  the  results  obtained  during  the  two  years  of 
funding  extending  from  1/88  to  1/90  in  support  of  our  project  entitled  "The 
problem  of  Robust  Compensation  for  Systems  with  Unmodeled  Dynamics."  The 
project  has  consisted  of  several  lines  of  research  which  are  quite  distinct, 
but  which  show  great  promise  toward  combining  them  into  a  single  comprehensive 
theory.  Our  goal  has  been  to  explore  the  problem  of  designing  feedback 
control  systems  that  are  insensitive  to  the  presence  of  high-frequency 
dynamics  not  accounted  for  explicitly  in  the  mathematical  model  of  the  plant. 
Some  of  our  previous  work  suggests  that  it  is  possible  to  design  controllers 
which  simultaneously  stabilize  a  given  nominal  system  as  well  as  a  large  class 
of  small  singular  perturbations  of  the  system.  Attached  are  six  papers 
summarizing  work  which  has  been  supported  all  or  in  part  by  the  present  grant 
and  which  either  have  appeared,  have  been  accepted,  or  have  been  submitted  for 
publ i cation . 
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This  report  describes  the  results  obtained  during  the  two  years  '6f, 
funding  extending  from  1/88  to  1/90  with  regard  to  our  project  entitled  "The 
Problem  of  Robust  Compensation  for  Systems  with  Unmodeled  Dynamics."  The 
project  has  consisted  of  several  lines  of  research  which  are  quite  distinct, 
but  which  show  great  promise  toward  combining  them  into  a  single  comprehensive 
theory.  Our  goal  has  been  to  explore  the  problem  of  designing  feedback 
control  systems  that  are  insensitive  to  the  presence  of  high-frequency 
dynamics  not  accounted  for  explicitly  in  the  mathematical  model  of  the  plant. 
Some  of  our  previous  work  suggests  that  it  is  possible  to  design  controllers 
which  simultaneously  stabilize  a  given  nominal  system  as  well  as  a  large  class 
of  small  singular  perturbations  of  the  system.  We  say  such  a  controller  is 
"robust"  with  respect  to  the  certain  class  of  unmodeled  dynamics. 

Attached  are  six  papers  summarizing  work  which  has  been  supported  all  or 
in  part  by  the  present  grant  and  which  either  have  appeared,  have  been 
accepted,  or  have  been  submitted  for  publication.  Two  have  already  appeared 
in  the  IEEE  Transactions  on  Automatic  Control,  one  is  scheduled  to  appear  in 
the  same  journal  in  May,  one  has  appeared  in  the  Proceedings  of  the  27th  IEEE 
Conference  on  Decision  and  Control,  and  one  is  presently  under  review.  In 
addition,  a  summary  page  is  included  describing  a  paper  presented  at  the  SIAM 
Conference  on  Control  in  the  90’s  (5/89,  San  Francisco).  We  also  are  in  the 
process  of  writing  a  journal  article  describing  the  work  carried  out  in  the 
the  final  stage  of  the  project.  The  manuscript  should  be  available  within  the 
next  two  months.  Part  of  the  research  effort  over  the  past  two  years  has  been 
carried  out  by  my  graduate  student  Mingde  Tan  as  part  of  his  Ph.D.  research. 
His  dissertation  should  be  available  sometime  this  spring. 


Reference  11]  contains  work  which  was  partially  carried  out  during  the 
funding  period  and  which  was  instrumental  in  establishing  the  direction 
outlined  in  the  original  grant  proposal.  Our  basic  idea  was  to  explore  the 
role  of  parasitics  in  the  performance  of  automatic  control  systems,  without 
having  to  resort  to  explicit  representations  of  specific  parasitic  effects. 
It  was  our  desire  to  develop  a  comprehensive  theory  of  compensator  design 
which  would  guarantee  performance  in  the  presence  of  a  large  class  of  possible 
parasitic  effects.  Reference  [1]  contains  preliminary  results  directed  at 
this  goal. 

Along  somewhat  different  theoretical  lines,  but  with  a  similar  class  of 
engineering  problems  in  mind,  references  [2]-[4]  constitute  initial  attempts, 
carried  out  with  the  aid  of  a  colleague  here  at  Wisconsin,  Professor  Chris 
DeMarco,  to  characterize  the  geometric  structure  of  the  class  of  system 
perturbations  under  which  a  control  system  retains  stability  and  perhaps  other 
performance  characteristics.  Contained  in  [2 ] - [4 ]  is  a  thorough  treatment  of 
the  case  where  system  order  is  constant  (the  nonsingular  case).  Our  more 
recent  work  may  be  viewed  as  extending  these  results  to  the  case  where 
parameter  variations  can  cause  changes  in  system  order  (the  singular  case). 

Reference  [5]  summarizes  the  work  carried  out  by  my  graduate  student, 
Mingde  Tan,  and  me  over  the  initial  phase  of  his  dissertation  research  (the 
first  year  of  the  funding  period).  Our  principal  idea  was  to  study  the 
effects  of  small  system  perturbations  on  internal  closed-loop  stability. 
Ultimately,  we  wish  to  examine  many  of  the  recent  robust  control  theories  of 
Vidyasagar,  Zames  and  Francis,  Stein  and  Doyle,  and  others  in  terms  of 
internal  closed-loop  system  behavior.  For  example,  suppose  a  plant  model  is 
given  and  a  corresponding  compensator  is  designed  using  some  methodology  such 
that  the  closed-loop  configuration  is  input-output  stable.  The  design  methods 


of  the  researchers  just  mentioned  typically  guarantee  that,  for  a  certain 
class  of  perturbations  of  the  plant  (and  sometimes  the  compensator),  the 
corresponding  perturbed  closed-loop  system  is  also  stable  in  an  input-output 
sense.  It  would  be  highly  desirable  to  know  whether  the  perturbed  closed-loop 
system  is  internally  stable  as  well.  From  our  perspective,  it  is  especially 
crucial  that  the  compensator  be  insensitive  to  perturbations  corresponding  to 
unmodelled  high-frequency  dynamics. 

In  order  to  formulate  this  problem  precisely,  one  must  first  develop  an 
understanding  of  how  internal  system  structure  is  affected  by  perturbations  of 

the  transfer  function.  To  this  end  we  have  developed  a  "perturbational" 

analogue  of  the  standard  state-space  realization  theory  for  rational  matrices 
In  our  theory,  families  of  rational  matrices  are  considered  (either  convergent 
sequences  or  continuous  parametrizations) ;  it  is  desired  to  find  corresponding 
(convergent)  families  of  state  equations  which  realize  the  given  transfer 
matrices.  We  have  obtained  results  that  show,  for  example,  that  every  such 

family  of  transfer  matrices  has  a  realization  and  that  "minimality"  of  a 

realization  can  be  related  to  both  the  degrees  of  the  given  :ational  functions 
as  well  as  to  controllability  and  observability  of  the  realizations 
themselves.  In  short,  we  have  succeeded  in  developing  the  perturbational 
analogues  of  the  standard  results  concerning  realization  of  a  fixed  rational 
matrix.  This  body  of  results  is  among  our  main  accomplishments  in  the  past 
two  years. 

Based  on  our  understanding  of  the  fundamental  issues  surrounding  internal 
realizations  of  perturbed  transfer  functions,  we  have  spent  a  large  portion  of 
the  past  year  exploring  the  relationships  between  robust  input-output 
stability  and  robust  internal  stability.  We  have  succeeded  in  establishing 
simple  conditions  under  which  internal  as  well  as  input-output  stability  is 


( 


robust  to  parasitic  model  uncertainty.  These  results  form  the  second  half  of 
Mingde  Tan’s  dissertation  and  are  presently  being  organized  for  journal 
publication . 

In  addition  to  the  work  with  Mingde  Tan,  I  personally  have  continued  the 
effort  I  began  in  1986  which  addresses  questions  similar  to  those  described 
above  related  to  internal  behavior  of  robust  closed-loop  systems.  The 
culmination  of  this  effort  so  far  is  the  paper  [6]  which  shows  that  robust 
design  methodologies  necessarily  must  incorporate  some  internal  system 
information;  the  exact  form  of  the  "minimal"  internal  information  is  as  yet 
unknown,  A  recent  breakthrough  has  been  summarized  in  [7]  which  was  presented 
in  May  1989  at  the  SIAM  Conference  on  Control  in  the  90's  in  San  Francisco. 
(There  were  no  published  Proceedings  for  this  conference.)  For  a  certain 
large  class  of  linear  systems,  it  is  now  possible  to  precisely  characterize 
the  family  of  singular  perturbations  under  which  closed-loop  stability  is 
retained. 

Overall,  our  two-year  research  effort  has  been  fruitful,  answering  many 
important  questions  and  leading  to  new  ones.  In  the  final  analysis,  the 
entire  body  of  results  we  have  generated  is  not  as  coherent  as  we  original 
hoped  for;  however,  we  feel  that,  given  the  time-limit  set  for  the  project, 
the  progress  made  in  each  research  direction  more  than  compensate  for  this 


fact. 
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Toward  a  Theory  of  Robust  Compensation  for 
Systems  with  Unknown  Parasitics 

J.  DANIEL  COBB,  member,  ieee 


Abstract— We  consider  the  problem  of  designing  a  robust  compensator 
based  on  a  plant  model  with  order  uncertainty.  The  uncertainly  is 
characterized  mathematically  as  a  class  of  generalized  singular  perturba¬ 
tions  of  the  plant.  This  paper  considers  the  case  of  static  compensation.  A 
necessary  and  sufficient  condition  is  established  under  which  actual 
closed-loop  behavior  is  close  to  that  predicted  by  the  plant  model  under 
sufficiently  small  singular  perturbations.  The  condition  is  shown  to  be 
generic. 


I.  Introduction 

THE  problem  of  robust  compensation  may  be  roughly  stated  as 
that  of  designing  a  good  controller  for  a  given  physical  system 
on  the  basis  of  a  model  which  contains  less  than  complete 
information  about  that  system.  The  resulting  closed-loop  configu¬ 
ration  should  exhibit  reasonable  performance  in  spite  of  the 
uncertain  aspects  of  the  system.  In  the  strictest  sense,  every  model 
contains  uncertainty;  hence,  any  good  controller  design  should 
address  the  issue  of  robustness. 

Among  the  many  types  of  robust  control  theories  appearing  in 
the  literature  is  the  asymptotic  approach.  Typical  results  in  this 
area  guarantee  reasonable  closed-loop  performance  under  suffi¬ 
ciently  small  perturbations  of  a  nomial  model  (e.g.,  variations  in 
the  coefficients  of  a  single  differential  equation).  Although  only 
local  in  nature,  such  results  are  often  a  first  step  in  developing  a 
global  theory  where  an  explicit  characterization  is  attained  for 
classes  of  systems  which  can  be  simultaneously  compensated.  The 
results  of  this  paper  fall  into  the  asymptotic  category. 

It  is  possible  to  view  most  asymptotic  robustness  theories  within 
a  common  mathematical  framework.  Let  (?,  Q,,  and  3  be 
topological  spaces,  and  let  <R  C  (P  x  Q,  inherit  subset  topology. 
(P,  Q.,  and  3  correspond  to  the  sets  of  all  possible  models  of 
plants,  compensators,  and  closed-loop  systems,  respectively.  The 
topologies  on  (P  and  Q,  are  chosen  so  that  small  perturbations 
characterize  measurement  error  inherent  in  developing  each 
model;  small  perturbations  in  the  topology  of  3  reflect  tolerable 
closed-loop  performance  error.  If  (R  is  interpreted  as  the  class  of 
all  plant-compensator  pairs  which  lead  to  closed-loop  systems  that 
are  well-defined  and  which  satisfy  any  additional  constraints 
present  in  the  design  problem,  we  may  naturally  define  the  loop- 
closing  map  Q:  (R  -*  3  which  takes  each  plant  and  compensator 
into  their  corresponding  closed-loop  configuration.  Many  robust¬ 
ness  questions  then  reduce  to  that  of  finding  the  points  of 
continuity  of  Q.  In  other  words,  we  wish  to  characterize  the  class 
of  alt  plant-compensator  pairs  such  that  small  perturbations  of 
each  pair  result  in  small  perturbations  in  the  closed-loop  system. 


Manuscript  received  July  6,  1987:  revised  May  25.  1988.  This  paper  is 
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We  now  examine  various  existing  theories  which  lie  within  the 
asymptotic  framework.  The  most  obvious  body  of  such  results 
centers  around  the  well-known  fact  that,  for  state-space  models, 
the  parameters  of  the  closed-loop  system  are  continuous  functions 
of  the  open-loop  plant  und  compensator  parameters.  For  example, 
if  we  let  (P  be  the  set  of  all  matrix  triples  £  =  {A,  B,  C)  arid  Q, 
consist  of  all  feedback  matrices  K,  and  if  we  combine  £  and  /fin  a 
standard  way,  then  (R  =  (P  x  Q.and  3  consists  of  triples  C  ({,  K) 
=  (A  +  BKC,B,Q.  Adopting  Euclidean  topology  oh  ,<J\  Q,,  and 
3,  it  follows  that  6  is  continuous  everywhere,  i.e.,  every 
compensator  is  robust  relative  to  every  plant,  Orie  immediate 
consequence  of  this  observation  is  that  closed-loop  eigenvalues 
are  continuous  functions  of  plant  and  compensator  parameters; 
hence,  every  stable  closed-loop  configuration  remains  stable 
under  sufficiently  small  parameter  variations.  These  facts  are  used 
routinely  in  many  control  system  analyses  without  explicit 
mention.  It  should  be  noted,  however,  that  the  perturbations 
considered  here  do  not  alter  either  plant  nr  compensator  order. 
Therefore,  this  approach  alone  is  inadequate  when  dealing  with 
order  uncertainty. 

The  main  body  of  existing  results  that  does  deal  with  order 
uncertainty  in  an  asymptotic  setting  can  be  broadly  termed 
singular  perturbation  theory  (see  [11-13]).  Here  a  typical  analysis 
treats  a  parametrized  system  of  the  form 

[^#/H£  aHi]  ■ 

y=[Q  c2 1*  0) 


with  A 22  stable  and  seeks  to  achieve  some  closed-loop  perform¬ 
ance  criteria  for  all  sufficiently  small  e  a  0.  (In  this  case,  we 
might  take  (P  =  |0,  »).)  A  major  drawback  with  this  approach  is 
that  explicit  knowledge  of  the  parasitic  structure  giving  rise  to 
order  uncertainty  is  assumed.  If  more  than  one  perturbation  (1) 
need  to  be  considered,  serious  problems  may  develop.  For 
example,  the  system 
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is  nominally  (e  =  0)  unstable,  but  can  be  stabilized  with  the  static 
compensator  u  =  2 y.  The  perturbed  system  (e  >  0)  is  also 
stabilized  by  the  same  compensator  for  sufficiently  small  e. 
Setting  e  =  0,  premultiplicaiion  of  (2)  by  the  matrix 


M= 


I  0  0 
0  1  2 
0  0  1 


yields  an  equivalent  system  equation  which  may  in  turn  be 
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perturbed  according  to 
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In  this  case  the  compensator  u  =  ly  yields  a  perturbed  closed- 
loop  system  having  a  pair  of  eigenvalues  and  \2l  with 
Re  X„  -»  +  oo  as  e  -*  0 ' .  Such  divergent  behavior  does  not 
coincide  with  any  reasonable  definition  of  small  perturbations  in 
3.  We  may  therefore  conclude  that  examination  of  a  single 
parasit: ;  effect  is  in  general  not  sufficient  to  guarantee  robustness 
of  a  compensator  with  respect  to  other  order  uncertainties. 

Additional  singular  perturbation  results  include  the  multiple 
time-scale  extensions  [7}  and  [8]  and  the  robust  compensation 
theorems  of  [5],  Multiple  time-scale  techniques  suffer  from  the 
same  drawback  as  single  time-scale  analyses  based  on  (1)  in  that 
they  assume  an  explicit  knowledge  of  parasitic  structure.  Also, 
much  less  is  known  about  the  £-dependence  of  the  time  response 
of  multiple  time-scale  systems  than  in  the  single  time-scale  case. 

In  (5]  it  is  shown  that  any  compensator  having  a  strictly  proper 
transfer  function  matrix,  which  stabilized  (1)  with  t  =  0,  also 
stabilizes  (l)  when  e  >  0  is  sufficiently  smull.  Furthermore,  it  is 
shown  that  the  corresponding  family  of  closed-loop  transfer 
matrices  converges  uniformly  on  compact  subsets  of  the  right-half 
complex  plane  as  t  -»  0  * .  These  results  thus  provide  a  means  for 
robustly  compensating  a  system  in  the  presence  of  a  large  class  of 
possible  perturbations.  One  drawback  to  this  theory  is  that  only 
single  time-scale  systems  (1)  are  treated.  In  practice,  a  much 
larger  class  of  perturbations  may  be  required  to  model  all  relevant 
effects.  Additional  problems  are  that  the  results  of  [S]  do  not  take 
into  account  uncertainties  in  the  compensator  model  and  that 
uniform  convergence  on  compact  sets  in  <g  is  difficult  to  relate  to 
time-domain  performance  of  the  system. 

Another  notable  asymptotic  robustness  theory  is  that  of  [6] 
where  the  graph  topology  is  introduced.  Let  <P  and  Q.  each  be  the 
space  of  all  rational  matrices,  3  the  space  of  strictly  proper  and 
stable  rational  matrices,  equipped  with  the  Hm  norm,  and  (R  = 
e~'  (3).  The  graph  topology  is  the  weakest  topology  on  (P  and  Q, 
under  which  6  is  continuous.  We  have  shown  in  [9],  however, 
that  singularly  perturbed  systems  generically  do  not  converge  in 
the  graph  topology;  hence,  in  this  sense,  robust  compensation  in 
the  presence  of  order  uncertainty  is  unattainable. 

In  view  of  the  shortcomings  of  the  existing  asymptotic 
techniques,  we  wish  to  propose  a  framework  as  well  as  some 
preliminary  results  for  an  alternative  robustness  theory  which  will 
be  taken  into  account:  l)  multirate  and  other  relatively  unexplored 
classes  of  singular  perturbations;  2)  the  necessity  of  dealing 
simultaneously  with  a  large  class  of  system  perturbations,  each 
corresponding  to  a  possible  higher  order  model;  and  3)  time- 
domain  behavior  of  the  closed-loop  system.  Although  treatment  of 
1)  and  2)  seems  on  the  surface  to  be  a  formidable  task,  we  will  see 
that  it  is  possible  to  approach  the  problem  in  a  roundabout  way, 
thus  avoiding  having  to  explicity  characterize  all  possible  parasitic 
phenomena.  We  feel  that  the  inclusion  of  3)  is  a  desirable  feature 
for  any  good  robustness  theory,  since  the  goal  of  system  design 
must  ultimately  be  satisfactory  closed-loop  time  response.  In  view 
of  this  fact,  a  time-domain  approach  has  certain  advantages  over 
frequency  domain  techniques,  since  the  relationship  between  time 
response  and  frequency-domain  behavior  can  be  rather  complex. 

Before  becoming  too  engrossed  in  technicalities,  we  will  briefly 
describe  (in  rough  terms)  the  problem  we  wish  to  address. 
Consider  the  system 

Ex=Ax+Bu 
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where  £,  A,  B,  are  real  matrices  with  E  and  A  square.  We 
assume  that  (4)  exhibits  existence  and  uniqueness  of  solutions  for 
each  initial  condition  x0  and  each  input  function  u;  from  [17]  we 
know  that  this  is  equivalent  to 

|s£-/4|  *  0. 

Such  systems  have  been  studied  extensively  (e.g.,  see  [14]-|  I6|), 
and  are  referred  to  as  singular  when  £  is  singular  and  regular 
otherwise.  The  polynomial 

A(s)  =  |s£-/t|  (5) 

may  be  considered  the  characteristic  polynomial  of  (4)  and  its 
roots  the  eigenvalues  of  £.  An  important  property  of  singular 
systems  is  that  smull  perturbations  in  the  entries  of  £  and  A  can 
change  the  system  order;  one  example  of  this  phenomenon  is  (1). 

Suppose  we  wish  to  find  a  compensator  of  the  form  u  =  Ky  - 
a  which  is  robust  with  respect  to  perturbations  in  £,  A,  B  and  C. 
Since  we  are  inevitably  interested  in  time  response,  we  might  ask 
which  compensators  result  in  a  closed-loop  system  whose  time 
response  varies  continuously  with  £,  A,  B,  and  C,  regardless  of 
the  perturbation.  Unfortunately,  it  is  easy  to  show  that  for  any  K 
there  exist  perturbations  in  the  system  matrices  that  yield 
divergent  behavior  in  the  closed-loop  system  trajectories  for  some 
initial  conditions.  A  more  meaningful  problem  can  be  formulated 
by  first  observing  that  not  necessarily  all  perturbations  in  the 
matrix  entries  of  (4)  are  physically  realistic.  For  example,  a 
simple  RC  circuit  consisting  of  a  single  resistor,  capacitor,  and 
voltage  source  muy  be  modeled  as 

-x+u 

y=x  (6) 

where  x  is  the  capacitor  voltage,  R  »  I,  and  C  =  e.  Positive « 
makes  perfect  physical  sense,  and  it  seems  reasonable  to  try  to 
design  a  compensator  based  on  the  low-order  model  correspond¬ 
ing  to  t  =  0.  On  the  other  hand,  if  e  is  negative,  the  system 
engineer  could  not  expect  to  produce  a  robust  compensator 
without  first  being  aware  of  the  negative  capacitance  and  then 
using  an  appropriate  higher  order  (in  this  case,  first-order)  model. 

A  simple  way  to  characterize  physically  meaningful  perturba¬ 
tions  in  the  plant  is  to  look  at  their  effect  on  plant  trajectories  for 
various  inputs  and  initial  conditions.  For  example,  in  (6)  an  initial 
condition  x0  =  1  yields  x(t)  =  e~,/l  which  converges  on  compact 
subintervals  of  (0,  oo)  as  t  -»  0+ ,  but  diverges  as  t  -*  0" .  Strictly 
speaking,  we  are  really  not  saying  as  much  about  perturbations  of 
(4)  which  can  occur  in  the  physical  world  as  we  are  about  those 
perturbations  which  are  consistent  with  the  measurements  taken 
while  formulating  our  plant  model;  a  system  model  is  good  only  if 
it  is  capable  of  predicting  the  behavior  of  the  actual  physical 
system. 

We  may  now  state  our  definition  of  asymptotic  robustness  more 
precisely.  For  a  given  plant  of  the  form  (4),  a  compensator  is 
robust  if  all  perturbations  in  both  the  plant  and  compensator, 
which  bring  about  only  small  variations  in  the  trajectories  of  each 
system  individually  under  alt  inputs  and  initial  conditions,  result 
in  only  small  variations  in  the  closed-loop  system  trajectories.  The 
meaning  of  the  phrase  small  variations  will  be  precisely  defined  in 
Section  III.  In  the  same  section  we  will  see  that  our  approach 
implicitly  incorporates  the  idea  that  small  system  variations 
should  correspond  to  only  small  changes  in  system  parameters. 

II.  Preliminaries 

i 

In  this  section  we  summarize  the  constructions  of  [10],  [II], 
[I3|,  and  1 14]  which  urc  pertinent  to  subsequent  developments. 

Let 


V- 


y=Cx 


(4) 


l(n,  m,p)»{(E,A,B,  C )  €  rf  0} 
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and  let  £(rt,  m,  p)  be  the  corresponding  quotient  manifold  (see 
1 18])  determined  by  the  equivalence 


(£|>  At,  B i,  C,)*(£2l  /tj,  Bj,  Ci)  iff  C|  =C2  and 
3  nonsingular  M  s.t.  A/£t  =  £j,  MAi  =  Ai,  and  A/Bi  =  £2.  (7) 

(The  arguments  /t,  m,  and  p  will  be  dropped  when  clear  from 
context.)  We  choose  the  equivalence  relation  (7)  because  pre- 
multiplication  by  M  has  no  significant  effect  on  the  system 
representation.  Indeed,  premultiplication  by  M  merely  performs 
elementary  row  operations  on  the  system  of  scalar  equations  (7). 
Hence,  we  are  merely  identifying  systems  formed  front  each  other 
by  reshuffling  the  equations.  We  do  not  wish  to  identify  systems 
which  are  related  by  a  coordinate  change  on  the  state  variable  x, 
since  this  would  reduce  the  system  space  to  one  consisting  of 
input-output  descriptions.  Our  intention  is  to  produce  results 
which  exploit  internal  information. 

The  equivalence  class  containing  <s  =  (£,  A,  B,  C)  is  denoted  { 
=  [£,  A,  B,  C).  In  this  case,  we  say  a  represents  {.  Let 

r=ord  cr  =  ord  £  =deg  A 

where  A  is  the  characteristic  polynomial  (5)  of  {,  and  note  that  a 
unique  matrix  C  is  determined  by  each  £  6  £.  A  sequence  £*  6 
£  converges  weakly  to  $€£({*  ^  £)  if  {  in  manifold 
topology.  Since  £  is  a  quotient  manifold,  the  natural  projection 
(£,  A,  B,  C )  [£,  A,  B,  C]  is  continuous  with  respect  to  weak 
convergence.  Conversely,  we  have  shown  in  [10]  that,  for  each 
convergent  sequence  £*  *+  £  in  £,  there  exists  a  sequence  (£*,  Ak< 
Bk,  Ck)  -*  (E,A,  B,C)  €  S  such  that  [£,  A,  B,  C]  =  £  and  [£*, 
Ak,  Bk,  C*]  =  £*  for  every  k. 

Let  £*  ^  £  =  [£,  A ,  B,  C ]  with  £  singular.  In  { 1 1 )  is  shown 
that  there  exist  nonsingular  matrix  sequences  A/*  -*  M  and  Nk  -* 
N  such  that 

[o  £].w.-[t  ,°,]  m 

where  r  =  ord  |,  Ajk  -*  A„  and  A#  -*  /4/with /1/nilpotent.  For 
sufficiently  large  k,  the  matrices  A/k  and  Atk  are  unique  up  to  a 
similarity  transformation.  For  a  constant  sequence,  the  decompo¬ 
sition  (8)  reduces  to  the  Weicrstrass  decomposition  for  mulrix 
pencils  (see  [17]) 

men= D  a]  •  M“= [ « •  <’» 

The  matrices  M  and  N  may  also  be  used  to  decompose  (4), 
yielding 

AfB=  ,  CN=[C,  C/J. 

Referring  to  [14],  we  say  that  (4)  is  slow  controllable  if  and 
only  if 

rank  [X£-/l  £]  =  «  (10) 

for  every  X  €  ©  and  fast  controllable  if  and  only  if 

rank  [£  B]  =  n.  (11) 

The  system  is  controllable  if  and  only  if  both  (10)  and  (1 1)  hold. 
In  addition,  we  say  that  (4)  is  impulse  controllable  if  and  only  if 

Im  A/+Ker  A/-¥lm  B/=W'~'.  (12) 

(All  four  system  properties  can  also  be  defined  directly  in  terms  of 
the  solutions  of  the  differential  equation  (4),  but  we  find  the  linear 
algebraic  characterizations  more  useful  in  the  context  of  this 
paper.)  Controllability  and  observability  imply  impulse  controlla¬ 


bility  and  impulse  observability,  respectively.  The  corresponding 
definitions  for  observability  are  dual  to  (10),  (II),  and  (12)  (see 
[1.4|).  Since  each  of  these  definitions  is  invariant  under  the 
equivalence  transformation  (7),  we  may  also  consider  the  subsets 
£jt*t  £/o  £c  £/r.  £».  ^/o*  £0.  £f0  C  £  determined  by  (10), 
(11),  (12),  and  their  duals,  as  well  as  the  controllable  and 
observable  systems  £,„  =  £,.  H  £„.  Various  properties  of  these 
spaces  are  studied  in  1 1 3|;  for  example,  £/r  and  £/„  urc  open,  and 
£,f  and  £,-„  arc  dense  in  £. 

Other  important  subsets  of  £  are  the  singular  subspace  £*, 
consisting  of  all  points  [£,  A,  B,  C|  with  £  singular,  the  regular 
subspacc  £"*»£-  £*,  and  the  subspuee  of  unit  index  systems 

£i  =  {£  €  £  |deg  |s£-/$|  =  rank£]. 


In  [10]  it  is  shown  that  £"  is  open  and  dense  in  £;  from  [13],  £c0 
O  £*  is  dense  in  £*. 

Let  D  be  the  set  of  all  C“  functions  Ifi  -♦  Ift  with  compact 
support  and  let  33+  be  the  space  of  distributions  with  support  in 
[0,  oo)  (see  [19]).  To  define  convergence  in  3D+,  we  adopt  the 
weak*  topology:  A  sequence/*  6  S,  is  said  to  converge  to /if 
(/,  4>)  for  every  €  30,  where  (/*,  <t>)  denotes  the  functional/* 
evaluated  at  the  point  d>. 

Associated  with  each  initial  condition  x0  €  IR"  and  each 
piecewise  continuous  input  u  there  exists  a  unique  solution 
(£*)  6  3D"  of  the  system  £*  (see  [17]).  From  linearity  it  follows 
that  the  solution  can  be  decomposed  into  natural  and  forced 
response 

'*mttk)  =  *x0o(tk)  +  *ouVik)- 

Letting 

[jj]  =MkBk,  [Q*  Cy*]  =  C*N*.  [jjjl  =N;'x r0 
we  have  from  [17]  that 

%„(£*)  =  exp  (/t,*)xft+cxp  (/!,*)  *  Bsku  (14) 
where  exp  {A)  G  D"+  2  is  defined  by 

exp  [A)(t)=e,A 


and  denotes  convolution.  Each  Sk*  satisfies  several 
properties  of  continuity.  Indeed,  convergence  of  A,k  guarantees 
uniform  convergence  of  exp  ( Ask )  on  compact  intervals  and, 
hence,  weak*  convergence.  Continuity  of  convolution  with 
respect  to  both  types  of  convergence  assures  that  each  sequence 
VsXQiMk)  converges  weak*  and  uniformly  on  compact  intervals 
whenever  £*  ^  £.  Furthermore,  since  (£)  *  ♦*„„(£)  for  £ 
e  £",  yxou  satisfies  the  same  properties  when  restricted  to  £". 

To  aid  in  writing  a  general  expression  for  ¥{„„(£*),  we  note 
that  there  exists  a  nonsingular  matrix  sequence  7*  (not  necessarily 
convergent)  such  that 


(15) 


r* 1  A*T‘-  [t  l\ 

where  Ajk  is  nonsingular  and  A/k  is  nilpotent.  Then  from  [17], 
exp  (Ajk)  0 


♦yw-r* 


vk-\ 

s 


/•  i 


Tk'X o/k 


+  r* 


exp  (Ajk')*Ajk'BJku 

»k- 1 

-  2  APBAU‘ 


/•o 


(16) 


COBB:  ROBUST  COMPENSATE, 


where 


»/* 


and  61  and  u'  denote  the  /th  distribution  derivatives. 


III.  Problem  Formulation 

We  are  now  in  a  position  to  precisely  state  the  basic  problem 
under  consideration.  For  two  reasons  we  are  forced  to  select  a 
rather  abstract  mathematical  framework  for  our  constructions. 
First,  since  perturbations  leading  to  changes  in  order  require  the 
use  of  systems  of  the  form  (4)  and  since  such  systems  can  have 
impulsive  solutions,  the  space  5D+  of  distributions  and  its 
associated  weak’"  topology  underlie  all  analyses.  Second,  it  will  be 
seen  that  in  order  to  meaningfully  incorporate  the  idea  that  small 
system  perturbations  should  lead  to  only  small  changes  in  the 
entries  of  the  matrices  E,  A ,  B,  and  C,  it  is  necessary  to  identify 
systems  according  to  the  equivalence  relation  (7).  Hence,  we  must 
work  with  the  non-Euclidean  system  spaces  £(n,  m,  p). 

We  consider  the  problem  of  compensating  the  piant  model  (4) 
with  a  static  system  of  the  form 

u  =  Ky  +  v  (17) 

where  K  is  a  matrix  and  u  is  an  external  input.  Let  (P  -  £(«,  m, 
p )  and  Q,  =  ifi"lp,  and  note  that  the  closed-loop  system  takes  the 
form 

Ex=(A+BKC)x+  Bu 

ett.KY-  (18) 


y~Cx. 


In  general,  the  system  (18)  may  not  exhibit  existence  and 
uniqueness  of  solutions  or  may  respond  to  certain  initial  condi¬ 
tions  with  impulsive  transients  (see  [15],  (17);.  Since  we  are  only 
interested  in  choosing  a  compensator  such  that  the  resulting 
closed-loop  system  does  not  suffer  from  either  of  these  defects, 
we  restrict  attention  to 

(R  =  {(^,  /O  6  (PxQ.  |deg  |j£-(/l  +  BKC)\  =  rank  E).  (19) 

Adopting  (19)  is  equivalent  to  assuming  that  (18)  has  unit  index; 
hence,  we  may  set  3  =  £t(n,  m,  p).  Note  that  the  loop-closing 
map  C  is  continuous  with  respect  to  manifold  topology  on  £ 
(weak  convergence);  i.e.,  small  changes  in  the  entries  of  E,  A,  B, 
C,  and  K  bring  about  only  small  changes  in  the  closed-loop 
system  matrices. 

We  say  that  a  sequence  £*  in  (P  converges  weakly  to  £  G  (P 
(I*-4  £)  if  *x0u(Zk)  -  *xou(()  weak*  for  every  x0  and  u  and  if  Ck 
-*  C.  On  the  other  hand,  we  say  that  a  sequence  £*  €  3 
converges  strongly  in  3  if  each  '?,„«({*)  converges  uniformly  on 
compact  subintervals  of  (0,  oo)  and  C  -*  C*.  Uniform  conver¬ 
gence  of  solutions  is  meaningful  for  systems  in  3  =  £,  only 
because  unit  index  systems  have  no  impulsive  components  in  their 
solutions  (see  [17]).  We  have  shown  in  [10]  that  strong  conver¬ 
gence  in  (P  implies  weak  convergence  in  <P. 1  It  is  easy  to  verify 
that  strong  convergence  in  3  implies  convergence  of  each 
'if i0 „(£*)  in  the  weak*  sense  (see  [19]);  hence,  strong  convergence 
implies  weak  convergence  in  3  as  well. 

Although  strong  convergence  of  £*  does  not  necessarily  imply 
that  the  entries  of  the  system  matrices  Ek,  Ak,  Bk,  and  C* 
converge  regardless  of  the  representation  (4)  of  £*,  such  a  strict 


1  Actually,  it  is  shown  in  |  IO|  that  convergence  of  for  every  x„,  u 

guarantees  convergence  of  f,  in  manifold  topology  when  u  ranges  over  2D”.  It 
is  easy  to  show,  however,  that  the  same  result  holds  when  u  is  restricted  to  be 
piecewise  continuous. 


A 
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requirement  would  not  be  particularly  meaningful,  since  premulti¬ 
plication  of  (4)  by  any  nonsingular  matrix  produces  an  entirely 
equivalent  representation.  From  [10]  it  does  follow  that  strong 
convergence  of  £*  implies  convergence  of’  some  representative 
sequence  (Ek,  Aki  Bkt  Ck)  €  2(/t,  m,  p).  We  are  therefore 
justified  in  interpreting  weak  convergence  in  <R  and  3  as 
convergence  of  system  parameters  and  stating  that  t  perturbation 
of  a  system  which  yields  only  small  changes  in  system  trajectories 
also  results  in  only  small  variations  in  system  parameters; 

A  plant-compensator  pair  (£,  K )  €  (R  is  asymptotically  robust 
(or  A  is  a  robust  compensator  for  £)  if  6(£*,  Kk)  **  6(1’,  K )  for 
every  .Vp,  u  whenever  £*  **  £  and  Kk  -*  K.  This  is  equivalent  to 
continuity  of  the  loop-closing  map  6  at  (£,  K)  with  respect  to 
strong  convergence  in  <31  and  3.  It  is  routine  to  verify  that  our 
definition  of  robustness  can  be  couched  in  term*  of  topologies  on 
(P,  Q,,  and  3  simply  by  imposing  on  each  set  the  weakest  topology 
that  makes  each  map  continuous  (see  [20]).  Our  main 
problem  of  interest  is  to  characterize  the  class  of  all  robust  plant- 
compensator  pairs  of  (£,  K )  €  <R  for  any  given  values  of  n,  m, 
and  p.  Equivalently,  we  seek  to  describe  the  class  of  all 
compensators  K  which  are  robust  with  respect  to  a  given  plant 
model  £. 

IV.  The  Cuss  of  Asymptotically  Robust  Compensators 

We  begin  by  presenting  a  result  which  formalizes  the  intuitive 
idea  that  robustness  can  fail  to  hold  only  when  the  plant  model  (4) 
is  singular.  Note  that  when  (4)  is  regular  (£  6  £"),  (R  =  I Rmp. 

Proposition  3.1:  If  £  €  £",  every  K  €  lflmp  is  robust. 

Proof:  Choose  K,  x0,  u,  and  sequences  Kk  -»  K  and  £**♦  £. 
Then  Ck  C.  Since  £"  is  open  in  manifold  topology,  £*  6  £" 
for  sufficiently  large  k.  From  [  10],  £*  £  so  continuity  of  6  with 

respect  to  weak  convergence  implies  6(£*,  Kk)  *+  6(£,  K).  Since 
each  4\0„  is  continuous  on  £"  and  6(£"  x  !Wmp)  C  £",  we  have 

*x0v(6(£*,  Kk))-+*xou( 6(£,  K)). 


Before  starting  our  main  result  on  robustness,  we  need  to 
consider  one  more  algebraic  system  property  of  (4).  We  say  that  a 
system  (4)  is  fast  cyclic  if,  in  the  Weierstrass  decomposition  (8), 
the  nilpotent  matrix  Af  is  cyclic.  If  A/  is  in  Jordan  form,  fast 
cyclicity  is  equivalent  to 


0  1 


:j 

Hence,  from  (9),  a  system  (4)  is  fast  cyclic  if  and  only  if  rank  E  = 
n  -  1.  Note  that  fast  cyclicity  is  independent  of  the  choice  of 
representation  for  £. 

In  order  to  prove  that  a  certain  algebraic  condition  on  the 
compensator  K  is  well  defined,  we  next  present  a  pair  of  lemmas. 
It  will  eventually  be  proven  that  this  condition  is  necessary  and 
sufficient  for  robustness. 

Lemma  4.1:  Let  N  and  T  be  any  n  x  n  matrices  with  N 
nilpotent  and  having  index  q.  Then  Ker  N  is  N  ^T-in variant. 

Proof:  Since  A'(<V‘?_I7')  =  NqT  =  0,  Im  Nq~'T  C  Ker 
N.  Hence, 


(N9~lT)  Ker  Af  c  KerN, 


□ 


Note  that,  if  N  is  cyclic,  Ker  N  is  one-dimensional. 

Lemma  4.2:  If  £  6  £J  is  fast  cyclic,  impulse  controllable,  and 
impulse  observable,  and  ( A /,  B{,  C/)  is  obtained  from  the 
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Weierstrass  decomposition  of  any  representation  of  £,  then 

Ay'-'U'.'+B/KC,)  IKcr/tyX)  (20) 

determines  a  nonempty  open  affine  half-space  in  'St""1  which  is 
independent  of  the  representation.  (The  vertical  bar  denotes  the 
restriction  of  the  linear  operator  to  the  subspace  Ker  Af.) 


From  [12,  Lemma  4.3],  pj  can  be  factored  as 

Pj(s) m +  b„. ,.ys- '■  +  •••  +  boj)  n  (*#»“  »)  (25) 

i-i 

where  each  cr,,  -»  0  and  bn’,  i  =  0,  ••*,/*  -  1  all  converge. 

n _ .  _ rr..  L  _ I  » _ J  .a 


restriction  or  me  linear  operator  to  tne  suospace  Ner  /!/,;  v  *7*  uyt  ■  -  ♦  t*  *  .  ^ 

Proo/-  In  view  of  Lemma  4. 1 ,  (20)  is  well-defined.  For  the  Equating  the  coefficients  in  (24)  and  (25)  of  s'  and  s»  yields 
case  r  *  /»  -  1,  /ly  =  0  and 


zi;"-|(/+S//fC/)=l  +  a/A-C/.  (21) 

For  r  <  n  -  1,  choose  a  nonsingular  T  so  that  T~[AjT  is  in 
Jordan  form.  Letting 


-r-'ty,  tc,  •••  cH.r]  =  C,T 


we  know  from  [14]  that  impulse  controllability  and  impulse 
observability  guarantee  b„ *  0  and  ct  *  0.  Also, 


j,  b'-rKc,  ■•*  b„.rKcn-r. |  1  +  bn-,Kcn.r  ^ 

. oj  1 


Hence, 


Anrt~'(ln-r+BfKCf)  |Ker  A/=b„.rKc\ .  (23)  Our  main  result,  Theorem  4.4  ct 

robust  static  compensator  gains  K. 

Setting  (21)  and  (23)  positive  determines  nonempty  open  affine  Theorem  4.4:  Let  £  G  £>’. 
half-spaces.  1)  A-  robust  K  6  (Rpm  exists  ifl 

From  [11],  (A/,  fly,  Cf)  is  unique  up  to  similarity  transforma-  controllable,  and  impulse  observable 


«„*  =  (-  l)'""  lim  <pj. 

For  sufficiently  large  j  it  follows  that 

sgn  fl  a,j  =  sgn  ctj  sgn  <pj 

/« i 

=  -  sgn  sgn  lim  y>, 

s(_l)«-X+|. 

Hence,  for  each  sufficiently  large  j  there  must  exist  an  /  such  that 
ay  G  08,  <T/j>0. 

Since  o,j  0,  there  exists  j  >  k  such  that  l/ay  >  k.  Set  X*  =  1/ 
ay  and  ek  -  <xj.  □ 

Our  main  result,  Theorem  4.4  completely  characterizes  the 
robust  static  compensator  gains  K. 

Theorem  4.4:  Let  £  G  £s. 

1)  A  robust  K  G  llpm  exists  iff  £  is  fast  cyclic,  impulse 


tion  for  a  given  £.  Clearly,  similarity  transformation  does  not  alter 
(21),  so  the  resulting  half-space  is  unchanged.  To  see  how  (23)  is 
affected  by  similarity  transformation,  note  that  (23)  means 

A}-'- '(/„.,+ BjKC,)x= b„„Kctx 
for  any  x  G  Ker  Aj.  Let  z  -  T~'x.  Then 

T- 'A}-'~' '(/„-,+ BfKCf)  Tz Ac,  z 


2)  Under  the  conditions  of  part  1),  K  is  robust  iff 


A'-'-'V+BjKC/)  |Kcr  A/>0. 


Proof: 


i)  (Necessary):  Let  r  =  ord  £.  We  need  only  consider  the  case 
r  <  n  -  1,  since  r  =  n  -  I  implies  £  is  fast  cyclic,  impulse 
controllable,  and  impulse  observable  (see  [14]),  Suppose  r  <n  - 
1  and  choose  a  representation  (£,  A,  A,  C)  for  £.  Invoke  the 
Weierstrass  decomposition  (9),  select  a  similarity  transformation 
to  put  Aj  in  Jordan  form,  use  the  notation  (22),  and  let 


(T-'AfTr-'-'Vn-'HT-'Bj)  0  yi 

*  K(CfT))\{T-'  Ker  A,)-bm-rKcx.  mT~'A,T.  (27) 

But  T*1  Ker  Aj  -  Ker  (T~'A  T)  so  the  resulting  half-space  is  •  •  yn_r.x 

again  unchanged.  □ 

A  final  technical  lemma  is  needed  to  prove  our  main  robustness  0 

theorem. 

Lemma 4.3:  Let  a,*;  /  *0,  p  be  convergent  sequences  in  (Each  yt  is  either  0  or  1.)  If  £  is  not  fast  cyclic,  impulse 

88  with  a  a,*  #  0  for  every  k,  and  let/**:  88  -*  (8;  i  =  1 ,  •  •  • ,  v  -  controllable,  or  impulse  observable,  then  either  yi  =  0  for  some  /, 

1 ;  k  -  1 ,  2,  •  •  •  be  continuous  at  the  origin  and  satisfy  /*( 0)  *  b„.r  =  0,  or  ct  =  0  (see  [  14]).  Choose  nonzero  sequences 

0,  where  v  >  p.  Then  there  exists  a  sequence  in  IK  such  that  for  yh  b„.r,k  -*  b„.„  C\k  -*  cu  and  Kk~*  K  such  that 

each  k\ 

1)  0  <  |e* |  <  l/k  n-r-i 

2)  sgn«*=  -sgn  a,*  ,  w_.  ^  b„.,.kKk ctk  U  ylk<0  (28) 

3)  the  polynomial  eks"  +  f,.,,k(ek)s’ +  •••  +  f., 

fn*l.kUlt)Sl>*i  +  +  f»k((k))SM  +  •••  +  (flu  +  /u(«*))$  + 

flo*  has  at  least  one  real  root  X*  with  X*  >  k.  for  every  k,  and  define 

Proof:  Fix  k,  let  <xj  =  -  l/j  sgn  fl„t,  and  consider  the 
sequence  (in  j)  |  a,  | 


bn-r,kKkC\k  jj  yik< 0 


for  every  k,  and  define 


Pj(s)=*ajs’+f,.  i  ,k(aj)s’- 1  +  •  •  • 

'  +  Kk+f,k{aj))s*  + 
+  (Oik+f\k(.<Xj))S+a<)k. 


Btk-T 


bn-r-l 

bn-r.k 


C»  =  [CuC2  *"  Cn-r\T~ 
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A/k(x)™T 


.  Ifn-r-Uk 


Now,  we  may  uniquely  define  sequences  aik\  i  ~  0,  •••,/»  -  1 
and  polynomials  plk',  i  -  I.  •••»  n  -  1,  with  Pik(0)  =  0, 
according  to 

**-'**+ (<*„_!  .*+/>„-!  ,*(*))*"-'+ •••+(fli*+Pi*(*))*+ffo* 

_  st-iA.  +  B^C,)  -B.KtCfK 

-BftKkC,  sAjk(x)-U+B/kKkC/k) 

By  elementary  matrix  arguments, 

rt-r-i 

r.*AifcC|*  n  7*.  (31) 

/-I 

Letting /»(*)  *  p(*( i/  =  n,  and  n  =  n  -  1,  we  may 
select  a  sequence  e*  satisfying  the  properties  in  Lemma  4.3  and 
define 


Ek=M 


M-'[o 

Bk*=M~l  ,  C*  =  [C,  Cjk\N- 


-',A„=A 


we  have 


♦  y*+r?it)s  f, 


_  v  _  z -  i  «</* 

~^nHds‘s  ,--i 
»«0 


0  .  71* 


7»-/- 1.*  1 


Consider  the  matrix 


Since  tk  -»  0  and  sgn  e*  =  -sgn  (28)  and  (31)  guarantee 
that  a*  -*  +oo.  If  we  set 


1  * 
a* 

S*  =  *.  '1 


A  routine  calculation  shows  that  the  (/,  y  )th  entry  of  SA_I  exp 
(/  S;')  is 

_ _  '*'• 

It  was  shown  in  the  proof  of  {10,  Theorem  4)  that 


det  (sEk -  (Ak  +  BkKkCk)) 

-  0(«*s',+ (a„.  ,.*+/,.  t,kUk))s"- 1  +  •  •  • 

+  (0i*+/i(e*))s+a<>*)  (33) 

for  some  constant  0.  From  Lemma  4.3,  (33)  has  at  ieast  one  real 
root  X*  >  k  for  every  k.  Since  (33)  is  just  the  characteristic 
polynomial  of  the  closed-loop  system  6({,  K ),  (15)  shows  that  X* 
must  be  an  eigenvalue  of  the  closed-loop  A^’for  sufficiently  large 
k.  Thus,  exp  (Ajkl)  cannot  converge  uniformly  on  compact 
subintcrvals  of  (0,  «),  since  this  would  imply  uniform  conver¬ 
gence  of  its  eigenvalue  exp  (X*).  Letting  u  =  0,  it  follows  from 
(13)  and  (16)  that  ♦*p*(6(?*»  K ))  does  not  converge  for  every  Xo. 

In  order  to  prove  that  K  is  not  robust,  we  have  only  left  to  show 
that  ♦,.,($*)  -»  *,„,({)  in  the  weak*  sense,  where  =  [£*,  Ak, 
Bk,  Ck  j.  To  do  so,  we  note  that 

*^-«p  (a.  (-j)'1)  *+*»  (-5)” 


'I  £ 


S;'exp(S;')-»- 


Therefore,  tfy*  -*  -6J~‘  for  each  j  &  1,  and  ^fmk~*  for 
any  xnu.  It  follows  from  (13),  weak*  continuity  of  convolution, 
and  continuity  of  Sf‘IQU  that  *,„„({*)  -*  ♦,„*($). 

(Sufficient):  Let  f*  be  any  sequence  in  £  such  that  ♦««(?*)  -♦ 
¥ ,„„(£)  for  every  x0u  and  with  Ck  -»  C.  Then,  from  f  10],  {*  -*  £ 
and  the  decomposition  (8)  may  be  invoked.  It  follows  that 

\Mk\  \sEk-Ak\  \Nk\.\sfr-A*\  |*4>»-/wl 

for  any  convergent  representation  (Ek,  Ak,  Bk,  Ck).  Suppose 
(-  \)n~'\Afk\  <  0  for  infinitely  many  k.  Since  \Ajk\  is  just  the 
product  of  the  eigenvalues  of  A/k,  there  exists  a  subsequence  of 
A[kj  with  at  least  one  real,  positive  eigenvalue  for  each  kj.  It 
follows  that  Aj)  has  an  eigenvalue  X;  -»  ».  Let  u  ■  0  and 
observe  that  ‘ 


,(?*/)»  $)a?kjX0/kj 
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!*,«{ *  me{tA ,1)  dt. 


Applying  the  decomposition  (8)  to  the  closed-loop  system  yields 
nnnsingulur  transformations  Mk  ami  Nk  such  that 


A  function  <f>  6  SD  can  always  be  chosen  such  that  an  eigenvalue 
of  F kj  satisfies 

1  exp  (>  'ty(f)  dt-*<».  where  all  sei 

J°  be  applied  tc 

Hence,  || Fkj )|  -♦  ».  It  follows  that  r*.*b/*  is  unbounded  for 
an  appropriate  choice  of  Xq  and  that  ^^(f*)  and  [from  (13)] 

Vx0u(kk)  are  not  convergent.  This  contradiction  leads  us  to 
conclude  that  (- >  0  for  sufficiently  large  k.  From  Lemm 

Appealing  to  the  notation  of  (22)  and  (27),  we  have  b„.„  ct  it 
0  and,  if  r  <  n  -  1,  =  •  •  •  =  y»-r- 1  =  1.  Choose  K  to 

satisfy  the  condition  (26).  For  r  =  n  -  1,  (21)  indicates  that  1  + 
b\Kc\  >  Oj  for  r  <  n  -  1,  (23)  implies  that  b„-fKcs  >  0.  ^  = 

Defining 


_  k-  0 

~  °  °k  J  * 


tik{Ak+BkKkCk)Nk-. 


%  k  0 

0  l. 


f 

KJ  ’ 


CkNk  =  \C\k  Cy*] 


we  have 


A*(s)  =  |s£*  -  (Ak  +  BkKkCk)  | 


|M*|A*(s)|N*| 

sir  ~(Ask+ BskKk  Csk)  -B,kKkCjk 

-BfkKkCsk  sA/k-tfn-r+  BfkKkCfk) 

=  My»|s"  +  (f.-i  - \A'fk\tr{A,k+BikKkCsk))s'"'  +  •  •  • 

(3 

where  ak  is  defined  by 

\sAjk- (/„-,+/»y*A*Cy*)|  =  \Ask\s"-' >+•••. 
From  elementary  matrix  arguments  we  have,  for  r  <  n  -  1, 
\sA,-(In.r+B/KC/)\=(-\y-rb„-rKcis"-'-'  +  -’- 


where  all  sequences  converge.  The  decomposition  (13)  rhiy  also 
be  applied  to  the  closed-loop  system  yielding 

.»-/».  [^t!]  • 

From  Lemma  4.3,  ’^(l*)  -*  ijoB($),  and 

(  exp  (  —  )  Ar0y*  +  —  exp  (  —  )  *  5y*«  ifo*<0 

*/  #t.wl  W  °k  W 


-B/kU 

so,  as  in  the  necessity  proof  of  part  1),  $1  u(£*)  -♦  1i '({JX)  for 
any  x0,  u.  Hence,  ^0u(G((k))  -»  %gU(G(b)  and  K  is  robust. 

2)  (Sufficient):  This  part  has  already  been  treated  in  the 
Sufficiency  section  of  1). 

(Necessary):  Invoke  the  Weierstrass  decomposition  (9).  If  (26) 
fails  and  r  *  «  -  1,  we  have  1  +  BfKC )  £  0  so  Bj  *  Oand  Cf 
*■  0;  hence,  there  exists  a  sequence  Kk~*  K  such  that 

\+B,KkCf<0  (39) 

for  every  k.  Now  define  alk,  a„.ly,  0i*,  Bn- i,t 
according  to 

xs"  +  (art.i,*  +  i8n.|,*x)s""l+  •  •  •  +(fl|*  +  0i*x)s+aojk 

sln.{-(As  +  BsKkC,)  -BtKkCf 

-BfKkC,  xs-O+B/KkC, )  ' 


lim  a*=  j  7}\tb'rHCl\r  If  nn~-r>  I  (35)  Letting  fk(x)  =  0,kx,  we  can  find  a  sequence  «*  satisfying  the 
v'  '  1  *  properties  in  Lemma  4.3;  define  a„*  =  -  1/t*.  Since  e*  0  and 

From  our  choice  of  A:  it  follows  that  the  closed-loop  system  C«,  Jn  (k  =  ~  s«n  >•*’ (39)  and  (40)  *uarantee  that  "  +  °°- lf 
K)  exhibits  no  impulsive  behavior  in  its  natural  response,  i.e. ,  56 

ord  6(£,K)  =  rank£  -  n  -  1.  Hence,  from  [12,  Lemma  4.3]  r  *i 

we  know  that  ,  I  /"'1  0  I  .  .  . 

£*  =  A/-1  .  N-',Ak~A 


On-\,k~  “(1  +B/KkC/). 


if  n-r=l 
if  n-r>  1. 


we  set 


I Mk  I  A*(J)|N* |  =  f*(a* s- 1)  II  Is ~  x«> 


0  — 


N-',Ak~A 


Bk  =  B,  Ck-C 


where  <pk,  ak,  and  X*  all  converge  and  lim  ak  =  0.  Matching  we  have  that  det  (sEk  -  (Ak  +  BkKkCk))  has  at  least  one  real 
coefficients  in  (34)  and  (36)  yields  root  X*  >  k  for  each  k.  As  in  the  sufficiency  proof  of  part  1), 

¥,  „(©({*,  Kk))  does  not  converge  for  some  jr0,  «•  Since 

lim  ak=  -lim  *,*.  (38)  K^k)  =  cxf>  <-“*>*■>-“*  CXP  *  B*u  (4,) 

From  (35)  and  our  choice  of  K,  (-l)""r  lim  or*  <  0.  Hence,  in  the  open-loop  system,  ^(f*)  -  **<>«($)  and  K  15  not 
fmm  CTRI  robust. 

ii uni  \Jvft  tr  ^  M  _  l  u.a  Mm,  /lftnni  tkd  nAiafinii  /OQ\  anri 


(- 1)'*"''  lim  v>*  =  (- 1)'*-'*1  lim  a*<0. 
Thus,  (-  1)"-V*  <  0  for  sufficiently  large  k,  and 

•»*  Hisr51- 


robust. 

If  (26)  fails  and  r  <  n  -  I ,  we  may  adopt  the  notation  (29)  and 
(30)  and  observe  that  b„.,Kc,  s  0.  Since  fast  cyclicity,  impulse 
controllability,  and  impulse  observability  guarantee  that  b„.r  *  0 
and  c,  *  0,  a  sequence  Kk~*  K  may  be  chosen  so  that  b„.rKkc{ 
<  0  for  every  K.  The  remaining  arguments  are  the  same  as  in  the 
necessity  part  of  1)  with  b„.r,k  =  b„.„  C|*  *  Cj,  and  yik  =1.0 
Theorem  4.4  is  somewhat  pessimistic  in  that,  in  the  strictest . 
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theoretical  sense,  robustness  can  only  be  guaranteed  when  at  most 
one  degree  of  singularity  is  present  in  the  plant  (4)  (rank  E  =  n  - 

1).  In  physical  terms  this  can  be  interpreted  as  meaning  that  a 
static  compensator  can  handle  only  a  first-order  unmodeled 
dynamic  element.  In  our  opinion,  this  indicates  that  some  basic 
assumptions  which  are  as  yet  not  well  understood  are  convention¬ 
ally  placed  on  system  models  in  engineering  practice. 

For  a  mathematical  explanation  of  how  nonrobust  compensa¬ 
tors  may  fail  to  stabilize  a  system,  consider  the  matrix  condition  in 
part  2)  of  Theorem  4.4.  This  condition  determines  an  open  affine 
half-space  in  the  set  ‘Si""’  of  compensation  gains  K.  Examination 
of  the  proof  of  Theorem  4.4  reveals  that,  for  systems  which  are 
fast  cyclic,  impulse  controllable,  and  impulse  observable,  a  static 
compensator  results  in  positive  feedback  either  for  all  admissible 
perturbations  simultaneously  or  for  none  at  all.  The  half-space  of 
robust  feedback  gains  is  simply  the  set  of  all  K  with  the 
appropriate  sign  to  guarantee  negative  feedback  for  all  perturba¬ 
tions  of  the  system  (4).  The  system  (6)  illustrates  this  point.  The 
robust  gains  are  simply  those  satisfying  K  <  1.  On  the  other 
hand,  part  1)  of  Theorem  4.4  maintains  that  unless  the  plant  is  fast 
cyclic,  impulse  controllable,  and  impulse  observable,  the  class  of 
admissible  perturbations  is  so  broad  that  any  compensator  results 
in  positive  feedback  with  respect  to  some  perturbation;  hence,  no 
compensator  is  robust.  This  is  illustrated  by  (2)  and  (3). 

Another  important  point  to  note  at  this  stage  is  that,  although  alt 
definitions  and  technical  arguments  until  now  have  been  couched 
in  terms  of  sequences,  each  statement  applies  equally  well  to  nets 
in  the  various  topological  spaces.  This  observation  is  important, 
since  the  space  of  distributions  D ,  does  not  satisfy  the  first  axiom 
of  countability  (see  (201), 

To  conclude  this  section  we  compare  our  results  to  those  of  |5j. 
Specifically,  (5,  Theorem  1 1  shows  that,  for  any  system 

x=Ax+Bu 

y=Cx  (42) 

and  any  compensation  matrix  K,  there  exists  a  singular  perturba¬ 
tion  of  (42)  of  the  form  (1)  which  destabilizes  the  closed-loop 
system.  (The  result  of  (5|  is  somewhat  more  general  in  that  it 
applies  to  all  dynamic  compensators  which  are  proper  but  not 
strictly  proper.)  According  to  Theorem  4.4,  if  we  take  such  a 
perturbation  and  set  e  =  0,  we  obtain  a  nominal  system 

[*  xHa  iMJ]  ■ 

y-\C\  C2(x  (43) 

which  must  either  fail  to  be  fast  cyclic,  impulse  controllable,  or 
impuise  observable.  For  example,  setting  t  -  0  in  (3)  yields  a 
system  of  the  form  (43)  which  can  be  shown  to  be  not  fast  cyclic. 
While  the  result  of  (5|  illustrates  that  a  specialized  cluss  of 
parasitics  can  lead  to  closed-loop  destabilization,  our  results 
characterize  the  same  phenomenon  but  in  the  context  of  a  broader 
class  of  perturbations  and  a  larger  family  of  nominal  systems.  For 
example,  our  Theorem  4.4  applies  to  systems  of  the  form  (43) 
with  An  singular  (as  long  as  |s£  -  A  |  #  0  is  satisfied),  while  [5] 
considers  only  the  case  of  An  nonsingular.  Our  result  also  shows 
when  destabilization  can  occur  as  a  result  of  perturbations  to  a 
given  order;  the  perturbed  order  required  to  destabilize  the  closed- 
loop  system  in  |5(  is  not  specified. 

V.  Genericity 

We  now  consider  the  class  of  systems  (4)  for  which  there  exists 
a  robust  compensator  K.  The  sets  of  impulse  controllable  and 
impulse  observable  systems  were  shown  in  1 13|  to  lie  dense  in  the 
system  space  £.  The  next  result  characterizes  those  systems 
which  are  also  fast  cyclic. 


ib?  • 

Proposition: 

1)  £  is  fast  cyclic  iff  {  €  £"  U  £"'1. 

2)  (£"  U  fi"-')  n  £/f  fl  £„  is  open  in  £. 

3)  £*~'  fl  £/(n  £,0  is  dense  in  £*. 

Proof: 

1)  Let  (£,  A,  B,  C)  be  any  representative  of  £.  If  r  ■  ord 
the  Weicrstrass  decomposition  (9)  shows  that  rank  E  *  r  +  rank 
A/.  But  $  is  fast  cyclic  if  and  only  if  either  r  -  not  rank  A;  *  n 
-  r  -  1 .  Hence,  £  is  fast  cyclic  iff  rank  E  *  not  rank  E  -  n  - 
1. 

2)  Let 

e=[£,  a,  b,  c ]  e  n=(£»  u  £"-»)  n  £,c  n  £,0 

and  apply  the  decomposition  (9).  Then  A/  is  cyclic,  and  Ker  A/  C 
Im  Aj.  Since  |  is  impulse  controllable  and  impulse  observable. 

!m  A/+Im  fl/=Im  /t/+Ker/t/+Im 

Ker  Af  ft  Ker  C/=  Ker  Aj  ft  Im  A/  D  Ker  C/=  0 

so  £  €  £/c  ft  £/0  (see  (14J).  It  follows  that 

n=(£B  n  £"-')  n  £/c  u  £/0. 

We  know  from  (13)  that  £/c  and  £/„  are  both  open,  so  0  is  the 
finite  intersection  of  open  sets. 

3)  It  was  shown  in  [13]  that  £„  O  £"_1  is  dense  in  £*.  Our 

result  follows  immediately,  since  £/c  ft  £,„  D  £„.  □ 

Note  that  part  3)  is  stated  in  terms  of  the  singular  subspace  £*. 
Since  every  point  in  the  regular  subspace  £"  is  necessarily  fast 
cyclic,  impulse  controllable,  and  impulse  observable  (see  [14]) 
and  since  £"  is  dense  in  £,  density  of  (£"  U  £"" ')  ft  £(c  fl  £(0 
in  £  is  trivial.  Part  3)  is  a  much  stronger  result. 


VI.  Discussion  and  Conclusions 


In  this  section  we  discuss  some  of  the  implications  of  our  theory 
and  use  these  to  suggest  further  research.  Theorem  4.5  shows  that 
a  generic  class  of  systems  can  be  robustly  compensated  using 
static  compensators  K.  This  does  not  mean,  however,  that  the 
complement  of  the  open  and  dense  subset  (£”  U  £"_1)  D  £tc  f) 
£lo  does  not  contain  interesting  systems.  On  the  contrary,  it  is 
easy  to  show  that  all  systems  of  the  form  (4)  with  r  <  n  -  I  and 
A  n  nonsingular  lie  outside  the  generic  class  described  by 
Theorem  4.5.  Another  interesting  observation  is  that  even  a 
system  which  does  lie  in  the  generic  set  can  be  trivially  augmented 
so  that  it  sits  outside  the  generic  set  in  a  higher  dimensional  system 
space.  For  example,  the  dimension  of  (4)  may  be  increased  simply 
by  defining  a  new  (scalar-valued)  state  variable  z  =  0  and  noting 
that 


[o  d[(Ho  ?][;]♦[:]■ 


y-\c  Q| 


X 

z 


(44) 


System  (44)  is  a  member  of  £(n  +  1,  m,  p).  It  is  easy  to  show 
that  (44)  is  not  fast  cyclic  and,  hence,  cannot  be  robustly 
compensated.  The  latter  point  can  be  countered  by  arguing  that 
only  variables  of  interest  should  be  included  in  a  well-devised 
state-space  model;  therefore,  the  variable  z  would  never  be 
present. 

There  are  at  least  a  couple  of  avenues  of  research  which  might 
eventually  resolve  these  issues.  Dynamic  compensation  is  still 
relatively  unexplored  in  the  context  of  singular  perturbations.  One 
promising  result  is  |5,  Theorem  2|  which  suggests  that,  when 
parasitics  are  present,  strictly  proper  compensators  are  more 
robust  than  nonstrictly  proper  ones.  Since  [5]  treats  only  the  single 
time-scale  case,  more  work  needs  to  be  done  to  see  whether  this 
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result  stands  up  to  a  larger  class  of  perturbations.  As  pointed  out 
in  Section  I,  the  issue  of  which  class  of  perturbations  is 
meaningful  in  a  given  system  anuiysis  is  of  fundamental  impor¬ 
tance.  Our  main  results  can  in  fact  be  proven  under  a  somewhat 
more  general  definition  of  system  perturbation  than  the  one 
provided  here  (strong  convergence).  However,  preliminary  work 
suggests  that  even  such  a  generalization  might  be  too  restrictive  to 
allow  a  coherent  robustness  theory  to  be  developed.  We  intend  to 
explore  these  issues  more  fully  in  the  future. 
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Abstract 

We  consider  the  problem  of  determining  whether  each  point  in  a  polytope 
of  nxn  matrices  is  stable.  Our  approach  is  to  check  stability  of  certain 
faces  of  the  polytope.  For  n£3,  we  show  that  stability  of  each  point  in  every 
(2n-4) -dimensional  face  guarantees  stability  of  the  entire  polytope. 
Furthermore,  we  prove  that,  for  any  kin2,  there  exists  a  k-dimensional 
polytope  containing  a  strictly  unstable  point  and  such  that  all  its 
subpolytopes  of  dimension  min{k-l,2n-5)  are  stable. 

lThis  work  was  supported  in  part  by  NSF  Grants  ECS-8612948  and  ECS-8611728  and 
in  part  by  AFOSR  Grant  No.  88-0087. 
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In  this  paper  we  consider  the  problem  of  ascertaining  whether  certain: 
subsets  of  Rn*n  consist  entirely  of  stable  matrices.  (Here  we  take  stability, 
of  a  matrix  to  mean  that  all  its  eigenvalues  are  in  the  open  left  half  planed)/ 
First  we  need  some  definitions.  A  (convex)  polytope  ?  in  a  vector  space  V  is 
the  convex  hull  of  any  nonempty  finite  subset  of  V.  The  dimension  of  9  is  the 
dimension  of  the  affine  hull  aff(9»)  of  9.  The  relative  boundary  of  9  is  the 
boundary  of  9  as  a  subset  of  the  topological  space  aff(?).  A  face  of  9  is  any 
set  of  the  form  Iln?,  where  n  is  a  supporting  hyperplane  of  9.  A  vertex  of  9 
is  a  O-dimensional  face.  An  edge  of  9  is  a  1-dimensional  face.  A  subpoly.tope 
of  9  is  the  convex  hull  of  any  set  of  vertices  of  9.  Finally,  a  k-dimensiohal 
half -plane  in  V  is  any  nonempty  set  of  the  form  JMtns,  where  R  is  a  closed 
half-space,  S  is  a  k-dimensional  affine  subspace,  and  SfltR.  (Note  that  this 
implies  that  the  affine  hull  of  #  is  simply  S.) 

In  the  robust  control  literature,  considerable  interest  has  been 
generated  by  the  problem  of  determining  whether  stability  of  a  polytope  in 
either  IRn  or  IRnxn  can  be  guaranteed  simply  by  checking  stability  of 
low-dimensional  faces.  (Stability  of  a  vector  x«Rn  means  simply  that  the 
polynomial  so+XnS0-1*. . .+Xt  is  Hurwitz.)  We  first  note  that  the  cases  n=0  and 
n=l  are  trivial;  stability  of  the  vertices  always  guarantees  stability  of  the 
polytope.  Several  recent  papers  consider  the  case  n£2.  For  example., 
polynomial  polytopes  of  a  particularly  simple  structure  ("interval 
polynomials")  were  addressed  by  Kharitonov  [1];  he  showed  that  only  four 
specially  constructed  vertices  need  be  checked.  A  more  recent  result  of 
Bartlett,  Hollot,  and  Lin  [2]  demonstrates  that,  for  an  arbitrary  polynomial 
polytope  checking  all  edges  is  sufficient  to  guarantee  stability  of  9.  With 
respect  to  polytopes  in  R*1*11,  Fu  and  Barmish  [3]  have  shown  that  stability  of 


J 


all  1  -dimensional  subpolytopes  is  insufficient  to  guarantee  stability  of  9. 
DeMarco  [4]  has  shown  that,  for  n£3,  (n-2) -dimensional  faces  are  insufficient, 
but  2n-dimensional  faces  are  sufficient. 

In  this  paper  we  refine  the  bounds  of  [4J  and  show  that  stability  of  all 
m-dimensional  faces  is  sufficient  to  guarantee  stability  of  9,  where 

/  *  fl.  n*2 
m(n)  =  \2n-4 ,  n>2 

Furthermore,  we  show  that  for  any  n  and  kSn*  there  exists  a  polytope  of 
dimension  k,  containing  a  strictly  unstable  point  (a  matrix  with  an  eigenvalue 
X  satisfying  Re^>0),  and  such  that  all  its  min{k-l ,m-l}-dimensional 
subpolytopes  are  stable;  hence,  in  this  sense,  m  is  minimal, 


2.  Sufficiency  of  m 

Throughout  our  analysis,  we  will  make  extensive  use  of  the  fact  that  any 

affine,  one-to-one  map  f:R^  — ♦  Rn*  determines  an  affine  isomorphism  between  R^ 

1c  k 

and  f(R  ).  Among  other  things,  this  implies  that,  for  any  polytope  ?cR  ,  t(9) 

is  also  a  poiytope  of  the  same  dimension  as  9;  furthermore,  f  sets  up  a 

one-to-one  correspondence  between  q-dimensional  faces  of  9  and  q-dimensional 

R 

faces  of  {(9).  In  addition,  f  maps  each  k-dimensional  half-plane  in  R  into 
another  k-dinensional  half-plane  (e.g.,  see  [5]).  Finally,  we  note  that  every 
polytope  is  compact  and  that  the  set  {xeR^J  }  is  a  polytope  whose 

q-dimensional  faces  are  generated  by  fixing  k-q  entries  of  x  at  either  +1  and 
letting  the  remaining  q  entries  vary  independently  over  [ 1 , 1 J . 

With  these  observations  in  mind,  we  prove  a  result  characterizing  the 
affine  structure  of  the  set  of  unstable  points  in  Rn*n. 


-2- 


Lemma  2.1  For  each  unstable  AelR  there  exists  an  luf ml  dimensioha! 

half-plane  #clRn*n  such  that  1)  A e#  and  2)  Bs#  implies  B  is  unstable. 

Proof  Case  I  --  A  has  a  real  eigenvalue  <*650:  Let  T=[v  W] ,  where  v  is  ah 
eigenvector  corresponding  to  /lo  and  W  is  chosen  to,,, make  T  nonsinguMr. 
Clearly,  the  map  f:Rn*  n+1  —*  IRn*n  determined  by 

fu.v.z,»T[j 

is  affine  and  one-to-one.  Let  »  be  the  {n2-n+l)-dimensional  half-plane 

% 

»={f(^,y,Z)|^o.y«IR1><n"1.Z«Rn‘lyn“1} 


Then  A «#  and  every  matrix  in  M  is  unstable.  Since  n2-n+l£n2-a,  we  need  only 
select  any  (n2-m) -dimensional  half-plane  H  satisfying  A«#cjf. 

Case  II  —  A  has  a  complex  eigenvalue  pair  a<j±i0o  with  «o>0;  Let 
T=[v  w  X],  where  v+-iw  is  an  eigenvector  corresponding  to  ao+i/3o  and  X-  is 
chosen  to  make  T  nonsingular.  Let  H  be  the  (n2-2n+-4) -dimensional  half-plane 


tr  Ufc2a0,Y«R2xn_2 


ZeR 


,n-2xn-2 


} 


(trU2t2ao  describes  a  4-dimensional  half-plane,  since  trU*<U,I>.)  H  contains 
only  unstable  points,  since  tr  U£2ag  implies  U  has  at  least  one  eigenvalue  A 
with  Re^>oo-  Also,  As#,  since  our  choice  of  T  guarantees  that  A  has 


L’~ 


Finally,  n2-2n+4>n2-m,  so  the  desired  HcH  exists.  /9 

Next  we  prove  an  easy  result  concerning  the  intersection  of  affine  sets. 

Lemma  2 . 2  Let  V  be  a  p-dimensional  Euclidean  space.  HcV  a  k-dimensional 
half-plane,  and  r  a  q-dimensional  affine  subspace  with  k+q>p.  Consider  any 
vector  x(J«#r>r.  There  exists  a  (k+q-p) -dimensional  half-plane  »  such  that 
x<j«#c#nr. 

Proof  By  definition.  tf=Rns,  where  R  is  a  closed  half-space  and  S  is  a 
k-dimensional  affine  subspace  satisfying  S<R.  There  exists  an  affine  subspace 
Scsnr  with  dim  s"*k+q-p  and  x0«S.  If  ScR,  let  XcRns  be  any  (k+q-p)-dimensional 
half-space  containing  Xo-  Then  JfcRnsnr=#nr.  If  S<R,  let  H=RnS.  Then  xp«#i 
since  x0«»^rcR.  Also,  diatflk+q-p,  since  Jt  is  nonempty.  y 

We  are  now  in  a  position  to  prove  our  first  main  result. 

Theorem  2 . 3  Stability  of  every  matrix  in  every  a-dimensional  face  of  9 
guarantees  stability  of  every  matrix  in  9. 

Proof  Our  arguments  here  are  similar  to  those  used  in  [2, Lemma  1].  Suppose  9% 
is^ir  unstable-  polytope  of  dimension  k>m.  Then  there  exists  an  unstable 
matrix  Ai«?k.  From  Lemma  2.1,  there  is  a  (n2-m)-dimensional  half  plane  #j, 
consisting  entirely  of  unstable  points  and  containing  Aj.  Since  H i  is 
unbounded,  there  exists  an  A2«#j  lying  on  the  boundary  of  ?k  and,  hence,  in 
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one  of  i  *  s  (k  i -dimension.!  faces  Pk-j.  I'l'-im  l.oinnia' 2 . 2 .  MU? 

Xi  'af f  C?k- 1 )  contains  a  (k~m-l)-dimensional  half-plane  it 2  such  that  A  3  •*>&... 
Proceeding  inductively,  we  find  that  there  exists  an  m-dimensional  face  f &  and; 
a  point  Ak-m«?m  .such  that  Ak-m  is  unstable.  .1 


3.  Minimality  of  m 

Our  next  task  is  to  show  that  m  is  the  smallest  integer  such  that 
stability  of  all  m-dimensional  faces  of  9  guarantees  stability  of  9. 

theorem  3 . 1  For  each  integer  n£2  there  exists  an  m-dimensional  polytope  9»clRn*n 
containing  an  unstable  point  and  such  that  all  its  (m-1) -dimensional  faces  are 
stable. 


Proof  Case  I  —  n=*2 :  Consider  the  affine,  one-to-one  map 


and  the  corresponding  i-dimensional  polytope  ?*{f (x) | |xjSl} .  Each  matrix  in  9 
has  characteristic  polynomial  A(s)»s2+s+x2 ;  hence,  each  vertex  of  9  (x*±l)  is 
stable,  but  the  point  corresponding  to  x*0  is  unstable. 

Case  II  —  n*3:  The  2-dimensional  polytope 


[-? 


x 


1 

0 

y 


-x 

-y 

-1 


si) 


has  characteristic  polynomial  A(s)=sJ+s2+(l*x2+y2 )s+l .  Each  coefficient  of 
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Ms)  -,s  t>us!iivc.  and  rhe  corresponding  2  2  Hurwiu:  matrix 

principal  second -order  minor  equal  to  Mg(x, y) =x2*y* .  thus,  each  ->dge  is 

stable,  but  the  matrix  corresponding  to  x=y*0  is  unstable. 

Case  III  —  n£4:  Consider  the  (2n-4) -dimensional  polytope 


0  1 
-1  0 
x  y 


-x 

1 

-y 

-i 


A  routine  calculation  shows  that  9  has  characteristic  polynomial 
p(s)»(s+l)n  4A(s),  where 


40  TT2  TT  TT  T2 

A(s)*s  +2s  +(2+x  x+-y  y)s  *(2+x  x+y  y)s+l+x  xy  y-(x  y) 


From  the  Schwartz  inequality,  it  is  clear  that  all  coefficients  of  ^  are 
strictly  positive.  The  corresponding  4x4  Hurwitz  matrix  has  its  leading 
principal  third-order  minor  equal  to 

M3 ( x , y ) *4xTx+4yTy+4 ( xTy ) 2+ ( xTx-yTy ) 2 

Clearly,  MgSO  with  equality  if  and  only  if  x*y=0.  Thus  the  (2n-5)-dlmensional 
faces  of  9  are  stable,  but  the  point  corresponding  to  x=y=0  is  unstable.  J 

It  is  interesting  to  note  that  Theorem  3.1  also  implies  that  the 
half-planes  considered  in  Lemma  2.1  are  maximal  in  the  sense  that  there  exists 
an  unstable  matrix  A  in  IRnxn  such  that  every  half-plane  of  dimension  greater 
than  n2-m  containing  A  must  also  contain  a  stable  matrix,  Indeed,  if  this 


wore  no t  r.ho  .wsa.  the  arguments  in  Theorem  2.3  could  bo  used  to  prove  unit  m 
is  not  minimal . 


4.  A  Stronger _Ve rsion  of  the  Minimality  Theorem 

The  construction  in  the  proof  of  Theorem  3.1  is  weak  in  three  respects;. 
1)  The  polytope  ?  contains  only  a  single  marginally  unstable  matrix  (i.e.  a 
matrix  having  all  eigenvalues  X  satisfying  Re  -*50  and  at  least  one  with 
Re  >1=0).  2)  The  construction  yields  only  a  polytope  of  dimension  m.  3) 
Arbitrary  subpolytopes  are  not  considered;  thus  it  is  not  clear  that  checking 
all  suopolytopes  of  dimension,  say,  m-1  would  not  guarantee  stability.  The 
minimality  proof  would  be  more  convincing  if  it  could  be  extended  to  give  a 
family  of  polytopes,  each  1)  containing  a  strictly  unstable  point  (and,  hence, 
infinitely  many  unstable  points),  2)  having  arbitrary  dimension  k.  and 
3)  having  all  min{k-l .m-1 }-dimensional  subpolytopes  stable. 

Theorem  4.2  shows  that  such  improvements  over  Theorem  3.1  can  be  made. 
The  proof  requires  a  simple  lemma.  For  any  normed  linear  space  V,  subset  QcV, 
and  point  rev  consider  the  distance  function 

d(7,Q)*inf  ’y-wii 
««0 

Let  conv(O)  denote  the  convex  hull  of  0. 
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J 


Lemma  i  1  Suppose  OrV  is  convex.  3>0.  and  L*"V  is  any  set:  stich  '  ha:,  t s . ■.*,]  ■ 
for  every  feT.  Then  d(j),0)<e  for  every  7)econv(Quf ) . 

Proof  Each  r?econv(Qur)  is  of  the  form  j?=t*7? where  ae(0.1]  and 
r}i,rj2«nur.  There  exist  (ult<u2en  such  that  1  >?  1  -cux :: <«  and  !j7j-»2  <e.  Let 
0*au>i+-(  i-a)<u2  •  Then  <u«0.  and 


d(t7,n)S!ir)-<i»8 

*i|a(r}i-«(l)  +  (l-a)(j)2-0a)  1 
<a||rji-a»ti|  +  (l-a)  Srjj-wail 
<£  2 

Theorem  4.2  For  any  integers  n,k  with  n£2  and  15kSn2,  there  exists  a  polytope 
of  dimension  k  containing  a  strictly  unstable  point  and  such  that  each 
min{k-l ,m-l}-dimensional  subpolytope  is  stable. 

Proof  Suppose  a  marginally  unstable  polytope  of  dimension  k  is  constructed 
such  that  all  its  min{k-l,m-l}  subpolytopes  are  stable.  Then,  since  the  set 
of  stable  points  in  Rnxn  is  open  and  the  union  of  all  min{k~l ,m-l}-dimensional 
subpolytopes  of  9^  is  compact,  the  subpolytopes  of  of  the  same 

dimension  are  stable  for  sufficiently  small  e,  but  is  strictly  unstable. 
Thus,  it  suffices  to  construct  any  k-dimens tonal  unstable  9^  with  stable 
subpolytopes. 
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If  u  -2.  I  ft  f(x.y)  - 


i -x 


x' 

•1, 


where  x.y  range  over  IR;  otherwise,  let. 


f(x.y)= 


0  1 
-1  0 
x  y 


where  x .  y«IR 


n-2 


fk:R 


,nxn 


We  consider  two  cases;  first. 

) ,  where  the 


according  to  f^( 


V 

‘z 

0, 

* 

0, 

partitioned  in  any  way  such  that 
one-to-one,  the  set 


pv. 

z, 


Since 


assume 
vectors 
each  f^ 


k<m.  Define 
x  and  y  are 
is  affine  and 


V{fk(w,z) 


51} 


is  a  k-dimensional  polytope.  As  in  the  proof  of  Theorem  3.1,  each  matrix  in 
flk  is  stable  except  for  the  point  corresponding  to  w=z=0.  The  union  of  the 
(k-l)-dimensional  subpolytopes  of  is  compact  and  nowhere  dense  in  f (IR^ ) ; 
hence,  there  exist  vectors  Wo.Zosf^)  such  that 


SV*fk(w+v,°,z+z°) 


51} 


is  unstable,  but  has  all  its  (k-1) -dimensional  subpolytopes  stable. 

Next  consider  the  case  k>m.  The  union  of  the  (m-l)-dlmensional 
subpolytopes  of  the  m-dimensional  polytope  9  (defined  in  Theorem  3.1)  are 
compact  and  nowhere  dense  in  HR*);  hence,  there  exist  xo.yo  such  that 
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J 


c={f(x-x0,y-yo) 


x 

.y 


N 

4 


SH 


has  all  its  (m-1 )-dimensional  subpolytopes  stable,  but  Q  is  unstable.  If  n»2, 
let 


g(x,y,z)=f (x,y) 


Otherwise,  define 


g(x,y,z)=f (x,y)+ 


Zn+1  Zn+2  Zn+3 


0 

0 


'2n+l 


"2n 


z 


3n-2 


Zn2-3n+7  zn2-2n+4 


In 


each  case,  zelRn  m.  Also  let  (x,y,w)=g(x,y,  ) ,  where  welR^  2n+^ 


and  e>0.  Note  that  each  g^£  is  affine  and  one-to-one. 
Now  consider  the  k-dimensionai  polytope 


SI} 

If  we  choose  the  matrix  norm  |MS=max|mij | ,  it  follows  that  for  every  vertex  A 

of  there  exists  a  vertex  A  of  Q  such  that  i A-AjJ <e .  Furthermore,  every 

(m-l)-dimensional  subpolytope  of  9^  can  be  expressed  as  a  disjoint  finite 

union  uA,  ,  where  each  A.  is  the  convex  hull  of  m-1  vertices  A . A_  ,  of 

v  v  l  m-i 

Suppose  A, . A  «fl  and  A  . . Am  let 

K*  1  q  q+l  m-1 


?ke’(ekt<Xty’") 
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0=conv{A . A  A  . A  .},  where  Pitch  A.  is  n  vertex  of  c  satisfying 

l  q  q-1  m-l  i 

A. -A.  ■!<£,  and  let  I’=  { A  . A  From  Lemma  4.1.  every 

11  q*l  m-l 

Beconv{A, . A.,  .  }cconv(Our)  satisfies  d(B,0)<c.  Hence,  for  sufficiently 

l  m-i 

small  e,  each  (m-l )-dimensional  subpolytope  of  9  is  stable.  2 

Kc 


5.  Conclusions 

Our  results  demonstrate  to  what  extent  the  techniques  for  checking 
polytope  stability  proposed  in  [2]  can  be  extended  to  the  case  of  nxn 
matrices.  We  have  shown  that,  without  further  information  describing  the 
particular  structure  of  a  polytope.  (2n-4) -dimensional  faces  must  be  checked 
for  stability.  Since  testing  even  one  such  face  can  be  a  formidable  task  when 
n  is  large,  and  since  the  number  of  (2n-4)-dimensional  faces  grows 
exponentially  with  n.  more  work  needs  to  be  done  before  a  computationally 
tractable  algorithm  can  be  devised  for  checking  stability.  It  is  our  hope, 
however,  that  our  work  will  be  useful  as  an  integral  part  of  some  future 
coherent  theory  of  robust  stability. 
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Abstract 

We  consider  the  problem  of  determining  whether  a 
polytope  9  of  nxn  matrices  is  D-stable  --  i.e.  whether 
each  point  in  9  has  all  its  eigenvalues  in  a  given 
nonempty,  open,  convex,  conjugate-symmetric  subset  D 
of  the  complex  plane.  Our  approach  is  to  check 
D-stabillty  of  certain  faces  of  9.  In  particular,  for 
each  D  and  n  we  determine  the  smallest  integer  m  such 
that  D-stability  of  every  m-dimens ional  face 
guarantees  D-stability  of  9. 

l.  Introduction 

Let  DcC  be  nonempty.  open.  convex.  and 
conjugate-symmetric  (symmetric  about  the  real  axis), 
and  define  an  nxn  real  matrix  M  to  be  D-stable  if  each 
eigenvalue  1  of  H  satisfies  4cD;  otherwise.  M  is 
D-unstable.  We  consider  the  problem  of  determining 

whether  certain  subsets  of  IRnxn  consist  entirely  of 
D-stable  matrices.  To  facilitate  discussion  we  begin 
with  some  definitions. 

A  (convex)  polvtope  9  in  a  vector  space  V  is  the 
convex  hull  conv(O)  of  any  nonempty  finite  subset  OcV. 
The  dimension  of  9  is  the  dimension  of  the  affine  hull 
affl?)  of  9.  A  face  of  9  is  any  set  of  the  form  (T'VIP. 
where  (I  is  a  supporting  hyperplane  of  9.  Finally,  a 
k -dimensional  half-plane  in  V  is  any  nonempty  set  of 
the  form  X=RoS.  where  R  is  a  closed  half -space.  S  is  a 
k-dimensional  affine  subspace,  and  S«R.  (Note  that 
this  implies  that  aff(W)  is  simply  S.) 

In  the  robust  control  literature,  considerable 
interest  has  been  generated  by  the  problem  of 
determining  whether  a  family  of  linear  systems  can  be 
shown  to  consist  entirely  of  D-stable  systems  by 
checking  D-stability  of  certain  representative  members 
of  that  family.  In  many  cases,  such  problems  can  be 
reduced  to  that  of  determining  whether  a  polytope  or 

other  subset  of  Rn  or  Rnxn  consists  entirely  of 
D-stable  points  (1).[2).  (D-stability  of  a  vector 

xeR11  means  simply  that  the  polynomial  sn-xnsn*. . ,-xi 
has  all  its  roots  in  D.)  We  are  primarily  interested 
in  the  technique  of  checking  D-stability  of  lower 
dimensional  faces  of  a  polytope  in  order  to  guarantee 
0-stability  of  the  entire  set. 

Most  "facial"  results  pertain  to  continuous-time 
(CT)  stability  --  i.e.  where  D  is  the  open  left  half 
complex  plane.  The  seminal  result  (3)  for  polynomial 
polytopes  motivates  the  approach.  In  (3)  it  is  shown 
that  a  polynomial  polytope  of  a  particular  simple 
structure  (an  "Interval  polynomial")  is  CT  stable 
whenever  four  specially  constructed  vertices  are  CT 
stable.  A  more  recent  result  (1)  demonstrates  that, 
for  an  arbitrary  polynomial  polytope,  checking  all 
edges  is  sufficient  to  guarantee  CT  stability.  With 

nxn 

respect  to  polytopes  in  R  .is  has  been  shown  [4] 
that  1)  an  arbitrary  polytope  is  CT  stable  if  all 
(2n-4)-dimenslonal  faces  are  CT  stable  and  2)  there 
exist  CT  unstable  polytopes  such  that  all 


(2n-D)  dimensional  faces  are  CT  stable:  hence,  the 
value  2n-4  is  minimal.  In  this  paper  we  extend  the 
results  of  [4]  to  D-stability  where  D  may  be  any 
nonempty,  open,  convex  conjugate-symmetric  subset  of 

C 

We  note  that  for  the  cases  n=0  and  n=l  out- 
problem  has  a  trivial  solution-  D-stability  of 
vertices  guarantees  D-stability  of  the  polvtope,  To 
handle  n22  we  need  to  partition  the  family  of 
stability  sets  D  according  to  the  following  two 
assumptions. 

Assumption  A.  0  is  of  the  form  D*{s*c|a<Re  s<b). 
where  -»<a<b£“ 

Assumption  B:  D  is  a  nonempty,  open,  convex, 

conjugate-symmetric  set  not  satisfying  Assumption  A. 

In  addition.  we  define  m.(n)»/^‘  .  n"?  and 

A  |_2n-4.  n>2 

mB(n)=2n-2.  We  Intend  to  show  that  mA  and  mB  are  the 
values  of  m  that  we  seek  for  cases  A  and  B. 


2.  Sufficiency  of  m.  and  mn 

Throughout  our  analysis,  we  will  make  extensive 
use  of  the  fact  that  any  affine,  one  to -one  map 

f:*^  — •  Rn  determines  an  affine  isomorphism  between 
k  k 

R  and  ftR  ).  Among  other  things,  this  implies  that, 
v 

for  any  polytope  9r.IR  .  f!'F)  is  also  a  polytope  of  the 
same  dimension  as  ?:•  furthermore,  f  sets  up  a 
one-to-one  correspondence  between  the  q-dtmensional 
faces  of  9  and  the  q-dimensional  faces  of  f (9*1 .  In 

g 

addition,  f  maps  each  k  dimensionul  half -plane  in  IR 
into  another  k-dimens lor.u ’  half -plane  (e  g  see  (5)1 
Finally,  we  note  that  every  polytope  is  compact  and 

that  any  set.  of  the  form  (xeR^Jiixn^r)  •  where  r>0.  is 

a  poiytope  whose  q-dimensional  faces  are  generated  by- 
fixing  k-q  entries  of  x  at  either  £r  and  letting  the 
remaining  q  entries  vary  independently  over  t  -r . r ] - 

With  these  observations  in  mind,  we  prove  a 
result  characterizing  the  affine  structure  of  the  set 

of  D-unstable  points  in  Rnxn. 

Lemma  2  1  I)  If  D  satisfies  Assumption  A  then  for 

each  D-unstable  M«Rnxn  there  exists  an 

(n1-mA) -dimensional  half-plane  #cRnxn  such  that  a)  M *H 

and  b)  N left  Implies  N  is  D -unstable. 

2)  If  D  satisfies  Assumption  B.  then  for  each 

D-unstable  M«Rnxn  there  exists  an  (n2-mB)-dimensional 

half-plane  #cRnxn  such  that  a)  Me#  and  b)  X«#  implies 
X  is  D-unstable. 


1  This  work  was  supported  In  part  by  NSF  Grant  So. 
ECS-8612948  and  by  AFOSR  Grant  No.  AF0SR-88-0087 . 
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Proof,  if  p  oontu.  ts  -in  unstable  A.  time  exists 
an  (n*'2n-4) -dimens  tonal  half-plane  H  consisting 
entirely  of  unstable  points  and  containing  A. 
From  dimensionality  arguments,  such  a  plane  must 
Intersect  a  (2n-4)  -dimensional  face  of  9.  iSee 
(4)  for  details. )  ? 


3.  Minimality  of  m*2n-4 

In  this  section  we  show  that,  for  every  integer 

n.  there  exists  a  polytope  T’eR0*11  containing  an 
unstable  point  and  such  that  all 
(2n~5)-dimensional  subpolytopes  of  9  are  stable. 
Hence,  we  conclude  that  checking  stability  of 
k-dimensional  subpolytopes  of  9,  for  any  kc2n-4 
is,  in  general,  not  sufficient  to  guarantee 
stability  of  9. 

Consider  the  polytope 


*■<  -1 


0  -y 


l  x  y  -I  J’ 

A  routine  calculation  shows  that  9  has 


n  **4 

characteristic  polynomial  p(s)«(s*l)  Ms), 
where 

M s )  *s4 *2s3 » ( 2  t-xTx*yTy ) s2 

T  T  T  T  T  2 

-(2-x  x-y  y)s-l*x  xy  y-(x  yi 

From  the  Schwartz  inequality,  it  is  clear  that 
all  coefficients  are  strictly  positive.  The 
corresponding  4x4  Hurwitz  matrix  has  its  leading 
principal  3x3  minor  equal  to 

M3  ( x ,  y )  *4xTx  *-4yTy  *4  ( xTy ) 2  ♦  ( xTx-yTy ) 2 

Clearly.  M3S0  with  equality  iff  x«y*0.  Thus.  9 

consists  entirely  of  stable  points,  except  for 
the  relative  interior  point  corresponding  to 
x»y*0.  We  conclude  that  checking 

(2n-5) -dimensional  faces  (in  this  case  the  entire 
boundary  of  9)  is  insufficient  to  guarantee 
stability. 


Comments  1)  The  preceding  example  can  be 
strengthened  by  adding  si  to  9,  where  e  is 
sufficiently  small.  This  yields  a  polytope  with 
stable  (2n-5)-dlmensional  boundary,  but 
containing  a  ball  of  strictly  unstable  points. 

2)  Since  the  union  of  all  (2n-5)-dimensional 
subpolytopes  is  nowhere  dense,  shifting  the 


parameter  set 


SI  by  an  arbitrarily  small 


vector  yields  an  unstable  polytope  with  all 
(2n-5)'dimenslonal  subpolytopes  stable. 

3)  The  poly  tope  9  described  above  can  be 
transformed  into  a  similar  example  with  any  given 
dimension  either  by  removing  parameters  or  by 
using  9  as  a  face  of  a  higher  dimensional 
polytope. 

Note  that  the  constructions  described  in  1). 
2).  and  3)  can  be  carried  out  simultaneously  to 
giveh  a  stronger  but  algebraically  messy  version 
of  the  minimality  proof  offered  above. 


'I  -■  -:o« 

no  tv  -noun  tha’  a  j 

l"r“5|‘;  sue!*,  thar  st.ibil.tv  uf 


tr.<*  saai.ivs.t 
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subpoly  topes  ot  a  given  ...oiytotin 

guar.mto.s  .-,tib‘iiry  ot  9.  turtnereere.  wv  na-e 
demonstrated  that  chocking  m  dimensional  >> 

always  sufficient.  Tnis  reduces  the  task 
determining  whether  a  poiytop*  is  stable  to  cna* 
ct  deciding  whether  several  low-dimensionai 
polytopes  are  stable.  Our  result  has  certain 
theoretical  significance:  however,  more  work 
needs  to  be  done  before  it  can  be  decided  whether 
the  result  will  help  to  reduce  the  computational 
burden  inherent  in  robust  system  design. 
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IThis  work  was  supported  in  part  by  NSF  Grant  No.  ECS-8612948  and  by  AFOSR 
Grant  No.  AFOSR-88-0087. 


1.  I  lit  rmliict  mu 

Let  0c£  be  nonempty,  open,  convex,  and  conjugate-symmetric:  (symmetric 
about  the  real  axis),  and  define  an  nxn  real  matrix  M  to  be  p-^table  if  each 
eigenvalue  1  of  M  satisfies  otherwise.  M  is  D-unstable .  We  consider  the 
problem  of  determining  whether  certain  subsets  of  (Rnxn  consist  entirely  of 
D-stable  matrices.  To  facilitate  discussion  we  begin  with  some  definitions. 

A  (convex)  polytope  ?  in  a  vector  space  V  is  the  convex  hull  conv(fl)  of 
any  nonempty  finite  subset  OcV.  The  dimension  of  9  is  the  dimension  of  the 
affine  hull  aff(9>)  of  9.  The  relative  boundary  of  9  is  the  boundary  of  9  as  a 
subset  of  the  topological  space  aft(9).  A  face  of  9  is  any  set  of  the  form 
nn7>,  where  FI  is  a  supporting  hyperplane  of  9.  Finally,  a  k-dimensional 
half-plane  in  V  is  any  nonempty  set  of  the  form  W=Rns,  where  R  is  a  closed 
half-space,  S  is  a  k-dimensional  affine  subspace,  and  S£R.  (Note  that  this 
implies  that  aff(#)  is  simply  S.) 

In  the  robust  control  literature,  considerable  interest  has  been 
generated  by  the  problem  of  determining  whether  a  family  of  linear  systems  can 
be  shown  to  consist  entirely  of  D-stable  systems  by  checking  D-stability  of 
certain  representative  members  of  that  family.  In  many  cases,  such  problems 
can  be  reduced  to  that  of  determining  whether  a  polytope  or  other  subset  of  !Rn 
or  IRn*n  consists  entirely  of  D-stable  points  [1],[2].  (D-stability  of  a 
vector  xelRn  means  simply  that  the  polynomial  sn+xnsn*l«-. .  . +xx  has  all  its 
roots  in  D,  where  Xi  is  the  ith  entry  of  x.)  We  are  primarily  interested  in 
the  technique  of  checking  D-stability  of  lower  dimensional  faces  of  a  polytope 
in  order  to  guarantee  D-stabillty  of  the  entire  set. 

Most  "facial"  results  pertain  to  continuous-time  (CT)  stability  --  i.e. 
where  D  is  the  open  left  half  complex  plane.  The  seminal  result  [3]  for 
polynomial  polytopes  motivates  the  approach.  In  [3 ]  it  is  shown  that  a 
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polynomial  poly  tope,  of  a  particular  simple  structure’  (an  "  interval 
polynomial")  is  CT  stable  whenever  four  specially  constructed  vertices  are  (if 
stable.  A  more  recent  result  £  1 j  demonstrates  that,  for  an  arbitrary 

polynomial  polytope,  checking  all  edges  is  sufficient  to  guarantee  CT 
stability.  With  respect  to  polytopes  in  Rn'n.  it  has  been  shown  [4]  that 
1)  an  arbitrary  polytope  is  CT  stable  if  all.  (2n-4)-dimensional  faces  are  CT 
stable  and  2)  there  exist  CT  unstable  polytopes  such  that  all 

(2n-5)-dimensional  faces  are  CT  stable;  hence,  the  value  2n-4  is  minimal.  In 
this  paper  we  extend  the  results  of  [4]  to  O-stability  where  D  may  be  any 
nonempty,  open,  convex,  conjugate-symmetric  subset  of  C. 

We  note  that  for  the  cases  n=0  and  n*l  our  problem  has  a  trivial 
solution:  O-stability  of  vertices  guarantees  D-stability  of  the  polytope.  To 
handle  n22  we  need  to  partition  the  family  of  stability  sets  D  according  to 
the  following  two  assumptions. 

Assumption  A:  0  is  of  the  form  D={s«c|a<Re  s <b} ,  where  -»£a<b£*». 

Assumption  B:  D  is  a  nonempty,  open,  convex,  conjugate-symmetric  set  not 
satisfying  Assumption  A.  — 

In  addition,  we  define  %(n)3^2n-4  n>2  an<*  mi)(n)=2n~2,  We  intend  to  show 

that  m.  and  m_  are  the  values  of  m  that  we  'seek  for  cases  A  and  8. 


2.  Sufficiency  of  m.  and  m. 


Throughout  our  analysis,  we  will  make  extensive  use  of  the  fact  that  any 

affine,  one-to-one  map  f — »  IRn  determines  an  affine  isomorphism  between 

k  k 

and  f(R  ).  Among  other  things,  this  implies  that,  for  any  polytope  ?cR  ,  f(?) 


is  also  a  polytope  of  the  same  dimension  as  ?;  furthermore,  f  sets  up  a 
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otio  to  ono  correspondence  between  the  q-dlmensional  faces  of  P  .nm  i.i'.o 

q-dimonsional  faces  of  f(9*).  In  addition,  f  maps  each  k- dimensional 
u 

half -plane  in  R  into  another  k-dimensional  half-plane  (e.g.  see  15 J ) , 
Finally,  we  note  that  every  polytope  is  compact  and  that  any  set  of  the  form 
{xalR^I  nxil^Sr} .  where  r>0,  is  a  polytope  whose  q-dimensional  faces  are 
generated  by  fixing  k-q  entries  of  x  at  either  +7  and  letting  the  remaining  q 
entries  vary  independently  over  [ -r . r 1  * 

With  these  observations  in  mind,  we  prove  a  result  characterizing  the 
affine  structure  of  the  set  of  D-unstabie  points  in  Rnxn. 

Lemma  2 . 1  If  D  satisfies  Assumption  A  (respectively,  Assumption  B),  then  for 

each  D-unstable  M*IRnxn  there  exists  an  (n2-m. ) -dimensional  (respectively, 

A 

(n2-mD) -dimensional)  half-plane  #cRnxn  such  that  a)  MeJf  and  b)  N«#  implies  N 

D 

is  D-unstable. 

Proof  Suppose  Assumption  A  holds.  If  a*-*,b=«,  the  statement  is  vacuously 
true;  otherwise,  we  need  to  consider  two  cases. 

Case  I  --  M  has  a  real  eigenvalue  ,*o<D:  Let  T-[v  WJ ,  where  v  is  an 
eigenvector  corresponding  to  /*o  and  W  is  chosen  to  make  T  nonsingulaf. 
Clearly,  the  map  Rnxn  determined  by  f(A,y,Z)»T^  is  affine 

and  one-to-one.  If  -*o<a,  let  i  =  ( -«*, ]  and  let  it  be  the  (n2-n+l)-dimensional 
half-plane  »*{f(^,y,Z)|yl«I,y«R1Xn"1,Z«Rn“lxn“1>.  If  a=  — ,  then  >*0>b  so  we  set 
I*[^o.w)  and  construct  it  in  the  same  way.  In  either  case,  Me#  and  every 
matrix  in  it  is  D-unstable.  Since  n2-n+l?.n2-mA,  it  remains  to  select  any 
(na-mA)-dimensional  half-plane  #  satisfying  Me#c#. 

Case  II  --  M  has  a  complex  eigenvalue  pair  aotiflo  with  «o>0:  Let 
T*(u  v  W],  where  u+iv  is  an  eigenvector  corresponding  to  ao+ij3o  and  W  is 
chosen  to  make  T  nonsingular.  If  ao<a,  let  it  be  the  (n2-2n+4)-dimenslonal 
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half  plane 


fX  Y1 

«-(T[o  ‘]t 


1  trS<2a0.Y«R2'n  ^.ZelR11  '3'n  ‘} 


l  Lr  X<2a0  describes  a 


4-dimensional  half-plane,  since  tr  X-<X,I>.)  H  contains  only  D-unstable 
points,  since  trXS2ao  implies  X  has  at  least  one  eigenvalue  A  with  Re  A&xq. 
If  ao>b,  #  is  defined  by  tr  X£2a0,  and  the  same  reasoning  holds.  In  either 
case,  M «#. 

Now  suppose  Assumption  B  holds.  We  again  consider  two  cases.  Case  I  •- 
M  has  a  real  eigenvalue  Again  let  T*[v  W 1 ,  where  v  is  an  eigenvector 
corresponding  to  >*o-  Since  D  is  convex,  either  (-»,^o]nD*0  or  [;»o.*)nDa0‘  In 
the  former  case,  let  it  be  the  (n2-n+l  )-dimensional  half-plane 
3»-{t[q  ^Jt  1  j^A0,yelRlxn  ^ZelR0  lxn  1).  For  the  latter  case,  alter  the 
definition  of  #  by  substituting  for  "ASA o" .  In  either  case,  M and 
every  matrix  in  it  is  D-unstable.  Since  n2-n+l>n2-mB,  it  remains  to  select  any 
(n2-i»B) -dimensional  half -plane  H  satisfying  M«#c#. 

Case  II  —  M  has  a  complex  eigenvalue  pair  ao±i0o*D:  Let  T=[u  v  Wj, 
where  u*-iv  is  an  eigenvector  corresponding  to  ao«-i0o-  Since  D  is  convex, 
there  exists  a  half -space  IIcC  such  that  ao+i0o«n  and  IlnD*0.  Let  #  be  the 
(n2-2n+2)-dimensional  half-plane 

»=•{ T  -0  a  y  T'1|a+i0«n;x,yelRlxn"2;Z«Rn"2xn"2} 

0  0  ZJ 

Clearly,  H  contains  only  D-unstable  points,  and  Me#.  Q 


Next  we  prove  an  easy  result  concerning  the  intersection  of  affine  sets 


Lemma  2.2  Let  V  be  a  p-dimensional  Euclidean  space,  #cv  a  k-dimensional 
half-plane,  and  r  a  q-dimensionai  affine  subspace  with  k+q>p,  Consider  any 
vector  XoS#nr.  There  exists  a  (k*q-p)-dimenslonal  half-plane  #  such  that 
xoe#c#nr. 
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Proof  By  defiuiMon.  tf-=R'\S,  where  R  is  a  closed  half-space  amt  3  is  ,t 
k-dimensional  affine  subspace  satisfying  S<R.  There  exists  an  affine  subspace 
Scsnr  with  dim  S*=k+q-p  and  xo«S.  If  ScR,  let  ScRns  oe  any  (kfq-p)  dimens iona! 
half-space  containing  xo.  Then  #cRnS''r«#or .  If  S<R,  let  H-Rns,  Then  Xo6#. 
since  xoStfnrcR,  Also,  dim jT=k+q-p,  since  H  is  nonempty,  " 

We  are  now  in  a  position  to  prove  our  first  main  result. 

Theorem  2.3  Under  Assumption  A  (respectively,  Assumption  B),  D-stability  of 
every  matrix  in  every  mA-dimensional  (respectively,  nig-dimensional)  face  of  9 
guarantees  D-stability  of  every  matrix  in  9. 

Proof  Suppose  Assumption  A  holds.  Our  arguments  here  are  similar  to  those 

used  in  [2, Lemma  l].  If  9*  is  a  D-unstable  polytope  of  dimension  k>mA>  there 

exists  a  D-unstable  matrix  Mi*?*.  From  Lemma  2.1,  there  is  an 

(nJ-m. ) -dimensional  half  plane  Hi,  consisting  entirely  of  D-unstable  points 
A 

and  containing  Mi.  Since  Hi  is  unbounded,  there  exists  an  M2«»i  lying  on  the 

boundary  of  9 *  and,  hence,  in  one  of  its  (k-1) -dimensional  faces  ?k-i-  From 

Lemma  2.2,  the  intersection  #tnaff(?k-i)  contains  a  (k-mA-l)-dimensional 

half-plane  Hz  such  that  M2*Hz.  Proceeding  inductively,  we  find  that  there 

exists  an  m  -dimensional  face  9m  and  a  point  Mk-ra«?ra  such  that  Mk-m  is 
A 

D-unstable. 

Under  Assumption  B,  the  same  proof  holds  if  we  replace  mA  by  nig.  jj 

3.  Minimality  of  m.  and  m„ 

-  A  D 

Our  next  task  is  to  show  that  mA  and  m8  are  the  smallest  integers  such 
that  D-stability  of  all  mA-dimensional  or  nig-dimensional  faces  of  9  guarantees 
D-stability  of  9  under  Assumptions  A  or  B,  respectively.  In  order  to  prove 
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this  wo  not'd  a  lomma  which  may  bo  interpret oci  as  a  multivariable  extension  ol 
L'Hospital's  rule.  For  any  k*k  matrices  Q  and  R,  we  use  the  notation  Q>0  and 
R<0  to  signify  that  Q  is  positive  definite  and  R  is  negative  definite, 
respectively. 


Lemma  3.1 


If 

Let  OeUclR4  with  U  open,  and  let 


(R2  be  C2 


functions . 


in 


addition,  suppose  ei(0)*e2(0)*0, 


de2 

-  0, 

a*ei 

=  o, 

a*e2 

ax 

ax* 

ax* 

x=0 

X=0 

X=0 

<  0 


x=o 


For  every  <5> 0  there  exists  an  e>0  such  that  0*iixii<e  implies  e a ( x ) < | e t ( x } ; . 
Proof  From  [6 , p . 340 J ,  for  every  Q>0  there  exists  an  e>0  si  a  that  iixiKi 
implies 


e i ( x )  -  ei(0)  - 


3£i  „  _  ixTa£ej. 

ax  x=Q  2  ax*  |  * 


Setting  Q  *  - 


T 

x  Qx 


<  £( - ) 

2 1 1+«5 ; 


yields 


x=0 


et(x) j  <  8  e2(x)+|xTQx 


T 

x  Qx 


T 

x  Qx 


<  6 


1  T 

and,  from  (1),  e2(x)  <  (5-^)x  Qx  <  0  for  x*0.  Hence,  for  x*0, 


ei(x) 


e2(x) 


iel(x? I 


(e2(x)t|xTQx)  -  |xTQx 


'ei(x)  I 
T 

x  Qx 


(1) 


ej(x)+jXTQx  1 


T 

x  Qx 


2V1+<T 

1  .  iiJL) 

2  2 


S 


Thus,  e2(x)<-j|ei(x) | 


Now  we  can  prove  our  second  main  result. 
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t 


riHiomn  .! .  1  Suppose  i)  satisfies  Assumption  A  (respectively,  Assumption  b). 
For  each  n  there  exists  an  m  dimensional  (respectively,  mQ-dimens iona  1 ! 
polytope  ?c|Rnxn  containing  a  D-unstable  point  and  such  that  all 
(m  -1 ) -dimensional  (respectively,  (m  -1 ) -dimensional )  faces  of  9  are  D-stable. 

A  B 

Proof  Under  Assumption  A  we  need  to  consider  nine  cases. 


Case  I  —  a>-«,b<»,n=2:  Consider  the  affine,  one-to-one  map 

b  x  ' 
a+b 

x  TJ 

and  the  corresponding  one-dimensional  polytope  9-{f (x) j |x|Sl> .  The  point  in  9 

corresponding  to  x=0  is  clearly  D-unstable.  It  suffices  to  prove  that  the 

characteristic  polynomials  A+  and  A-  of  f ( x ) -bl  and  al-f(x),  respectively,  are 

Hurwitz  for  all  x*0.  This  is  in  fact  true,  since  A  ( s )  =s +^jj“S+x  and 

2 

23  (b-a)  2 

A  (s)=s  +-(b-a)s+ — - —  +x  have  positive  coefficients  for  x*0. 


Case  II  —  a>-«,b<»,n=3:  Let 


f(x.y)  * 


b  1  -x 
-1  b  -y 
a+b 


It  is  sufficient  to  show  that  the  characteristic  polynomials  A  and  A  of 


f(x,y)-bt  and  al-f(x.y),  respectively,  are  Hurwitz  for  *0.  A 
straightforward  calculation  yields  A+(s)=s3+^~s2+(l+x2+y2)s+^~  and 
A  (s)=s3+|(b-a)s2+(l+x2+y2+2(b-a)2)s+^~(l+2x2+2y2+(b-a)2) .  Each  polynomial 
has  positive  coefficients  for  *j*0.  The  fact  that  they  are  Hurwitz  follows 

^  u  *  2  O 

from  positivity  of  the  second-order  leading  principal  minors  Mg  *  -~(x  ) 

—  K.a  22  2 

and  Mg  =  -j-(4+3x  +3y  +9 (b-a)  )  of  the  corresponding  3x3  Hurwitz  matrices. 

T  ' 

b  1  -x 

T 

-1  b  -y 

x  y  (^)lj 

where  x,y,«IRn  2.  A  tedious  calculation  shows  that  A+(s)  =  (s+^g^)n  4A+(s)  and 


Case  III  --  b>-<»,a<»,n>3:  Let  f(x,y)  = 
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f 


A  (s) 


,  b-a,  n 


A  (s).  where 


A*(s)  =  s4+(b-a)s3+(  l+xTx*yTy+(^p)2)s2+~^(2*xTx+yTy)s+xTxyIy-(xTyl2>-(‘;^-)2 

b  W  <X 


a"  ( s )  =  s4 +3 ( b-a ) s3+ ( l+xTx+yTy+~(  b-a ) 2 ) s2 + ( b-a )  ( 1  +|xTx+|yTy+|( b-a ) 2 )  s 

^xTxyTy-(xTy)2+{~)2(1^2xTx^2yTy>(b-a)2) 


From  the  Schwartz  inequality,  A+  and  A  have  positive  coefficients  when  [yj*°- 
Furthermore,  the  third-order  leading  principal  minors  of  the  corresponding  4x4 
Hurwitz  matrices  of  A+  and  A~  are 


M3(s)  *  (^)2(2(l^(~)2)(xTx^yTy)MxTx-yTy)2f4(xTy)2) 


M3(s) 


( 2+§xTx+|y  Vf  ( xTx-yTy ) 2+ ( xTy ) 2 ) ( b-a ) 2 


.27  63  T  63  T  .  ..  .4  81,.  ,6 

4 +  8X  x*Ty  y)(b”a)  +T(b“a) 


Since  and  Mg  are  positive,  A+  and  A  and,  hence,  A*  and  A  are  Hurwitz. 

The  remaining  six  cases  are  handled  similarly  by  choosing  all  eigenvalues 
in  the  interior  of  D,  except  for  one  or  two  on  the  boundary  of  D.  For 
example,  for  a>-»,  b»»,  n>3,  set 


f(x,y) 


a  1  -x 
-1  a  -y 


l  x  y  (a+i) Ij 

Adopting  Assumption  B,  suppose  0  is  not  of  the  form  {s|a<Res<b>.  Since  D  is 
convex,  there  exists  a  real  ao<D  such  that  the  line  L={ao+i/s|/3«R}  satisfies 
Li*D.  Since  D  is  conjugate  symmetric  and  open,  there  exists  a  (3o>0  such  that 
<xo±l(3o  are  boundary  points  of  D,  but  «o±13*D  when  ||3|<0o>  Furthermore,  there 
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i 


exists  a  6> 0  such  that  uoidel).  Again  invoking  convexity,  the  open  diamond 

d."int  conv{ao±<5^o .«o±i(3o }  is  contained  in  D.  To  simplify  the  problem, 
0 

consider  the  open  diamond  d.-— -(d.-ao )  =  int  conv{+i,+5}.  We  need  only 

o  Po  o 

construct  a  single  polytope  9  containing  a  matrix  with  a  pair  of  eigenvalues 
at  ±i  and  with  all  mB-dimensional  faces  consisting  of  matrices  with  all 
eigenvalues  in  d.:  then  |3o?+aoI  would  then  satisfy  the  desired  properties  with 


respect  to  d_. 

0 


Consider  the  (n2-2n+2)-dimensional  polytope 

T, 


-  < 


w  l+x  y 
l+x  -w  z 

L  y  z  o 


5  «} 


where  y,z«IR 


n-2 


Clearly,  9,  has  a  D-unstable  point  M  at  w=x*0 . y*z=0.  We  will 


show  that  for  sufficiently  small  e,  every  pint  in  9  except  M  is  D-stable. 

Hence,  9=9  satisfies  the  desired  properties. 

& 


Case  I  —  n=2: 


Each  point  in  9  has  characteristic  polynomial 

6 


2  2  2  2  2  1/2 
A(w,x,s)=s  +l-w  -x  and  hence  has  eigenvalues  ±i(l-w  -x  )  .  Let 


Case  II  —  n=3:  Each  point  in  9  has  characteristic  polynomial 

C 


A(w,x,y,z,s) 


=  s3+(l-w2-x2-y2-z2)s-(w(y2-z2) +2xyz) 


Let  g(w,x,y,z,«,0)  =  ^ 
polynomial  function  and,  hence,  analytic. 


Re  A(w,x,y,z,a+i)3) 
Im  A(w,x,y,z ,a+i0) 


..  ag... 

3(«.0) 


0 

0 

0 

0 

0 


It  is  easy  to  see  that  g  is  a 
A  straightforward  calculation  shows 


11 
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finis .  from  the  implicit  function  theorem,  there  exists  a  unique  analytic 


function  h:U  — *  !R  such  that  h(0) 
T 

[w  x  y  z]  eU. 


f° 

'[i. 


and  g(w,x,y,z.h(w,x,y.z) )  =o  for  every 


Next,  let  [el]=h~[i]’  A  tedious  computation  shows 

.2. 


3e 


3(w,x,y.z) 


3e2 


3(w,x,y,z) 


=  0, 


3  e  i 


3(w,x,y,z) 


=0. 


,2 

3  ea 


3(w,x,y.z)‘ 


-  I 


From  Lemma  3.1,  there  exists  an  e>0  such  that  e2(w,x,y,z)<-riei(w,x,y,z) 


whenever  0*n[w  x  y  z]lil<«.  Since  ei  and  e2  are  continuous,  we  may  also  assume 

; e i  j <1 ;  i*l,2.  Returning  to  it  follows  that 

h2(w,x.y,z)<l-||hl(w,x,y,z) | .  |h2(w,x,y,z) |<1 ,  and  (hi (w,x,y,z)-l |<1  for  all 

T 

[w  x  y  z]  *o.  Hence,  hed^. 

Case  III  —  n>3:  We  have 


4  22TT2  T  T  T  TT  T2 

A(w,x,y,z,s)  =  s  +•  ( 1-w  -x  -x  x-y  y)s  -(wy  y-2xy  z-wz  z)s+y  yz  z-{y  z) 


■  [*a  Again>  * u  a  pois'n°i,iai  tuncti°»: 


in  this  case 


Jg_ 


3(«,U) 


[0  -2 

[2  0 


»2n-2 


Thus  there  exists  an  open  UdR  with  OeU  and  h:U  — *  !R  such  that  h(0)  = 
and  g(w,x,y,z,h(w,x,y,z))*0  for  every  [w  x  y  z]T«U.  Let  l] «h- [®] .  Then 


3b  1 


d(w,x,y,z) 


3ea 


3{w,x,y,z) 


=  0, 


32et 


3{w,x,y,z) 


=  0, 


3(w,x,y,z) 


-10- 


Applying  Lemma  3.1  as  in  Case  II.  it  follows  that  hlw.x.y.zjen.  for  every 

d 

r 

[w  x  y  z]  !“0. 

Note  that  Theorem  3.2  also  implies  that  the  half-planes  constructed  In 
Lemma  2.1  are  maximal  in  the  sense  that  there  exists  a  D-unstable  matrix  M  in 
Rnxn  such  that  every  half-plane  containing  M  of  dimension  greater  than  n2-mA 
or  n2-r#B  must  also  contain  a  D-unstable  matrix.  Indeed,  if  this  were  not  the 
case,  the  arguments  in  Theorem  3.2  could  be  used  to  prove  that  mA  and  mg  are 
not  minimal., 

4.  Conclusions 

Our  results  demonstrate  to  what  extent  the  techniques  for  checking 
polytope  stability  proposed  in  [1]  and  [3]  can  be  extended  to  the  case  of  nxn 
matrices.  We  have  shown  that,  without  further  information  describing  the 
particular  structure  of  a  polytope,  either  (2n-4) -dimensional  or 
(2n-2) -dimensional  faces  need  to  be  checked  for  D-stability,  depending  on  the 
structure  of  D.  Since  testing  even  one  such  face  can  be  a  formidable  task 
when  n  is  large,  and  since  the  number  of  (2n-4)  -dimensional  and 
(2n-2)-dlmensional  faces  grow  exponentially  with  n,  more  work  needs  to  be  done 
before  a  computationally  tractable  algorithm  can  be  devised  for  checking 
D-stability.  It  is  our  hope,  however,  that  our  work  will  be  useful  as  an 
integral  part  of  some  future  coherent  theory  of  robust  stability. 
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Abstract 

In  this  paper  we  present  a  theory  which  characterizes  LTI  state-space 
realizations  of  perturbed  rational  transfer  function  matrices.  Our  approach 
is  to  model  system  perturbations  as  sequences  in  the  space  of  rational 
matrices.  First,  we  give  a  definition  of  convergence  in  the  space  of  rational 
matrices  which  is  motivated  by  the  kinds  of  parameter  uncertainties  occurring 
in  many  robust  control  problems.  A  realization  theory  is  then  established 
under  the  constraint  that  the  realization  of  any  convergent  sequence  of 
rational  matrices  should  also  be  convergent.  Next,  we  consider  the  issue  of 
minimality  of  realizations  and  propose  a  method  for  calculating  the  dimension 
of  a  minimal  realization  of  a  given  transfer  matrix  sequence.  Finally, 
necessary  and  sufficient  conditions  are  discussed  under  which  a  sequence  of 
state-space  systems  is  a  minimal  realization  and  under  which  minimal 
realizations  of  the  same  transfer  function  sequence  are  state-space 
equivalent.  Relationships  with  standard  algebraic  system  theoretic  results 
are  discussed. 


xThis  work  was  supported  by  AFOSR  Grant  No.  AF0SR-88-0087 . 


V  '*> 

1_.  I  n  t  ro  due  t  ion 

The  theory  of  state-space  realizations  fur  strictly  proper  rational 
matrices  has  been  thoroughly  studied  (e.g.  see  [H5j).  More  ••ecently. 
techniques  for  handling  improper  transfer  matrices  have  been  devised  (see 
[14|>.  In  this  paper  we  extend  those  ideas  to  the  perturbational  case  •  i.e. 
where  a  system  is  described  by  a  convergent  sequence  of  rational  matrices 
(possibly  improper).  A  realization  is  then  a  sequence  of  (generalized) 
state-space  systems.  The  problem  is  made  nontrivial  by  imposing  the 
constraint  that  the  matrix  entries  of  the  realization  sequence  should  also 
converge . 

Part  of  our  motivation  for  this  problem  comes  from  the  study  of  robust 
control  problems  --  specifically  from  those  dealing  with  order  uncertainty  and 
singular  perturbations.  For  example,  the  robustness  problems  addressed  in 
L  l  ]  - [ 3 ]  are  based  on  singularly  perturbed  system  models.  Physical  systems  are 
invariably  subject  to  some  variations  in  parameters,  often  resulting  in 
changes  in  model  order.  It  is  desirable,  therefore,  to  design  compensators 
which  meet  performance  criteria  independent  of  system  perturbations.  Many 
robust  control  theories  (e.g,  [4J)  emphasize  input-output  performance 
characteristics.  Our  intention  is  to  develop  some  fundamental  tools  for 
examining  robustness  problems  associated  with  a  system's  internal  structural 
properties. 

One  way  to  approach  this  problem  might  be  througli  the  application  of 
algebraic  system  theory  (see, e.g. . [17] ) .  In  this  setting,  the  transfer 
function  sequence  is  viewed  as  a  rational  function  over  the  ring  c  of 
convergent  real  sequences  using  pointwise  operations.  Unfortunately,  we  will 
see  that  existing  results  in  algebraic  realization  theory  apply  to  our  case 
only  marginally.  This  is  due  to  three  key  facts:  t)  The  ring  c  is  not  an 
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integral  domain  2)  Most  results  in  algebraic  realization  theory  deni  only 
with  the  case  of  proper  transfer  functions.  3)  An  abstract,  version  of  the 
Weierstrass  decomposition  for  matrix  pencils  over  a  ring  does  not  yet  exist. 
Nevertheless,  our  feeling  is  that  the  properties  of  sequences  of  transfer 
functions  are  sufficiently  important,  from  the  point  of  view  of  robust  control 


theory  that  they 

deserve 

separate 

treatment,  not  only 

for 

the 

sake  of 

mimicking  standard 

resit  Its 

from  the 

algebraic  theory,  but 

also 

in 

order  to 

obtain  deeper  insight  into  the  specific  structure  of  realizations  over  this 
particular  ring. 

From  an  analytic  perspective.  considerable  work  dealing  with 
perturbations  of  rational  matrices  has  appeared  (e.g.,  [4 ] , [7] -[ 12 ] ) .  In 
these  papers  various  rational  matrix  topologies  have  been  proposed,  motivated 
by  a  variety  of  control  problems.  The  closest  of  these  to  our  work  are 
[10] -[12],  where  a  singular  perturbation  theory  for  transfer  functions  is 
developed  and  a  specific  form  of  realization  is  given.  However,  [10]-(12]  do 
not  explicitly  address  those  problems  dealing  with  the  existence  of 
realizations  in  general  and,  in  particular,  the  minimal  realization  of 
perturbed  systems.  In  [4]  rational  matrix  convergence  is  characterized  in 
terms  of  the  "graph  metric"  which  is  used  to  address  certain  problems  in  local 
simultaneous  stabilization.-  It  is  easy  to  show  that  the  graph  metric  induces  a 
topology  which  is  very  different  from  that  corresponding  to  simple  system 
parameter  convergence.  The  work  of  [7]  and  [9]  also  treats  the  orobiem  of 
topologizing  the  set  of  rational  matrices  and  is  closely  related  to  ours,  but 
again  does  not  examine  the  realization  problem.  Our  work  is  motivated  solely 
by  realization  and  robustness  issues;  our  constructions  are  designed  to  yield 
the  simplest  definition  of  convergence  corresponding  to  convergence  of  system 
parameters . 
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We  are  mainly  concerned  with  the  interplay  between  two  types  of  f.n 


system  representations.  First,  let  IR(s)  be  the  set  of  all  rational  functions 

rm  $  **' 

over  IR.  and  let  R(s)  be  the  set  of  r*m  matrices  over  IR(s).  Next,  consider 


(generalized)  state-space  systems 


Ex  =  Ax  *  Bu 


y  =  Cx. 


where  E  and  A  are  nxn  real  matrices  satisfying  the  standard  regularity 
assumption  det(s£-A)*0.  B  is  nxm,  C  is  rxn.  and  E  may  be  singular.  For  the 
sake  of  brevity,  we  identify  the  system  (1)  with  the  matrix  4-tuple 
(T'(E,A.B.C)elRn(,"!l  m  l>.  The  transfer  matrix  of  (1)  is 

H(SI  •  C(SE-A)  8  .  •  '21 

Throughout  the  paper  we  assume  that  the  values  of  m  and  r  are  fixed;  we 
consider  n  to  be  a  variable. 

Definition  1.1  1)  A  state-space  system  aelRn*n^m+r  1  is  said  to  have  dimension 

n.  In  this  case,  we  write  dima=n. 

2)  If  a  rational  matr^*  H  is  of  the  form  (2).  we  say  that  (E.A.B.C)  is  a 


realization  of  H. 


With  regard  to  (1)  and  (2),  a  (nonperturbational )  realization  theory 
already  appears  in  [14].  We  now  summarize  the  main  results  of  this  theory. 
Theorem  l . 2  [14] 

1)  Every  rational  matrix  has  a  realization. 

2)  The  minimal  dimension  over  all  realizations  of  H,  denoted  * ) .  is 

jui  H(  s ) )  =  t>(Hg(s))  +  v(^Hf(i)). 

where  u(-)  is  MacMillan  degree,  and  Hg  and  H^.  are  the  unique  strictly 

proper  rational  matrix  and  polynomial  matrix,  respectively,  satisfying 

H=H  +H.. 
s  f 
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4 


3)  A  4  -tuple  a  is  a  minimal  realization  of  some  rational  matrix  H  if  and 
only  if  o  is  controllable  and  observable  (as  defined  in  [8J). 

4)  If  <r1  =  (E1  .Aj  ,B1  ,Cj )  and  o2  =  (  E0  .A0  .  B,, .  )  are  minimal  realizations  of  the 

same  rational  matrix,  there  exist  nonsingular  matrices  M  and  N  such  that 
E  =ME.N .  A  =MA.N.  B,,=MBt.  and  C  =C.N. 

The  results  of  our  present  work  may  be  considered  to  be  a  generalization  of 
theorem  1.2  to  the  case  of  rational  matrix  sequences  (H^). 

In  Section  2  we  will  choose  a  natural  definition  for  convergence  of 
rational  matrices.  Working  from  this  definition,  we  will  consider  sequences 
H^— *H  in  IR(s)rm  and  attempt  to  characterize  those  sequences  {o^}  in  Rn^2n  m  l* 
such  that  1)  converges  to  some  rr  in  the  matrix  sense.  2)  a ^  is  a 
realization  of  for  sufficiently  large  k.  and  3)  a  is  a  realization  of  H. 
We  view  this  approach  as  a  way  of  modeling  the  possible  perturbations  in  the 
internal  structure  of  a  system  corresponding  to  a  given  perturbation  in  the 
input-output  description  { ) . 

In  our  realization  theory,  we  will  see  that  Theorem  1.2.  part  1)  remains 
true  (Section  5).  Corresponding  to  the  expression  for  /u  in  part  2).  in 
Section  6  we  will  define  and  give  an  explicit  expression  for  a  degree  function 
which  equals  the  dimension  of  all  "minimal  realizations"  of  a  sequence  of 
transfer  matrices.  It  will  be  shown  that,  properties  3)  and  4)  do  not  hold  as 
stated  for  sequences;  however,  we  will  discuss  important  special  cases  where 
similar  statements  do  hold.  In  Sections  5  and  6  we  will  also  discuss  the 
connections  between  our  work  and  the  standard  algebraic  realization  theory 
(see  [17]). 
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2.  Convergence  in  the  Space  of  Rational  Matrices 

We  first  consider  the  problem  of  defining  a  topology  on  IR(s)  and  later  on 
IR(s)tm.  Convergence  in  IR(s)  will  be  defined  in  the  most  natural  way  such  that 
small  perturbations  in  IR(s)  correspond  to  small  perturbations  in  the 
coefficients  of  numerator  and  denominator  of  the  rational  function.  To  begin, 
we  must  define  convergence  in  the  set  IRf s )  of  all  polynomials  over  iR.  Suppose 
P  :k=l  .2.  .  .  .  and  P  are  polynomials  in  IRfs]. 

K 

Definition  2._1  We  say  converges  to  P  if  there  exists  an  integer  q<=«  such 

that  deg  Pk<q:  k=1.2......  deg  PSq,  and  a^— *a^:i=0 . q.  where 


Vs1  '  Vs  ‘  "•  *  aiks 
p<s>  '  Vq'  "■  "  V  '  ao ' 


a0k: 


Remarks 

1)  If  we  regard  P^elR^s ] ;  k»l,2.***  as  functions  over  C.  we  might  be  tempted 

to  define  P  — »P  when  lim  P. (s)=P(s)  for  any  set.  But  we  notice  that  in 
k  k-»»  k 

this  definition,  deg  Pk  may  not  be  bounded.  For  example,  let 
1  k 

Pk(s)=^s  ■‘■1.  This  observation  brings  us  to  a  crossroads  in  the  theory: 

If  we  were  to  allow  convergent  polynomial  sequences  to  have  unbounded 
degree,  the  same  would  be  true  for  sequences  of  rational  functions.  This 
would  result  in  an  undesirable  situation  where  state-space  realizations 
could  have  unbounded  dimension.  Hence,  we  insist  on  bounded  degree  based 
both  on  physical  intuition  and  on  a  desire  for  mathematical  elegance. 

2)  Definition  2.1  is  equivalent  to  the  following  two  conditions: 

a)  {deg  p  |k*1.2. •••}  is  bounded. 

b)  lim  P.  (s)=P{s)  for  every  s  e  E. 

k-»  K 

Indeed,  necessity  of  a)  and  b)  is  obvious.  On  the  other  hand,  if  {P^} 
satisfies  a).  { s )  and  P(s)  can  be  written  as  in  Definition  2.1.  Choose 
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q*l  distinct,  complex  numbers  {s^-  .  ..s  Then,  from  b), 

V!k-V! 

where 


K  • 

••  si 

r 

V 

'a 

q 

••  s2 

i 

•  ?k 

alk 

,  t 

ai 

•ir 

Vi 

i. 

■a0k- 

■a0' 

We  know  that  the  Vandermonde  matrix  V  satisfies  det  V#0  as  long  as 

s.?*8..i#j:  therefore  V  1  exists,  and  $. — ►£  as  k — *». 

+  J  K 

We  can  define  a  topology  on  IR[s]  which  is  consistent  with  our  notion  of 

convergence  in  Definition  2.1.  To  do  so,  identify  every  element  in  IR[s) 

00 

with  an  element  in  (R  according  to 

v' . V  •  "o  •—  'VP,' ••••'«■ 

and  let 


\+l  3  t(p0,pl’  '  ‘  ‘  'V0,0,  ‘  ’  ’>6R  I  Pi6R:  1-0,1. 2.  ••••■»}. 


Then 


00 

a  »  ui 

k.i " 


is  the  set  of  all  polynomials.  On  we  take  identification  topology 
(e.g.,  see  [19, p. 120])  with  respect  to  the  bijections  f^iiR"—^  defined 
by 


fk^al ,a2’ ’ *  * ’*V  =  <ai -a2' ‘ * ’ -ak-0. • ’ * ) • 

That  is,  a  set 


U  =  Ha^ .  •  •  •  .a^.O.  •  •  • )  |  (a^ .  •  •  •  ,a^)€V} 

v 

is  open  in  if  and  only  if  V  is  open  in  (R  On  H  we  impose  the 
inductive  limit  topology  [19. p. 420]  with  respect  to  the  —  i.e.  we 
impose  on  $  the  finest  topology  which  makes  the  natural  imbeddings 
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continuous . 


It  is  routine  to  prove  that  *P  in  the  sense  of 
Definition  2.1  if  and  only  if  converges  to  P  in  'X. 

4)  It  is  shown  in  [6. Lemma  4.3]  that,  if  {P^}  is  convergent  in  tR r s | .  thou 
there  exist  convergent  real  sequences  {a..}.  and  {  k)  .  with  «  —0 

IK  IK  K  j  k 

and  limr,,*0,  such  that 

K 

Pk(s)  *  ^k^ttiks-1)^s~^ik) '  (3> 

i  i 

In  particular,  if  the  roots  of  P  are  bounded,  then  deg  P,  =  deg  limp 

k  k  k 

In  order  to  define  convergence  in  IR(s),  we  adopt  a  standard  quotient 
space  construction  over  IR(s]x(IR[s]-{0>  )  (e.g..  see  [18] .  p.136)  and  identify 
each  rational  function  with  a  unique  equivalence  class  under  the  relation 

(a.b)s(c.d)  «  ad=bc. 

We  use  the  expression  a/b  to  denote  both  a  rational  function  as  well  as  its 
corresponding  equivalence  class  in  IR[s]x(|R[s]-{0}  ) .  Note  that  a  similar 
construction  may  be  employed  in  identifying  rational  matrices  N/deRrm(s)  with 
equivalence  classes  of  pairs  (N, d)elR[s  ]rmx(IR[s  ] -{0} ) . 

Adopting  ordinary  quotient  set  topology  on  IR(s).  we  arrive  at  the 
following  definition. 

PRX.i HA t ion  2_. 2  Suppose  h^:k  =1,2 .  and  h  are  in  IR(s).  We  say  h^  converges 

to_h  in  IR(s).  if  there  exist  nR— *n  and  dk-*d  in  IR[s],  with  dk.d*0.  such  that 

‘V  t1k’11k:  k  =  l  1 2 .  aad  n/d"h. 

Along  similar  lines,  we  now  give  three  alternative  definitions  for 

rm 

convergence  in  !R(s) 

i i SALi^AD  .2^3  Suppose  Hk;k  =1,2 .  and  H  are  rxm  rational  matrices  with 

components  Ik  ^k  and  h^,  respectively.  We  say  converges  to  J  in  IR(s)rm  if 

*1ijk”"thij  in  IR(S*  as  k- 


7 


Definition  2  3’  Suppose  H^:k  =1.2,  •••  and  H  are  rxm  rational  matrices.  We 

say  H.  converges  to  H  in  (R!s)rm  if  there  exist  N. — k\'  in  IRfs]rm  and  d,  — *d  in 

IRfs]  such  that  N^/d^R^k  =1.2 .  and  N/d=H.  (Here  we  assume  product 

nn 

topology  on  (R(s|  and  that  the  quotient  space  constructions  above  are  applied 
componentwise  on  IRfs  | 1  mx(IR[s  J  -{0 } )  .  ) 

pm 

Definition  2.3 ' '  We  say  converges  to  H  in  (R(s)  if  there  exist  N^— »N  in 
IR(s)rm  and  Dk— »D  in  •R[s]m  with  D^  and  D  nonsinguiar  such  that 

W*’***'1-2 _  and  ND-1=H. 

Remarks 

1)  It  is  easy  to  show  that  Definitions  2.3.  2.3'  and  2.3' '  are  equivalent. 
A  fourth  alternative  definition  is  the  same  as  2.3''  except  using  left 
instead  of  right  factorizations. 

2)  Note  that  a  sequence  which  converges  in  the  sense  of  Definition  2.3  also 

converges  in  identification  topology  with  respect  to  the  map 
W:IR  — «IR(s)  defined  by 


3) 


KfE.A.B.C)  = 


OadjjsE  -A)B 
det(sE  -A)  ' 


where  (E.A,B.C)elRnl^nhrnTr* .  The  construction  of  the  topology  on  IR(s)tm 


shows  that  X  is  continuous. 

If  — »H  and  Gk— >G,  then  »H+G  and  H^— *HG;  more  generally,  IR(s) 


rm 


is  a  topological  ring  with  respect  to  identification  topology  on  IR(s)  and 

.  rm 


the  corresponding  product  topology  on  IRfs) 


In  particular,  relative 


rm 

topology  on  the  subgroup  of  polynomial  matrices  IRfs]  is  the  same  as 

rm 

product  topology  with  respect  to  Definition  2.1.  Note  that  IRfs]  is 
rm 

closed  in  IR(s)  . 

We  will  show  in  Section  5  that  our  definition  of  convergence  is  the 

rm 

"right"  definition  for  the  realization  problem,  since  a  sequence  in  IR(s) 
converges  in  our  sense  if  and  only  if  it  admits  a  convergent  sequence  of 
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state-space  realizations.  One  view  of  the  results  of  this  paper  is  that  they 
characterize  local  properties  of  the  map  H. 


3.  Time-Scale  Decomposition  of  Transfer  Matrix  Sequences 

Clearly,  any  rational  matrix  H  can  be  uniquely  expressed  as 

where  H  is  strictly  proper  and  H.  is  a  polynomial  matrix.  We  now  generalize 
s  r 

the  decomposition  to  the  sequential  case:  this  must  be  carried  out  in  a  way 
that  preserves  convergence. 

rm 

pe_fij]ition„  3. 1  1)  We  say  a  convergent  sequence  {H^}  in  IR(s)  is  a 

slow  sequence,  if  H  is  strictly  proper  for  every  k  and  there  exists  a  bounded 

-----  ^ 

region  Ac£  such  that  all  poles  of  each  H.  lie  in  A. 

2)  A  convergent  sequence  { >  is  called  a  fast  sequence  if  for  every 

M<»  there  exists  a  K<°°  such  that  k>K  implies  that  each  pole  p  of  satisfies 

p  >M  (all  poles  tend  to  infinity). 

Remarks 

Pm 

1)  The  set  of  all  slow  sequences  in  IR(s)  forms  a  proper  subspace  of  the 

rm 

real  vector  space  of  all  convergent  sequences  in  IR(s)  The  same 

statement  holds  for  fast  sequences. 

2)  Any  slow  sequence  can  be  expressed  as  Kgk  =  N^/d^  where  d^  is  convergent 

and  monic  for  every  k  and  deg  N.<deg  d.  ,  where  deg  N =max{deg  n..)  for  any 

k  k  itj 

polynomial  matrix  N .  Thus  deg  limNk<deg  lim  dk>  This  shows  that,  the 
limit  of  every  slow  sequence  is  strictly  proper. 

3)  Since  the  limit  of  any  fast  sequence  can  have  no  finite  poles,  such  a 
limit  must  be  a  polynomial  matrix. 

4)  Every  convergent  sequence  of  polynomial  matrices  is  a  fast  sequence. 
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3)  If  a  sequence  is  both  slow  and  fast,  it  must  be  strictly  and  have  no 
poles  whatsoever  for  large  k:  hence,  the  sequence  must  be  identically 
zero  for  large  k, 

6)  A  sequence  ot  matrices  {H^}  is  slow  (fast)  if  and  only  if  each  component 
sequence  {h^k}  is  slow  (fast). 

Definition  3.2  1)  We  say  Hk=Hsl<-Hfl<  is  a  time-  sraie.  decomposition  of  { Hk } 

when  <Hs^}  and  (H^)  are  slow  and  fast  sequences,  respectively. 

2)  In  a  time-scale  decomposition.  {Hgk)  and  (Hfk)  are  called  the  slow 
part  and  the  fast  part  of  { ) . 

Note  that,  from  remarks  2)  and  3)  above,  if  fIk~Hsk'k'ek  is  a  time-scale 
decomposition  of  { H.  > .  then  H  ,  — «H  ,  and  H..— *H - .  where  H,  and  are  the 
strictly  proper  part  and  the  polynomial  part  of  H=Jim  H^.  Theorem  3.3  tells 
us  that  every  convergent  sequence  {Hk)  has  an  essentially  unique  time-scale 


.  rm 


decomposition. 

Theorem  3.3  1)  For  every  convergent  sequence  {H^}  in  Rlsr"'.  there  exist  a 
slow  sequence  { }  and  a  fast  sequence  {HfkJ  such  that  Hk=Hsk'Hfk  f°r  every 
k. 


2)  If  {Hgkl  and  {Hfk>  are  slow  and  fast  sequences,  respectively,  and 


A  /V 


Hk=Hsk’Hfk  for  every  tflen  Hsk=Hsk  and  Hfk=Hfk  for  suff iciently  large  k. 
Proof  1)  We  need  only  treat  the  case  r=m=l:  the  multivariable  case  can  then 

be  handled  componentwise.  If  hk— *helR(s).  we  can  find  nk— *n  and  dk— *d.  with 

nk/dk=hk  and  n/,<**h,  Since  dk— *d,  from  we  can  writR  dk  =  rkdskdfk '  where 

)  =  s"  *  b..  ,  ,  sM_1  -  ...  *  b. 


sk’ 


Ac-1  ,k' 


Ok 


dfk(s)  '  auks  *  aw- 1 . ks  ■  alks  *  1 1 

with  each  {b^i  convergent,  rk — *T* 0,  and  a.k — K)  as  k — ♦».  Let 


Vs’  =  VsP  ‘ 


7*lkS  *  Z0k’ 
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and  let  <i 

l-maxl  v .  p  ij\  .  We  will  show  that 

there 

exist 

sequences 

n  .  -  x  ,  .  f  .  . . 

sk  fj-  1 .  k 

*  XlkS 

X0k 

q-l 

*  ylkS 

U-.  -=  v  .  ,  s  -  .  . . 

fk  q-l.k 

*  V0k 

such  that  ti,  (l-ii  .  d  .  +n*.  tL.  .  Equivalently,  we  need  to  show  that 
k  l<  sk  sk  fk  rk 

V"skdfk  -  nfkdsk>  "  V 

Note  that  equation  (4)  may  be  written  in  matrix  form 


written 

Alk 

i8lk 

■'42k 

*2k 

=  V 


where 


=  rk  ‘ 


V  l,k 


>l,k 


/i-i.kj 


q-1 . kj 


with  Alk>  A2k.  ^1R.  ^2k>  and  /-k  having  dimensions  uxij,  q x/u.  A^d.  d*q.  and 
(^q)xi,  respectively.  Also  note  that  Alk— • >rl  .  where  M  is 

nilpotent  and  upper  triangular,  and  — K)  as  k — *».  Hence,  there  exists  a  K<» 
such  that  (5)  has  a  unique  solution  when  k>K.  For  k£K.  let  hsk  be  any 
strictly  proper  rational  function  and  let 

2)  We  have  Hsk‘Hfk=Hsk*Hfk  for  sufficiently  large  k.  Hence, 

".kWHrk-  (SI 

But  the  left  side  of  (6)  is  slow,  and  the  right  side  is  fast.  Hence,  both 
sides  are  identically  zero  for  large  k.  J 
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I'o  conclude  this  section,  we  note  that  a  time  scale  decomposition  of  any 


transfer  matrix  sequence  of  the  form  H.  =C.  (sE.  -A.  )  8  .  .fhere  (IF.  .8.  .0.  )} 

KKKK  K  KkkK 

is  convergent  in  fRn^2n  m  and  detls-limE^  -limA^JSO.  can  be  achieved  by 
invoking  the  perturbat ional  form  of  the  classical  Weierstrass  decomposition 
for  matrix  pencils  as  developed  in  [7].  Indeed,  from  [7.p.l47j.  there  exist 
convergent  nonsingular  matrix  sequences  and  {N'^l.  with  limM^  and  limN'^ 

nonsingular,  such  that 


Wk 


I  0 
0  A 


fk 


•  Wk 


Ask  0 

0  I 


(7) 


where  limA^  is  nilpotent.  Let 


sk 


fk 


“  W  t<:3k  ctkJ  “  CA- 


(8) 


Then 


Vs'  *  cskISIAk,'leSk  - 


(9) 


From  Definition  3.1  it  is  clear  that 


Hsk(S)  -  Csk(Sl-Ask>  Bsk 


(10) 


and 


’Vs1  ■  ctklsAfk-irlfifk  1111 

are  slow  and  fast  sequences,  respectively.  Hence.  (9)  is  a  time-scale 
decomposition  of  { > . 

4 .  The  Characteristic  Polynomial  Sequence 


In  this  section  we  Investigate  several  useful  properties  of  the  sequence 

of  characteristic  polynomials  corresponding  to  a  convergent  sequence  {H^>  in 

IR(s)tm.  We  first  extend  the  conventional  definition  of  the  characteristic 

polynomial  to  improper  transfer  matrices.  Recall  that  the  characteristic 

polynomial  A  of  a  strictly  proper  rational  matrix  is  defined  as  the  least. 

common  monic  denominator  of  all  minors  of  H  . 

s 
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Definition  i  1  If  H  is  a  rational  matrix  with  H-li  -He  for  soma  strictly 
proper  H,  and  polynomial  matrix  H..  the  characters  Uc  polynomial  A  of  H  is 

S  I 

defined  as  the  characteristic  polynomial  of  H  . 

Consider  the  sequence  of  characteristic  polynomials  corresponding  to 

{ .  Since  !(k  N'k  d^,  it  follows  that  Ak  divides  ^i^rm)  for  e-,^  ^ .  thus, 
boundedness  of  {deg  d^>  ensures  boundedness  of  {deg  A^}.  Let  /7=lim{deg  A^l. 


and  note  that  deg  Sr)  for  sufficiently  large  k.  For  all  such  k,  A^  can  thus 
be  .uquely  identified  with  a  point  <A^>  in  the  real  projective  space  IP77  (see. 
e.g..  f 20 ] )  according  to 


i  i-1 

s  -a  * 


aQ  t-f  (0 . 0.1. a.  . aQ )  6  IR 


In  fact,  there  is  a  one-to-one  correspondence  between  IP  and  the  set  of  monic 


polynomials  A  with  degA<//.  These  observations  lead  to  the  following 


definition. 

Def  .4. 2  Let  be  any  convergent  sequence  in  IR(s)im,  and  let  Ak  be 

the  characteristic  polynomial  of  Hk-  Set 

(*<Ak> .  deg  Ar  Sr/ 

Pk  =  1 

L<1>.  deg  Ar  >rj 

The  sequence  { pk >  is  called  the  characteristic  polynomial  (CP)  of  { Hk 1 . 

It  is  easy  to  show  that  { p k>  converges  if  and  only  if  there  exists  a  real 
sequence  (rk>  such  that  {rkAk>  converges  to  a  non-zero  limit  AslR[s],  In  this 
case.  limpk=<A>.  We  now  present  several  pathological  situations  that  can 
arise  in  dealing  with  the  CP. 


Example  4.3 

1)  The  following  example  illustrates  that  when  {l!k}  is  convergent,  the 
corresponding  CP  may  not  converge.  Consider  the  sequence 
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2) 


3) 


Hk(s)  -  <! 


(s-2) 


(s~l)(s-2-£) 

(s-3) 


k  even 


k  odd 


l(si-l)(s-3^) 


and  let  H(s!- 


s  - 1 


We  may  write  Hk-=Nk  dk-  where  ( s -2 )  (s-3)  and 
f(  s  - 1 )  (s-2“)  (sr3) ,  k  even 


( s-1 )  ( $+2 )  ( s  ‘•3"k ) 


k  odd 


Thus  H.  — »H .  hut 
k 


Vs) 


’( s  *-1 )  ( s-2*j^) .  k  even 

_( s-1 )  (  s -3 )  .  k  odd 
K 


Clearly,  \p.)  Ls  not  convergent.  Note,  however,  that  <  H.  )  can  be  divided 


( 1) 


into  two  subsequences  with  convergent  CP's  according  to  Hk 


and 


H*:‘ 1  =  11  .  . 
k  2k 


In  some  cases,  { >  may  converge  even  though  (Ak>  does  not.  Consider  the 
rational  sequence 


Hk(s, 


(£s-1)(s-2) 


In  this  case,  {  Hk }  has  CP  determined  by  Ak{  s  )  =  (  s*-k)  (  s*2 ) ,  but  { > 
converges,,  since 

i  i 

s-2 . 


Jar(s)  =  i(s-k)(sr2) 


Finally,  we  note  that  convergence  of  { Pk }  (or  even  { Ak  1 )  does  not 
guarantee  that  lim  Ak  is  the  characteristic  polynomial  of  lim  Hk.  For 
example,  let 


Hk(s)  - 


s  *-2 


ls-l)(s-2*i) 


Then 


Ak(s] 


(s  +  l)(s*2rjL)  —  (8-1)18*2). 
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But 


Hk(s) 


1 

s-t 


Next  we  examine  some  basic  properties  of  the  CP  with  respect  to  the 
time-scale  decomposition.  First  we  need  a  simple  result  for  individual 
systems . 

Lemma  _4  .J  Suppose  H=H ^  -Hp .  where  and  have  no  common  poles,  and  let  1, 
A.  ,  and  A  be  the  characteristic  polynomials  of  H.  H  .  and  H  ,  respectively 

Id  I  & 

Then  A=A^A^ ■ 

Proof  From  the  definition  of  the  CP  we  can  assume  without  loss  of  generality 
that  and  are  strictly  proper.  Suppose  (A^.B  .Cj)  and  ( A2 .  B^,  *  C^ )  are 
minimal  realizations  of  and  :  then  \  ( s 1 -det ( s i -A . ) .  If  we  let 


,  C  -  [cl  c2J 


then  (A.B.C)  is  a  minimal  realization  of  H  with  CP 

A(s)  =  detlsI-Aj IdetlsI-Ag) . 

In  particular,  for  any  time-scale  decomposition.  Lemma  4.4  implies  that, 
when  k  is  sufficiently  large,  we  have 

4k  ■  4skV  "2) 

where  Ak>  Agk.  and  Afk  are  the  characteristic  polynomials  of  Hk>  Hsk>  and  H^, 
respectively. 


'A1 

0  1 

*2. 

A  = 

.  B  = 

.0 

A2- 

Lemma  4.5  Let  Hic=H3[<+Hfi<  be  a  time-scale  decomposition. 


1)  If  Agk  is  the  characteristic  polynomial  of  Hgk  and  AeR[s]  is  moriic. 
then  <Agk>— ><A>  if  and  only  if  Agk— *A. 

2)  If  is  the  CP  of  then  — ►<!>. 

3)  The  CP  of  { >  is  convergent  if  and  only  if  the  CP  of  {HgkJ  is 
convergent.  When  the  two  CP's  converge,  their  limits  coincide. 
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Proof  t)  Sufficiency  is  obvious.  To  show  necessity,  observe  that  there  must 

exist  a  real  sequence  { 7 k)  such  that  TkAsk~“,A  and  that  {Agk>  has  bounded 

roots.  From  (3).  {7^}  converges.  Since  Agk  and  A  are  monic. 

2)  We  have  0  =<^f i<  *  ■  where 

]T(s^ik) .  H^k  is  not.  a  polynomial  matrix. 

Afk  '  '  * 

1.  Hfk  is  a  polynomial  matrix. 

Here  the  A.,  satisfy  the  property  that,  for  every  M<»,  there  exists  a  K<~  such 

1 K 

that  >M  for  each  i  and  each  k>K.  I.et 

ik 

fiJr:  is  not  a  polynomial  matrix, 

i  ik  * 


1,  Hfk  is  a  polynomial  matrix.- 


Then  rkAfk— ♦  !  in  IRfs) . 


3)  The  result  follows  immediately  from  (12)  and  part  2).  2 

The  final  result  of  this  section  focuses  on  the  observation  made  in 

rnt 

Example  4.3.  part  1).  that  the  CP  of  a  sequence  (H^)  in  IR  (s)  which  is  not 
convergent  can  sometimes  be  decomposed  into  convergent  subsequences.  We  can 
in  fact  demonstrate  that  a  finite  decomposition  of  this  sort  can  always  be 
achieved . 


Theorem  4^6  If  Hk— ♦H,  then  {pk)  has  finitely  many  limit  points. 

Proof  From  Lemma  4.5,  part  3).  we  need  only  consider  the  casp  where  (H.  )  is  a 

K 

slow  sequence.  Let  Hk=Nk/dk  and  H=N/d.  From  the  definition  of  the 
characteristic  polynomial.  Ak  divides  droin{r,m}  for  each  k  But 
^minfr.m} — ^min{r,m}  g0  tjie  uniqUe  monic  representative  of  each  limit  point 
of  (pk)  must  divide  The  result  then  follows  from  the  fact  that  any 

polynomial  over  IR  has  finitely  many  monic  divisors. 
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Corollary  4  J  It  H^— *H.  then  there  exist  f  initely  many  strictly  increasing 
sequences  (kj>:i  =  l . Jt  of  positive  integers  such  that. 

1)  {kj  |  i=l ,  2 . Jt;  j-1.2.  .  .  1  -  {1.2.3..  .  . 

2)  (kj  |  j  =  l .  2  .  .  .  \  fl  {k^i  j  =  l  .2.  . . . }  =  0  when  p*q. 

3)  each  {H  . }  has  convergent  CP. 

k5 

Proof  From  Theorem  4.6.  there  are  only  finitely  many  limit  points 

p1 . p^elp’1  of  { ,o^ > .  Since  IP^  is  a  compact  Hausdorff  space,  each  open  subset 

U  of  IP51  satisfying  {p* . p^KU  contains  a  tail  of  { .o^ } .  Indeed,  otherwise 

there  would  exist  a  subsequence  of  Ip  I  with  no  limit  point,  contradicting 

K 

compactness  of  IP7'.  Let  U, . U  be  nonintersecting  neighborhoods  of 

p* . p* .  respectively:  then  there  exists  a  K<“  such  that  {p^JcuUj  for  k>K. 

Let  k j  =  j :  j  =  l . K.  The  remaining  k^  may  then  be  defined  iteratively 

according  to  kj=min(  {k  j p^eU . }- {k*  j q<j  } ) .  If  V^cL'^  is  another  neighborhood  of 

p.,  then  by  compactness  of  IP1"*  there  must  be  a  tail  of  the  subsequence  {p  .} 
1  k1 

“j 

contained  in  V,;,  hence,  p  »p  .  * 

1  ki 


3.  Existence  of  Realizations 

We  are  now  ready  to  formally  define  realizations  of  a  given  transfer 
matrix  sequence  and  discuss  their  existence.  We  base  our  definition  of  a 
realization  of  {Hk)  on  the  standard  definition  of  a  realization  of  a  single 
rational  matrix  H  as  in  Theorem  1.2.. 
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mt 

Definition  5.1  1)  .Suppose  {H.  )  converges  in  IR(s)  .  We  say  a  sequence  {c/,.l 

K  K 

in  IRll^n^m  1  *  is  a  realization  of  {H^}.  if  there  exists  an  integer  K  and  a 

<7€lRnln*m  +  r)  such  that  a k  is  a  realization  of  when  k>K  and  <;k— *o  in 
^nian-mn' ) 

2)  A  realization  {o^}  in  (pn*nrm'hr*  said  to  have  dimension  n. 

Note  that  the  dimension  of  a  realization  {o^l  is  given  simply  by  dim 
for  any  k.  if  Hk~ *H.  then  continuity  of  U  implies  that  a  is  a  realization  of 
H  in  the  conventional  sense.  >ve  will  show  that  there  exists  a  realization  for 
any  convergent  sequence  { > :  this  generalizes  part  1)  of  Theorem  1.2  to 
sequences  and  demonstrates  that  the  definition  of  convergence  in  IR(s)im 
outlined  iri  Section  2  is  'he  correct  one  for  our  purposes. 

To  simplify  subsequent  discussion,  we  will  make  use  of  the  mapping 
'3 :IR( s ) rm — HR(s)rni  defined  by 

$(H)(s)  =  4H4>- 

o  o 

It  is  easy  to  see  that  i  is  an  isomorphism  on  IR(s)rm  and  that  <1  *='$.  Some 
elementary  properties  of  i  follow. 

Lemma  5.2  Let  HelR(s)im,  and  let  {H^}  be  convergent  in  (R(s)rm. 

1)  (E.A.B.C)  is  a  realization  of  H  if  and  only  if  (A.K.B.C)  is  a  realization 
of  ■*» (  H )  . 

2) 

3}  If  the  characteristic  polynomial  of  H  is  A(s)=sn-)7^_1sn  1  +  ...<-i70.  then 
the  characteristic  polynomial  of  'S(H)  is  c(s)=r(>?0sn- .  . .  jS*  1 )  for 
some  r*0. 

Proof  1)  Suppose  H { s ) -C ( sE -A )  Then 

i(H)(s)  =  -ic{iE-A)_1B  =  C(sA-E) _1B . 
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2)  From  !).  if  H  lias  a  realization  of  degree  n.  then  so  does  'i(ll).  The 
converse  follows  from  S(1{H))=H. 

3}  Let  (E.A.B.C)  be  a  minimal  realization  of  H;  then  (A.K.M.C)  is  a 


minimal  realization  of  *1(11)  Hence.  A(  s )  =  r jdet (sE- A)  for  some  r^O  (for 


details,  see  |14|).  and 


r(s)  «  r2det(sA-E) 

=  r2< -s)ndet(~R -A) 

r2.  .n-.l-n  . 1 . n-1 

•  -(-SI  Mjl  -Vl(j)  -..-*^1 


Lemma  5.3  If  {H^l  is  a  fast  sequence,  then  is  a  slow  sequence. 

Proof  Since  all  poles  of  tend  to  infinity,  we  can  write 


V8)  n  F 


NVs, 


Tf(«iks-1) 
l-l  1K 

where  { >  is  convergent  and  each  a^— -*0  as  k- 


Let  q=max{deg  n  ,p-l}, 
k  1JK 


p-3  S\(s} 

=  —i rr—  • 

«  K  H(s-«.k) 

Note  that  s^N^I^-)  is  a  polynomial  matrix,  each  of  whose  elements  has  degree  at 
most.  q.  Clearly,  '5 ( )  has  bounded  poles,  and  the  denominator  of  )  has 
degree  q+1,  so  ^(H^)  is  strictly  proper.  " 

Now  consider  a  time -scale  decomposition  Nj^sk^fk  of  an  arbitrary  {Hkl 


in  IR(s) 


Suppose  and  {’KH^ll  have  realizations  of  the  form 


{(I*Ask*B.k*Csk’»  and  {(I,Afk'Bfk’('fk,}'  Then  eaoh  has 

< Afk •  I  * Bf ^ ’ C^. k )  as  one  of  its  realizations.  Defining 
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f  «  ! 


°  Afk 


Ask  0 
0  I 


Ck  "  [Csk  !'fk! 


it  is  easy  to  check  that  { ( £R . AR , 8R , ) }  is  a  realization  of  {Hk>.  Therefore, 
we  only  need  to  prove  existence  of  realizations  for  slow  sequences. 

Theorem  5.4  Every  slow  sequence  has  a  realization  of  the  form  {( l  ,.-\R.BR.i 'R) !  . 
Proof  First  we  treat  the  case  r=m=l.  Let  H.=N’R/dR.  where  {N'R}  and  {dR|  are 


convergent  in  IRfs).  Then 


ci  1  .  k 


sql.- 


Nk  ’  ^q-l.ks  ^q-2.kS 


Ok 

«*“2 . , 


where  a.,  and  p.,  converge  as  k — >®.  To  obtain  a  realization  of  (H,  }  of  the 

ik  ik  K 

desired  form,  set 


•  Bk  =  0  •  Ck  =  [/}0k  ^lk 


Vl.k^ ' 


L  ~0k  '"Ik  q -1 . kJ 

Now  we  consider  the  general  case.  By  the  definition  of  convergence  in 


IR(s)  .  every  component  sequence  {h^R}  convergent. 

{( I  ,Ar’  .Br^  .CrJ  ) }  is  a  realization  of  {h^}.  Let 


AR=diag{ARj|j-1.2 . m>.  8R=diag{BRj|  j-1.2.  •  •  •  .m}.  Cj^tC*1 


Suppose 


AR  =  diag{AR| 1  =  1.2 . r > ,  BR  =  •  ,  CR=diag{Ck  i-1 .2.  •  •  •  ,r) . 


A  simple  calculation  verifies  that  { ( I .A^.B^.C^) }  is  a  reaiization  of  { > .  . 

Combining  the  time-scale  decomposition  with  Lemma  5.3  and  Theorem  5.4.  we 
arrive  at  the  following  result. 


i  ill 

Corollary  5.5  Every  convergent  sequence  in  (Rfs)  has  a  realization. 
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X 


A 


Theorem  5.4  (but  not  Corollary  5.5)  may  also  be  proven  in  an  abstrai  “. 

algebraic  framework  as  outlined  in  117. Chapter  4].  Briefly,  consider  the 

commutative  ring  c  of  convergent  sequences  in  IR  using  pointwise  operations. 

and  let  the  set  of  r^m  proper  rational  real  matrices  be  denoted  by  IR  (s),m.  A 

P 

rrn 

convergent  sequence  (H.  }  in  IR  (s)  may  then  be  viewed  as  a  formal  power 

k  p 

series  over  the  ring  of  r*m  matrices  with  elements  in  c.  Indeed,  we  may 
expand  each  element  of  each  about  s=<*>.  yielding  the  series 

Hk  ■  031 

where  the  sequences  (H..  >  are  convergent,  From  this  point  there  are  two  ways 

1 K 

to  proceed.  First,  one  can  prove  realizability  by  constructing  a  certain 

infinite-dimensional  Hankel  matrix  from  the  H..  .  It  must  then  be  shown  that 

lk 

the  span  of  the  columns  of  the  Hankel  matrix  is  a  finitely  generated  module 
over  c.  A  second  approach  is  to  show  that  the  formal  power  series  (13)  is 
’rational"  in  a  certain  algebraic  sense.  This  immediately  guarantees 
realizability.  Both  conditions  can  be  demonstrated  in  our  framework  fairly 
easily;  however,  our  proof  of  Theorem  5.4  is  more  direct  and  is  sufficient  for 


our  purposes. 


6.  Minimality 

rtn 

In  Section  5  we  showed  that  every  convergent  sequence  {H^}  in  R(s)  has 
a  convergent  realization  {^1-  Ih  this  section,  we  explore  the  issue  of 
minimality  of  a  realization. 

Definition  6.1  1)  If  n  is  the  smallest  integer  such  that  {H^}  has  a 

realization  of  dimension  n.  and  is  a  realization  of  (Hk)  with  dimok=n. 

then  we  say  {a,  }  is  a  minimal  realization  of  {H.  }. 

K  *  j\ 
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*  t  , 


2)  If  a  sequence  of  state-space  systems  {rr  I  is  a  minimal  realization  of 

K 

its  transfer  matrix  sequence,  we  say  <crk>  Is  minimal. 

Obviously.,  all  minimal  realizations  of  {H^}  have  the  same  dimension. 
This  fact  enables  us  to  define  a  degree  function  5  on  the  set  of  convergent 
rational  matrix  sequences  by  setting  SUI^}  equal  to  the  dimension  of  any 
minimal  realization  of  IH^J.  In  this  section  we  develop  a  simple  expression 
for  5iH^}  for  slow  sequences  and  then  extend  it  to  the  general  case.  Next,  we 
examine  a  natural  conjecture  for  determining  whether  a  sequence  < o- k >  is 
minimal  and  relate  minimal  realizations  of  the  same  (Hk)  in  a  manner  analogous 
to  Theorem  1.2.  part  4).  finally,  we  relate  our  results  to  the  realization 
theory  outlined  in  [17]  for  algebraic  systems  over  the  ring  c. 

In  our  development  it  will  be  helpful  to  exploit  various  properties  of 
the  mapping  which  takes  each  state-space  system  into  a  particular  choice  of 


numerator  and  denominator  of  its  transfer  function.  Specifically,  define 
i n : Ik  — HR  according  to 

r  ( A. B . C )  =  (C-adj(sJ-A)B.det(sI-A)). 

Here  we  have  identified  as  defined  in  Remark  3)  after  Definition  2.1.  with 

IR  .  Note  that  f^  is  continuous:  if  rn(A,B,C)  =  (N,d) ,  then  (I.A.B.C)  is  a 

realization  of  N/d,  Also  notice  the  distinction  between  1’  and  H,  as  defined 

n 

in  Section  2.  We  denote  Im  f  =1’  (IRn^ n+n,fct  1 ) . 

n  n 

The  following  series  of  Jemmas  leads  us  to  the  first  main  theorem  of  this 


section. 


LSMiiLJS.Jt  Consider  any  pair  (N.d)  where  d  is  monic,  degd^n,  and  N/d  is 
strictly  proper  with  characteristic  polynomial  A.  Then  (N,d)#Imr  if  and 
only  if  A  divides  d. 
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Proof  (Sufficient)  Suppose 


d(s)  =  A(s )  ff  ( s~/3, ) , 
i  =  l 

and  let  (I.A.B.C)  be  a  minimal  realization  of  N/d.  Define 


where 


.  C  =  [C  OJ. 


det(sI-A  )  =  det(sI-A)det(sI~i:}  -  Ms)  It  (sr/j.)  <  d(s). 

i=l 

Since  C  (sI-A  )  B  =C(sI-A)  B=N(s)/d(s) , 

*  *  *  *  M(s) 

C  adj(sI-A  )B  =  det(sI-A  =  N(s). 

d(s) 

9(B  9|C  9|C 

Hence,  (N.d)»r  (A  ,B  ,C  ). 

(Necessary)  Suppose  (N,d)»r  (A.B.C) ;  then  (I.A.B.C)  is  a  realization  of 
N/d.  From  (16, Theorem  5-181,  we  can  find  a  similarity  transformation  T  such 


TA  A  A  1  TB  1 

11  12  "13  1 

T_1AT  =  0  A22  A23  ,  T_1B  =  Bg  ,  CT  -  [0  1'2  J . 

0  0  A„„  0 


where  (A„0.B0,C.)  is  a  minimal  realization  of  N/d.  Note  that 

det(sI-A22)»A(s) .  Thus 

d(s)  =  det(sI-A)  =  dettsI-AjjJdetlsI-A  3)A(s). 

From  Corollary  4.7,  (H^)  can  be  decomposed  into  n  sequences 

{H*^);i  =  1.2 . it,  where  we  define  H^=H  ..  Each  sequence  has  convergent  CP 

K  J  j.1 

kJ 


satisfying 


<A*i}>  <A(1)>:  i-1.2 . n, 
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*  «! 


WllPt'e  \ 


( I ) 


( i>  t(i) 


is  monie. 
Let 


if  is  slow.  then,  from  Lemma  4.o,  part  1). 


A 


LCM 


where  I.i’M  denotes  the  least  common  multiple.  Also  define 


A,.1  -  A 


(DA 


k  tn 


:  i  - 1 . 2 , 


,  K. 


( u: 


Mote  that  each  A^ '  is  a  polynomial  and  that,  if  {H^}  is  slow.  A^' ' — »A. 

Lejma  6.3  Let  {H^}  be  a  slow  sequence  with  »H,  and  suppose  H  has 
characteristic  polynomial  A.  Then  A  divides  A. 

Proof  According  to  Corollary  4.7.  {H.J  can  be  decomposed  into  n  subsequences 

“  *  K 


with  convergent  CP’s  A^1— *A(  1 


[f  A  divides  A’*’  for  each  i.  then  A  divides 
A.  Hence,  it  suffices  to  treat  the  case  where  {H^i  has  convergent  CP  A^ — *A. 


r  in 

Let  p-=min{r .m} .  and  consider  the  sequence  { H.  }  of  lx  £  (.)(,)  rational 

K  i*l  1 


matrices,  each  Hk  consisting  of  the  minors  of  of  all  orders.  Obviously. 
%  %  * 

H. — >H.  where  H  is  defined  similarly.  It  follows  from  elementary  arguments 

K 

that  has  characteristic  polynomial  A^  (same  as  H^)  and  that,  for  any 

polynomial  q.  qH^  is  a  polynomial  matrix  if  and  only  if  A^  divides  q.  in 
particular.  AkHk  is  a  polynomial  matrix..  Since  iSfs]^  is  closed  in  IR(s)"1  for 
any  j,  AH  is  a  polynomial  matrix.  Thus,  the  characteristic  polynomial  A  of  H 
(and  H)  divides  A.  !■ 

Lemma  6.4  Pn  is  an  open  mapping. 

Proof  Note  that  T  is  multilinear;  thus,  it  is  a  composition  of  functions  on 
-  n 

Euclidean  spaces  IRP  of  the  form  f  ( x ^ . Xp’  *-xixj  anci  . xp)sxl*  ‘ ' '  Xp‘ 

Since  f  and  g  are  open,  compositions  of  open  maps  are  open,  and  products  of 
opens  sets  are  open,  it  follows  that  l'  is  open.  “ 
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Lemma  6.5  Let  X  and  Y  be  topological  spaces  with  X  first  countable,  and  let 
Q : X — »Y  be  an  onto.  open,  and  continuous  map.  For  any  convergent  sequence  (y^} 
in  V  with  y^— »yeY  and  any  xeQ  1 { y ) .  there  exist  x^eQ  1 ( ) :k-t .2. . . .  such  that 
xk-x. 

Proof  Let  {t.\  :  i;  1 ,2.  •  •  • }  be  a  countable  basis  of  neighborhoods  of  x  with 

U .  i  t  .  Since  Q  is  open,  each  V.=Q(l\)  is  a  neighborhood  of  y.  Hence,  for 

any  V..  we  can  find  an  integer  K.  such  that  y.  eV.  when  k>K . .  Furthermore. 
i  x  k  l  i 


there  must  exist  points  x/^eU.  :k*K.  *-1  .K.  *2, . . .  with  Q(x^l*)-y,. 

Kill  K  K 


For  k<K  . 


select  any  x^i!e  Q  l(y  ).  This  process  defines  sequences  {x! 1  * l : i«l , 2 . . . . 
n  ci  k 

Without  loss  of  generality.  we  may  assume  K^^K..  If  we  let 

x=x/l,;k=K.  where  Krt“0.  the  construction  shows  that  each  U, 

k  k  1-1  i  0  i 

contains  a  tail  of  the  sequence  {x^}.  Hence,  x^ — »x.  2 

Lemma,  .6 ..6  Suppose  { H^ >  is  a  slow  sequence  with  H^-h >H.  If  there  are  pairs 
(N’k.dk).  (N.d)elm  rn  such  that  Nk/dk=Hk.  N'/d=H.  and  ( Nk . dR )— *( N . d ) ,  then  {HR} 
has  an  n -dimensional  realization. 

Proof  Note  that  (pn(nfmft  ^  is  first  countable.  Thus,  if  we  restrict  the  range 


of  P  to  Im  r  ,  we  may  use  Lemmas  6.4  and  6.5  and  the  fact  that  P  is 
n  n  n 


continuous  to  conclude  that  there  exists 

0n(n+m+r ) 


convergent  sequence 


(Ak,Bk.C*k) — ►  (A.B.C)  in 


such  that  r[i(Ak,Bk,Ck)  =  (Nk.dk)  and 


I*  (A.B.CMN.d)  • 


Notice  that  ( I , )  is 


realization  of 


Nk/dk=Hk;k=l ,2 .  and  (I, A.B.C)  is  a  realization  of  N/d=H . 


l.em!na__6^_7  If  a  slow  sequence  { Hk >  has  an  n-dimensional  realization,  then  it 
has  ..n  n-dimensional  realization  of  the  form  { ( I , Ak , ) > , 

Proof  Suppose  {Hk>  has  a  realization  having  dimension  n.  Since  { >  is  slow, 
the  decomposition  ( 7 ) - ( 1 1 )  shows  that  (Hkl  is  of  the  form 

Vs1  •  <WSI-Askl'1V 
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c 


where  is  q-q  with  q<n.  Define 

TA  .  01  TH 

V  -  .  B„  ■  ! 

k  [  0  oj  k  [  < 

Theore_m_6^8  For  any  slow  sequence  { H.  > . 


ck  .  [<:sk  O] 


51  Hk)  -•  Ill's 

Proof  Let  no=deg  A  .  We  first  show  that  there  exists  an  no-dimensional 
realization  of  {H^}.  Let  and  M=AH:  then  and  N  are 

polynomial  matrices  with  (N^  *  .A^ 1 )— *(.\*. A) .  Since  all  poles  of  are 
bounded,  Remark  4)  after  Definition  2.1  shows  that  deg  A^'=deg  A=*n0.  Thus, 
from  (14)  and  Lemmas  6.2  and  6.3.  (N^ * .A^* * ) . (N, A)elm  T  .  Suppose  { , A^ ) > 

is  constructed  by  setting  N  =N’  ^  and  A  -A^*  whenever  H  where  the 

k*  J  k1  J  k1  J 

j  J  J 

kj  are  defined  by  the  decomposition  of  Corollary  4.7.  Then  (N^.A^jslm  . 

The  desired  result  follows  from  Lemma  6.6, 

It  remains  to  show  that  no  is  the  minimal  dimension  over  all  realizations 

of  { > .  Suppose  {Hk>  has  an  n-dimerisional  realization.  Then,  from 

Lemma  6.7.  it  has  an  n  dimensional  realization  of  the  form  { ( I  .A^.B^.C^M  . 

Let  ( N*r . dR ) ^ Tn ( . Bk . Ck )  and  {N,d)=rn<  l im ( AR . BR . CR ) ) :  from  Lemma  6.2.  AR 

divides  d.  for  every  k.  Letting  d^*=d  .,  a/^  divides  d/  ^ .  Since  I'  is 

K  J  ,  1  K  K  n 


m  .  |  — it 

k1  k 
J 


continuous,  d^— -*d;  thus,  closure  of  IR[sJc(R(s)  guarantees  that  each  A*1* 
divides  d.  Thus  A  divides  d,  and 

n  =  deg  d  >  deg  A  =  n o . 

The  following  result  offers  one  method  of  calculating  <5{Hk>  for  an 
arbitrary  convergent  sequence  {H^}  in  (R(s)rm. 

Theorem  6.9  If  is  a  time-scale  decomposition,  then 


«<Hk)  .  5<Ksk>  *  5(Hfk). 
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Proof 

Suppose 

( 7 ) - ( 11 ) .  it 

'V 

< Afk ■ 1 ' 8 

there 

p.xists 

dim  a  , 
sk 

<dim  crgk 

Ek  = 

^fk^^fk' 1  'Bfk’Cfk^  are  real izations  of  {Hgk)  and  { H f > .  Suppose 
there  p.xists  a  realization  ^sk^^sk'^sk’^sk'^sk*^  ^*sk^  wittl 
, ini  <dim  ogk  .  By  Lemma  6.7,  we  may  assume  that  Egk=I,  Let 


I  0 
0  A 


fA  .  01 
sk 

0  I 


•  c.  -  1C  ,  C-.  ]  . 

k  sk  fk1 


Then  { crk >  =  { (Lk.Ak.Bk,Ck)  )  is  a  realization  of  {H^}  with  dim  ok  <dim  ok  .  This 
is  a  contradiction.  A  similar  argument  shows  minimality  of  {o^}.  I 

Thus,  one  way  to  find  <5  { >  is  to  first  perform  a  time-scale 
decomposition  and  then  to  use  Theorem  6.8  to  find  <5{Hgk>  and 

5{Hfk}-d{'^(Hfk)} •  Fortunately,  the  next  theorem  simplifies  this  task  and 
shows  how  to  calculate  <5{Hk)  without  resorting  to  time-scale  decomposition. 

rm 

Recall  that,  for  any  He!R(s)  .  Hf  denotes  the  polynomial  part  of  H. 

Theorem  6.10  Suppose  Hk— *H.  Then 

«{Hk)  *  max  Tim  (deg  A*i}  *•  (H^1 1  )f ) ) . 

i  k  ‘  a 

Proof  Suppose  Hk=Hgk-H£k  *s  a  time-scale  decomposition  of  { > ,  and  let 

„(!}.„  H(i) kH 

•  "  »  **n  »  i  n  *  .  *’i»  *  a 

Skj  fkj 

[t  is  clear  that  "Hsk^*^fk^  is  a*so  a  time-scale  decomposition  and  that 
HgkJ  and  have  characteristic  polynomials  Ag*!  and  A^* ,  respectively. 


From  (12), 


*li]  - 

%  sk  fk 


From  (15)  and  Lemma  4.5.  part  2),  Agk*— »A^ 
5{Hgk)=deg  A.  Also,  since  each  Agk* 
deg  Agk*=deg  A(i>  for  large  k. 


Hence,  from  Theorem  6.8, 
is  convergent  and  monic. 
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4 


t 


*’  »  "-i 


From  (3).  --A...  s  is  of  the  form 

i  K 


<4fk>  ■  'nii«,k»-n>  • 


where  sik— *0  as  *».  Note  that  some  of  the  «  may  vanish:  so.  from 
lemma  5.2.  part  3).,  has  characteristic:  polynomial  of  the  form 

:k(s)  =  sq«lks‘ 

where  q  may  depend  on  k.  Therefore,  the  limit  of  any  convergent  subsequence 
of  {r^}  is  of  the  form  s<l .  Suppose  r  is  the  least  common  multiple  of  these 


limits:  then  r^s1*.  where 

q  =  lim  deg(?k)  =  lim  /uH( Hf ^ ) )  =  Tim  A*(Hfk). 
k  k  k 

The  last  equality  is  obtained  from  Lemma  5.2.  part  2).  Arguing  as  in 
Lemma  5.2,  part  2).  <5H(Hk)  >  =«5 ( >  for  any  { Hk > .  Hence,  from  Theorem  6.8.  we 
have 


«{Hfk>  =  «{'1(Hfk)}  =  deg  r  =  q. 

Theorem  1.2.  part  2)  and  Theorem  6.8  show  that  /n(Hfk)  =deg  Afk  +(j(  ( ) f ) . 
Theorem  6.9. 


From 


5{Hki  =  deg  A  *■  lim  (deg  Afk  +  ( Hk ) f ) » 

=  Urn  (deg  A  +  deg  Afk  +  ^<  ( Hk ) f ) ) 
k 

*  max  lim  (deg  A  *  deg  A^J  *  u(  (  1 ) f ) ) . 


It  remains  to  prove 

deg  Ak**  =  deg  A  *•  deg  A^.  (16) 

By  the  definition  of  Ak^  ., 

deg  Akj)  =  deg  A  +  deg  A^1*  -  deg  A(il.  (IT) 

Since  deg  A^k*=deg  A***,  it  follows  from  (15)  that 

deg  A^11  =  deg  A(i>  ►  deg  AfkU.  (18) 

Combining  (17)  and  (18).  we  obtain  (16).  ■“ 
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Corot  Ut'i  i>  11  Suppose  >H.  If  the  CP  of  is  convergent,  then 


0<Hk>  =  lim  #/(Hk). 
k  * 

prooj  In  this  case.  7t=l  and  so 

<5 { }  =  Tim  (deg  AR  p((Hk)f)l 
k 

-  TTm  M\)a  -  uU Hk)f )  1 

k  a 

•  TTm  u(H.  ). 
k  K 

Our  next  goal  is  to  generalize  part  3)  of  Theorem  1.2.  An  obvious 
conjecture  is  that  a  realization  (<rk)  of  (Hk)  is  minimal  if  and  only  if  each 
ok  is  controllable  and  observable  (as  defined  in  [8]).  While  controllability 
and  observability  for  every  k  (or.  indeed,  for  infinitely  many  k)  are  clearly 
sufficient  for  minimality,  the  next  examples  demonstrate  how  necessity  can 
fail. 


Example  .6.  ! 2  1)  Even  for  a  slow  lHk>.  minimality  of  { cr^  1  does  not  imply 
controllability  and  observability  even  at  a  single  point.  Consider  the 
sequence  { Hk >  in  Example  4.3,  part  1)..  We  can  decompose  {Hk}  into  two 
subsequences 


„(1)  _  s+2 

( s  1 )  (s»2*g^"M^  ) 

h(2)  .  s-*-3 

"w  i 

(s+l)(s*3*^) 

For  Hk^  and  Hk2*,  the  CP's  are 

A^!)  =  (s^Ms^+gj-y)  —  (S+-1)  (s-t-2)  -  A(t) 

A<2)  «  (S‘l )  ( 8+3+2^)  —  (s*l)  (s*3)  -•  Ai2). 

Thus 

A(s)  =  LCM  {A(1),A(2)}  =  (s+1 )  (s+2)  (s*-3) . 

Since  (Hk>  is  a  slow  sequence,  Theorem  6.9  indicates  that  5{Hk)=deg  A  =3 . 


•  l-  > 


Every  minimal  realization  {o^l  of  { must  have  dimension  3.  but  u{o^)=2  for 

each  k;  hence,  no  cr^  is  controllable  and  observable. 

2)  In  this  example,  the  CP  converges,  but  controllability  and 

observability  on  a  subsequence  is  the  most  that  can  be  achieved,  bet 

'1 .  k  even 

Vs>  -  — .  k  odd  ' 

l(is-n2 

A  simple  calculation  yields  d{H^>=2,  but  ^(^1  =  1  when  k  is  even,  so  any 
realization  must  contain  infinitely  many  terms  which  are  not  controllable  and 
observable. 

The  next  result  offers  a  weak  extension  of  Theorem  1.2.  part  3)  to  the 
sequential  case. 

Theorem  ,6.13  Consider  a  convergent  sequence  {cr.  )*{(E^.A.  . . C^ ) ) ,  and  let 

HksCk^sKk~Ak^”lBk‘  C:P  of  c:onver&es  and  is  minimal,  then 

there  exists  a  subsequence  (u  }  such  that  <t.  is  controllable  and  observable 

kj  kj 

for  every  j.  If.  in  addition.  {H^l  is  a  slow  sequence,  then  lim  is 

controllable  and  observable,  so  is  controllable  and  observable  for  every 
sufficiently  large  k. 

Proof  From  Corollary  6.11,  there  exists  a  subsequence  (Hk  }  of  { >  such  that 
${H.  }=ju(H.  )  for  all  j.  Therefore,  each  a  has  dimension  u(ti  )  and  must  be 

k  kj  k.  kj 

controllable  and  observable.  If  {Hk}  is  a  slow  sequence  and  the  CP  converges, 
Lemma  4.5,  part  1)  shows  that  Ak~ *A,  where  A  is  the  characteristic  polynomial 
of  H=limHk>  Since  each  Ak  is  monic,  for  large  k  we  have 

/^(Hk)=deg  Ak  =deg  A  =^(H) .  From  Corollary  6.11,  dim  a  ««5{ Hk Hk )  i  hence,  n  is 

controllable  and  observable.  j 

Restricting  attention  to  slow  sequences.  Example  6.12,  part  1)  has 

special  significance  from  an  abstract  algebraic  perspective.  It  is  easy  to 
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show  that  algebraic'  controllability  and  observability  over  the  ring  e  of 
convergent  real  sequences,  as  defined  in  [17. Chapter  21.  Is  equivalent  t,o 
controllability  and  observability  of  cr^  for  sufficiently  large  k.  Thus, 
linear  systems  over  c  do  not  satisfy  the  property  that  minimality  implies 
algebraic  controllability  and  observability. 

We  conclude  this  section  by  examining  the  problem  of  extending 
Theorem  1.2,  part  4)  to  the  sequential  case.  The  following  examples  show 
that,  in  our  case,  two  minimal  realizations  of  the  same  sequence  {H^}  may  not 
be  related  by  nonsingular  transformations  (cf.  [17, Theorem  4.19]). 

Example  6,14  1)  In  fact,  two  minimal  realizations  may  not  be  related  by 

nonsingular  transformation  for  any  value  of  k.  To  see  this,  let  {H^}  again  be 
the  slow  sequence  given  in  Example  4.3,  part  1).  In  Example  6.12,  part  1 ) .  we 
showed  that  <5{ Hk }^=3 .  Consider  the  two  minimal  realizations  { ( I } 
and  {(I,A21c.Bk.C2k)}.  where 


‘lk 


0 

0 


1 

0 


■(8*S)  -(11^) 


0 

0 


1 

0 


0 

1 

«4’ 

0 

1 


-C«*g)  -(6*£)J 


2k 


0 

0 


1 

0 


5  1 


7  1 


k  even 


k  odd 


0 

1 


i-<6w>  -<iiw)  -(6+eJ 


‘O' 

[6+! 

4 

1 ] ,  k  even 

8k  - 

0 

1 

•  cik  ■  5  »■  c2k  .. 

H 

4 

1 j ,  k  odd 

Suppose  there  exist  convergent  nonsingular  matrix  sequences  {M^}  and  {N^}  such 
that 


31 


Vk 


r'  MkAlkN lk  J  A2k 


Mk8ik  ’  82k'  ClkNk  =  C2k 

for  each  k.  Then  and  ^2k=^kAlk^k1 '  a  simple?  -calculation  shows 
that  A  k  and  A2k  have  different  spectra,  yielding  a  contradiction  for  each 
value  of  k. 

2)  When  the  CP  converges.  Theorem  6.13  implies  that  there  exists  a 
subsequence  on  which  every  minimal  realization  is  controllable  and  observable. 
Hence.  Theorem  1.2.  part  4)  guarantees  that  for  any  two  minimal  realizations 


there  exist  nonsingular  sequences  {M^>  and  {N^}  which  relate  the  various 
matrices.  However,  it  may  not  be  possible  to  find  {M^>  and  {N^}  which 


converge.  Consider  the  sequence  given  by 


1 

k 


"k(a)  ■  i  , 

with  realizations  (-,-1,^,1)  and  (•*;, -1 . 1  ,£) .  We  have  immediately  .  so 


,Vk. 


However,  there  is  still  one  interesting  case  where  a  result  is  possible. 
Xil*2£®i!LJL.J5  Suppose  {H^}  is  a  convergent  sequence  in  IR(s)  and  let 
{(E^.A^.B^.C^)}:  i=l,2  be  two  minimal  realizations  of  {H^}  with 

(Eik’Aik*Bik’Cik^  c°htrollable  and  observable  for  every  i,k.  FUi *.  ier  assume 

that  each  (E^  ,A^  ,B^  ,C^ )*lim(E^, A^.B^.C^)  is  controllable  and  observable. 
Then  there  exist  nonsingular  matrix  sequences  {M^}  and  (Nk)  with  Mk~ *M  and 
Nk— »N,  M  and  N  nonsingular,  such  that  Mk  ElkNk=E2k'  MkAlkNk=A2k’  Vlk^k*  and 

ClkVC2k  f0r  8Very  *' 

Proof  Applying  the  decomposition  (7), (8)  to  l  (Eii<,Aik,Bik'Cik^  yields 

/■*  •*  **  v 

decomposing  matrices  M^*— •M^  and  and  decomposed  system  matrices  Algk> 

Aifk'  Bisk'  Bifk’  Cisk*  and  Cifk'  with  limAifk  nilP°tent>  This  determines  in 
two  ways  the  same  time-scale  decomposition  Hk=Hgk+H^k  given  by 
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isk(aI"Aisk) 

•in 

Bisk 

-1 

ifk|sAifk'n 

Bifk 

fk 

Note  that,  for  sufficiently  large  k,  each  of  the  subsystems  ^ 1 ^isk'^isk’^^ak^ 

and  (A^^.  I  must  be  controllable  and  observable.  From  [16. p.  208  1. 

T  -IT 

the  similarity  transformation  ’^s^<::(V2skV2sk,  V2skVlsk‘  where  vis(<  is  tlle 
observability  matrix  of  the  pair  (A^.C^),  takes  ( I  •Algki8lsj<'Clsk)  iat0 
( I -A2sk'U2sk’C2sk* '  Furthermore ,  {Tsk>  converges  to  the  nonsingular  matrix 


T  -IT 
T  MV,  V„  )  VI  V,  . 
s  2s  2s  2s  Is 


A  similar  construction  yields  Tfk—+Tf, 


straightforward  calculation  shows  that  the  sequences 


Mk  =  M2k 


rTsk  0 


1  fk 


Mlk'  *k 


=  N 


lk 


t",1  0 

sk 


P-1 

fkJ 


2k 


yield  the  desired  result.  J 

Our  final  result  follows  with  the  aid  of  Theorem  6.13. 

Corollary  6.16  If  { >  is  a  slow  sequence  with  convergent  CP  and 
{(Eik.Aik.Bik,Cik)} ;  i*1.2  are  any  two  minimal  realizations,  then,  for 
sufficiently  large  k.  there  exist  nonsingular  matrix  sequences  {Mk >  and  {Nk> 
and  nonsingular  matrices  M  and  N  such  that  *M,  N^— *N,  MicEncl\'E2j<' 

•VlkVA2k’  Vlk*B2k’  and  ClkVC2k  t0r  every  k- 


7.  Concluding  Remarks 

The  problem  discussed  in  this  paper  is  the  realization  of  convergent 
transfer  matrix  sequences  with  convergent  generalized  state-space  sequences. 
Just  as  state-space  sequences  may  be  decomposed  according  to  time-scale 
behavior,  a  time-scale  decomposition  for  any  rational  matrix  sequence  may  also 
be  achieved.  We  have  shown  that  convergence  of  the  CP  of  a  sequence  of 
rational  matrices  Is  a  crucial  ‘••sue  in  the  minimal  realization  problem.  It 
was  proved  that,  when  the  characteristic  polynomial  of  a  rational  matrix 


sequence  is  not  convergent,  the  rational  sequence  can  be  decomposed  into 
finitely  many  subsequences  in  such  a  way  that  each  subsequence  has  convergent 
CP.  Our  results  demonstrate  that  the  general  problem  can  be  reduced  to 
finitely  many  subproblems,  each  of  which  can  be  handled  using  a  simpler 
theory.  It  is  hoped  that,  our  results  will  complement  the  robustness 

literature  at  large. 
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Technical  Notes  and  Correspondence 


Linear  Compensator  Designs  Based  Exclusively  on 
Input-Output  Information  are  Never  Robust 
with  Respect  to  IJmnodcIcd  Dynamics 

J.  DANIEL  COBB 

Abttract—Wt  investigate  the  effects  of  unmodeled,  higher  order 
dynamics  or  parasitics  on  the  stability  of  linear  control  systems.  We  first 
describe  a  class  of  perturbations  of  a  given  state  equation  which  cannot  be 
distinguished  from  the  original  on  the  basis  of  input-output  measure¬ 
ments  alone.  Then  it  is  shown  that,  given  any  plant-compensator  pair, 
such  perturbations  of  each  system  cun  always  tic  round  which  destabilize 
the  closed-loop  configuration,  finally,  the  effect  of  destabilizing  pertur¬ 
bations  on  output  behavior  is  explored. 

I.  INTRODUCTION 

The  effects  of  high-frequency  or  parasitic  phenomena  on  closed-loop 
system  performance  have  long  been  studied.  A  popular  framework  for 
addressing  this  issue  has  been  that  of  singular  perturbation  theory  (see, 
e.g.,  III.  |2|).  The  point  of  view  that  parasitics  are  ultimately  connected 
with  unmodeled  plant  dynamics  has  become  quite  popular  in  recent  years, 
sometimes  with  surprising  consequences.  For  example,  it  was  shown  by 
Rohrs  et  al.  [8]  and  Ioannou  and  Kokotovic  f3]  that  high-frequency 
phenomena  can  lead  to  instability  in  adaptive  control  schemes.  Adaptive 
controllers  being  highly  nonlinear,  u  natural  question  to  ask  is  whether 
parasitics  could  have  a  similar  destabilizing  effect  on  control  systems 
which  are  based  on  linear  compensators.  This  was  answered  in  the 
affirmative  hy  Khalil  in  [4|  and  |5).  A  notable  effort  to  circumvent  these 
difficulties  in  the  case  of  linear,  time-invariant  systems  was  made  by 
Vidyasagar,  culminating  in  the  results  of  |6|  and  |7|. 

Our  work  is  most  similar  to  |7|,  but  differs  primarily  in  that  wc 
investigate  the  stability  of  a  closed-loop  system  when  both  the  plant  and 
compensator  are  perturbed.  The  idea  of  perturbing  both  systems  has  been 
largely  neglected  in  the  literature  (with  the  notable  exception  of  |6|),  even 
though  one  can  easily  make  a  strong  case  for  considering  such 
perturbations.  Indeed,  one  need  only  recognize  that  a  compensator,  like 
the  plant,  is  a  physical  system  governed  by  a  mathematical  model  which  is 
inherently  subject  to  uncertainty. 

In  light  of  examples  such  as  those  contained  in  (4)  and  |5],  even 
arbitrarily  small  model  errors  arc  to  be  feared  since  such  effects  have  the 
capability  of  destabilizing  a  system  just  as  certainly  as  larger  errors  do.  In 
fact,  those  examples  illustrate  that  in  some  cases,  small  errors  can  cause 
greater  instability  than  do  larger  ones. 

In  this  paper,  we  intend  to  show  that,  when  uncertainties  in  both  plant 
and  compensator  are  taken  into  account,  even  strictly  proper  compensa¬ 
tors  are  subject  to  parasitic  destabilization.  Hence,  properness  of  the 
compensator  is  really  not  the  pivotal  issue  here  as  it  is  in  j7|.  Wc  will 
show  that,  if  only  input-output  information  concerning  the  plant  and 
compensator  is  available,  robust  compensation  can  never  be  achieved. 

The  results  of  this  paper  are  by  nature  primarily  negative.  We  do  not 
claim  to  have  a  clear  understanding  yet  of  exactly  what  constitutes 
sufficient  information  for  robust  compensation,  although  we  do  mention  a 
possible  approach  to  finding  an  answer  in  Section  V.  It  is  hoped  that  our 
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results  will  stimulate  further  discussion  in  an  area  which  has  been 
neglected  by  alt  but  a  handful  of  researchers. 

II.  PRELIMINARIES 

We  study  systems  characterized  by  the  linear,  time-invariant  state 
equations 

x=Ax+Bu,y  =  Cx+Du  (1) 


and  perturbations  of  (1)  given  by 


[::][;]•[*  s]  [;]•[!]• 


where  the  submatrices  in  (2)  satisfy 

An  -AhA  2j  Aji-A,  B,-AnA  (3) 

C|-Ca4  J2  Aj,*C,  -CiA^Bi^D  (4) 

and  4j2  is  nonsingular.  If  we  set  c  =  0  in  (2)  and  eliminate  €.  (1)  is 
obtained;  hence,  (2)  with  e  =  0  may  be  thought  of  as  a  state  augmentation 
of  { I ).  Setting  <  >  0  in  (2)  constitutes  a  perturbation  of  that  augmentation. 
For  the  moment,  wc  allow  Aa  to  be  either  stable  or  unstable. 

To  aid  our  analysis,  we  will  use  the  decomposition  for  singularly 
perturbed  systems  developed  in  [  10|  where  it  is  shown  that  there  exist  real 
matrix-valued  analytic  maps  t  M,  and  t  >-  N,,  defined  on  some 
interval  |0, 0),  such  that  M,  and  N,  are  square  and  nonsingular  for  every  * 
and 

«[: [i  X.] x]w-[t  °] 

(5) 

with  A/,  and  A„  analytic  and  A/0  nilpotent.  According  to  [10],  the 
matrices  M,  and  N,  are  unique  up  to  change  of  bases;  hence,  we  may  take 
iVf0  and  No  to  be  any  matrices  which  achieve  the  decomposition  (5)  at  <  = 
0.  For  example,  let 


Next,  define 


[j;]  [j;]  .ic„  q.i-ic,  c2]Af,. 

i  (3)  and  (6)  yield  the  decoupled  state  equations 

(-][:]•[*!][:]•[*]• 


Equations  (3)  and  (6)  yield  the  decoupled  state  equations 


y=C„x,+Cr,xf 


[;]-*■  [•]■ 
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We  now  present  a  series  of  technical  results  which  will  be  useful  in 
Sections  III  and  IV. 

Lemma  I:  Al0  =  A,  B,a  =  B,  C,0  =  C,  0  ®  Cjo  »  C2An', 
and  A/,  =  tF,  for  every  e  €  [0,  0)  where  F0  =  A  *2'. 

Proof:  From  (5)  and  (6),  we  have 


A, a  0 
0  / 


-"•[>  ?]• 

[fc]  -*  [s] 

IC,0  C/01  =  |C,  C,|AJ0=|C  Cj/tjj'l. 


We  thus  have  tMn,  =  A/,fin,  so  Ajt  =  tF,  where  Fo  =  - 

/45j.  □ 

From  Aj,  —  tF, ,  we  immediately  obtain  the  well-known  result  that  the 
eigenvalues  of  (2)  which  tend  to  infinity  as  t  -*  0*  are  "close"  to  those  of 
(l/«)Aa  (see,  e.g.,  (2,  Corollary  2.1]).  One  useful  way  of  slating  this 
result  is  the  following. 

Lemma  2:  If  y  is  an  eigenvalue  of  Ah,  y  >  0,  and  R  <  oo,  then  there 
exists  (0  >  0  such  that  (2)  has  an  eigenvalue  X,  satisfying  |X,|  >  /(and 
jarg  \  -  arg  (l/«)ul  <  7  whenever  0  <  e  <  t0. 

Proof:  From  (7),  the  eigenvalues  of  (I /t)F~l  are  also  eigenvalues 
of  (2).  Since  F0" 1  *  A21  and  F~ 1  is  continuous  in  e.  each  F* 1  has  an 
eigenvalue  y,  with  Choose  t0  so  that  { I /<) |  y,  |  >  R  and 

|arg  y,  -  arg  ^1  <  7  whenever  0  <  <  <  t0,  and  let  X,  =  (1 /«)/»«•  Then 
X,  is  an  eigenvalue  of  (2),  |X,|  >  R,  and  |argX,  -  arg  (i/«)g|  =  |arg  y, 
-  arg  g|  <  y.  □ 

Suppose  the  transfer  matrices  of  ( I )  and  (2)  are  P  and  P, ,  respectively. 
We  will  need  conditions  under  which  an  eigenvalue  of  (2)  is  also  a  pole 
of  P,. 

Lemma  3:  If  (Ah,  Bi,  Cj)  is  controllable  and  observable,  there  exists 
(o  >  0  and  R  <  »  such  that  every  eigenvalue  X,  of  (2)  satisfying  |X,|  > 
R  is  also  a  pole  of  P,  whenever  0  <  t  <  <0. 

Proof:  An  eigenvalue  X,  of  (2)  is  a  pole  of  P,  if 


f"\/-/4n  -An  B,  _  w-t  \l~A„  0  B„  1 

[  -A„  t\I-Au  B2  j  [  0  t\,l-F~'  B/,  J 


have  full  rank.  Choose  R  >  max  { |  X}  (X  is  an  eigenvalue  of  A}.  From 
Lemma  I ,  (F„"  Bt 0,  C/oF"  1 )  »  (Au,  B%,  C2).  Hence,  there  exists  t0  > 
0  such  that,  whenever  0  <  t  <  to,  (F‘ ',  B/, ,  C/,F~ 1 )  is  controllable  and 
observable  and  |X,|  >  R  implies  that  X,  is  not  an  eigenvalue  of  A„.  It 
follows  immediately  that  all  matrices  on  the  right-hand  sides  of  (8)  and  (9) 
have  full  rank.  □ 


III.  Input-Output  Equivalence 

In  this  section,  we  explore  the  relationship  between  the  nominal  and 
perturbed  systems  (I)  and  (2)  and  discuss  the  conditions  under  which  they 
are  indistinguishable  if  only  input-oulput  information  is  available. 
Consider  the  process  of  obtaining  or  veritying  an  input-output  model  of  a 
physical  system.  We  are  allowed  to  take  measurements  by  applying  an 
input  signal  starting  at  t  -  0  and  by  observing  the  output;  it  is  assumed 
that  no  direct  access  to  internal  states  is  possible.  Once  a  nominal  model  is 
obtained,  a  controllable  and  observable  realization  can  be  chosen, 
yielding  the  state  equation  ( 1 ).  Since  we  have  no  direct  control  over  initial 
states  except  through  the  input  ports,  and  since  /  «  0  presumably  occurs 
long  alter  the  system  was  built,  the  system  may  be  assumed  initially  at 
rest.  Hence,  we  choose  t(0)  =  0  and  £(())  =  0  in  (I)  and  (2). 

We  define  the  class  of  admissible  input  signals  It  to  be  all  C1  functions 
u:(0,  r]  -»  i Rm  satisfying  max  ||u(/)||  <  Xo,  max  ||t)(/)||  <  K\,  and  u( 0) 
=  0  where  the  constants  r  <  »,  Kn  <  ®,  and  K,  <  «*>  are  independent 
of  u.  Prom  an  engineering  standpoint,  it  is  not  unreasonable  to  place  such 
restrictions  on  u.  indeed,  in  any  real-world  scenario,  there  is  a  maximum 
length  of  time  one  would  be  witling  to  invest  in  collecting  data,  as  well  as 
a  maximum  amplitude  of  voltage,  force,  or  other  input  quantity  that  could 
possibly  be  generated  using  available  technology.  Furthermore,  there  is 
always  an  upper  bound  on  the  rate  at  which  u(t)  can  be  made  to  vary  (e.g., 
every  amplifier  has  a  maximum  slew  rate).  Thus,  the  constants  r,  K0,  and 
Ki,  although  possibly  very  large,  must  be  finite.  Since  no  input  is  applied 
prior  to  /  *  0  and  since  K,  <  »,  we  must  have  u(0)  =  0.  We  would 
surely  be  in  serious  trouble  if,  in  order  to  design  a  robust  compensator,  we 
needed  the  capability  of  generating  inputs  over  arbitrarily  large  intervals 
of  time  oi  with  arbitrarily  large  amplitudes  or  rates  of  change. 

Associated  with  any  real-world  measuring  device  is  a  minimum  error 
which  can  be  detected.  For  example,  if  a  function  y  represents  an  output 
voltage,  velocity,  or  other  physical  quantity  of  interest,  there  must  exist  a 
number  5  >  0,  characteristic  of  the  measuring  device  alone,  such  that 
another  output  y  cannot  be  distinguished  from  y  if 


sup  {||>’(0-3,(/)|||  0s/St] <4.  (10) 


For  the  remainder  of  the  paper,  we  assume  a  fixed  source  of  input  signals 
and  measurements  and,  consequently,  a  fixed  set  Tl  and  number  6  >  0. 

The  quantities  11  and  4  together  determine  an  equivalence  between 
systems:  two  systems  are  indistinguishable  under  input-output  measure¬ 
ment  if  for  every  u  6  'll .  the  output  functions  y  and  f  of  the  two  systems 
satisfy  ( 10).  The  next  .esult  applies  this  idea  to  the  nominal  and  perturbed 
models  (1)  and  (2). 

Theorem  I:  If  An  is  strictly  stable,  there  exists  <o  >  0  such  that, 
whenever  u  €  11  and  0  s  t  <  t,„  the  respective  outputs  y  and  y,  of  (1) 
and  (2)  satisfy  max  {||y(r)  -  y,(0||  |0  s  t  £  r}  <  4. 

Proof:  We  first  note  that  yf.1)  *  JJ  C,0  exp  (riA,o)Bs0u(t  -  rt )  rf») 
-  CjnBfou(t)  =  y(().  Hence,  we  need  only  show  that  there  exists  to  such 
that  |jy,(r)  -  y0(Oll  <  5  whenever  0  s  l  s  r  and  0  S  e  <  q. 
Decomposing  y,  =  y„  +  y;,  in  the  obvious  way,  we  have  ||y„(/)  - 
yj0(r)||  s  Kol’0\\C„exp(i)A„)B„  -  Cj0exp(qAj0)j5,ol|  dt). Choose t, 
>  0  such  that  0  £  t  <  et  implies  max  {||C„  exp  (vA„)B„  -  C, o  exp 
(<jAiD)5,o||  |0  s  i|  s  r)  <  4/(2A'0r).  Integrating  by  parts,  we  obtain 


IMO-y/olOllsW/.ll 


d  v 


■  II  #/<  II  +  Ao  ||  C/,  5/,  -  C/o  fl/o  || . 


There  exists  >  OandA  <  oo  such  that  l|exp(/F',)||  <  K,  ||C/,||  <  K, 
and  ||F/,|  <  K  whenever  t  z  0  and  0  £  t  <  <j.  Let  $  *  4/(4A|A,J(A' 
+  r)).  We  know  that  there  exists  t j  >  Osuch  that  ||exp  ((jj/<)Fi'i)||  <  S 
whenever  5  s  y  s  r  and  0  <  e  <  <j  (see,  e.g.,  [13)).  Finally,  there 
exists  t,  >  0  such  that  ||  C/,B/,  -  C/n5/n||  <  S/4Rn  when  e  <  <4.  Let  to 
=  min  {ei,  ti,  <3,  <4}.  Then  0  s  t  <  t0  implies  By.f/)  -  >o(08  <  h/2 
+  K,KHK&  +  tS)  +  4/4  =  4.  □ 

Wc  have  thus  established  that,  for  sufficiently  small  <,  (1)  and  (2)  are 
indistinguishable  on  the  basis  of  input-output  information.  Hence, 
although  the  physical  system  is  nominally  described  by  (I),  an  equally 
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valid  model  from  an  input-ouiput  perspective  is  given  by  (2)  with  t 
sufficiently  small  and  /422  strictly  stable. 

IV,  Closed-Loop  Destabilization 

We  arc  now  ready  to  investigate  the  effects  that  the  system  perturba¬ 
tions  in  Section  II  have  on  a  closed-loop  configuration.  Consider  the 
feedback  compensator  governed  by 

z  =  Fz  +  Gy,  u  =  Hz  +  u.  (II) 

We  consider  only  compensators  with  strictly  proper  transfer  matrices 
since  the  results  of  [7]  indicate  that  nonstrictly  proper  compensators  are 
never  robust  with  respect  to  unmodeled  dynamics.  Perturbations  of  (I  I) 
are  of  the  form 

[-][?]■[£  a] 

u  =  H,z  +  Hit+v  (12) 


(12) 

!j  =  G 

(13) 

=  0 

(14) 

-HtFtt  Gi  =  0  (14) 

and  Fa  is  nonsingular.  The  discussion  of  Section  III  applies  equally  well 
to  both  plant  and  compensator. 

Combining  (I)  and  (13)  in  a  standard  feedback  configuration  yields 

DH-  [;]*[»]■ 

y  =  Cx+DHz.  (15) 

Combining  the  perturbed  systems  (2)  and  (12)  gives 


/  0  0  0 
0/00 
0  0  (I  0 
0  0  0  </ 


^11 

B,H, 

A  12 

B,H j 

C.C, 

F„ 

Q1C1 

F,i 

An 

B>H, 

A 22 

BiH> 

GtC, 

F,  1 

g*c2 

F» 

y 

=  C,jr+C2z. 

[  An  BiH> 
[CjCj  Fl:  _  ' 


'-[*]• 


observable,  and  X  is  nonsingular  with  an  eigenvalue  in  the  sector  S.  Then 
there  exists  f0  >  0  such  that  H,  has  a  pole  p,  6  S  satisfying  |p,|  >  R 
whenever  0  <  t  <  <o- 

Proof:  Since  X  is  nonsingular,  the  closed-toop  system  (16)  is  of  the 
form  (2).  Let  p  6  S  be  an  eigenvalue  of  X.  There  exists  y  >  0  such  that 
s  t  S  whenever  |arg  s  -  arg  (1/«)mI  <  y-  The  result  then  follows  from 
Lemmas  2  and  3.  D 

Now  consider  behavior  of  the  output  y(t)  in  the  closed-loop  system 
(16).  Theorem  3  shows  that  under  certain  conditions,  the  instability 
described  in  Theorem  2  also  has  a  pronounced  effect  on  y(l).  Let  m  denote 
Lebesgue  measure. 

Theorem  3:  Suppose  R  <  »,  5,,  5,  >  0,  (X,  Y,  Z)  is  controllable  and 
observable,  and  X  is  nonsingular  with  an  eigenvalue  in  the  sector  S. 

1)  There  exists  e0  >  0  such  *hat  corresponding  to  each  t  6  (0,  <«). 
there  exist  vectors  Xu,  £  IH",  Zo.  6  /?*,  fo,  €  /?'*,  fa,  €  iW11  with  ||.roJ> 
II Zo,  II .llto.ll.llfo.il  <  ^andu  set  0,  C  [0,  r]  with  m 0,  <  62  such  that  the 
output  y,  of  (20),  subject  to  xfO)  =  *0,,  z(0)  =  Zo«.  t(0)  =  to..  f(0)  = 
fo,,  and  u  *  0,  satisfies  ||y,(r)||  >  R  for  every  t  6  [0,  r)  -  0,. 

2)  There  exists  t0  >  0  such  that  corresponding  to  each  e  G  (0,  f0),  there 
exist  a  continuous  function  w,:|0,  r]  -*  !Rm  with  ||u,(/)||  <  6|  for  all  t  G 
|0.  r)  and  a  set  fl,  C  |0,  r)  with  ml),  <  52  such  that  the  output  of  (20), 
subject  to  x{Q)  =  z(0)  =  {(0)  =  f(0)  =  Oand  u  ■  u„  satisfies  ||y,(/)|| 
>  R  for  every  t  G  (0,  r]  -  0, . 

Proof:  I)  Since  R  is  arbitrary  and  the  system  (16)  is  linear,  we  need 
only  prove  the  result  for  a  single  vector  norm,  say,  the  Euclidean  norm. 
The  decomposition  (7)  may  be  applied  to  (16),  yielding  real-valued 
analytic  matrix  functions  •  •  •,  F,  defined  on  an  interval 

[0. 01.  Since  F0  =  X~ 1  is  nonsingular,  F~ 1  is  analytic.  It  is  shown  in  ( 151 
that  there  exists  a  continuous  complex  unitary  matrix-valued  function  t  - 
U,  defined  for  sufficiently  small  values  of  e  that  puts  F~'  into  continuous 
upper  triangular  form— i.e., 


Mi.  “i2,  ‘  •  * 

0  • 


U-'F-'U,^  :  .  •  •  • 

I  ,  ,  .  I 

L°  •••  0  #w.,  J 

where  each  of  the  maps  e  -*  p„  and  e  -» cty,  is  continuous.  Additional  row 
and  column  interchanges  can  be  used  to  reindex  the  p„\  equivalently,  U0 
may  be  chosen  so  that  #i,0  G  S. 

Let 


■  mi*  ..  0 

• 


Let  (15)  and  (16)  have  transfer  matrices  H  and  H„  respectively. 

From  this  point  on,  we  assume  that  A22  and  Fu  are  strictly  stable 
matrices.  Thus,  according  to  Theorem  1 ,  (2)  and  (12)  are  equivalent  to  ( 1) 
and  (11)  for  sufficiently  small  e  in  an  input-output  sense.  The  perturbed 
closed-loop  system  (16)  is  also  of  the  form  (2);  no  obvious  conclusions 
can  be  drawn,  however,  concerning  stability  of  either  (16)  or  the  matrix 


In  view  of  Lemmas  2-4  as  related  to  ( 16),  we  see  that  the  properties  of  X 
as  well  as  those  of  the  matrices 


Since  N,  is  nonsingular  on  [0,  0),  ||Af,||  is  nonzero.  Standard  norm 
inequalitites  reveal  that  ||lf',||  <  5t.  From  (7),  it  follows  that  the  natural 
response  of  (16)  due  to  the  initial  condition  w,  is 


3’.(0*2jj^|(C/.y.)exp<it/;,F;,f/1)  ?  .  (17) 

6 

*  «i 

From  Lemma  1,  (F"1,  C/0)  =  (X,  ZX-').  This  pair  is  observable  since 
X  is  nonsingular;  the  corresponding  observability  matrix  is 


are  crucial  for  understanding  the  behavior  of  (16). 

We  are  ultimately  interested  not  only  in  the  eigenvalues  of  the  closed- 
loop  system,  but  also  in  the  poles  of  H,  and  the  behavior  of  (he  system 
output  y(t).  The  next  two  results  treat  first  the  closed-loop  poles  and  then 
output  behavior.  As  a  means  of  quantifying  instability,  let  a  e  (0,  w/2)  and 
consider  the  open  sector  S  =  {s  €  f)  -  {0J|  |arg  s|  <  a}. 

Theorem  2:  Suppose  R  <  «,  (X,  Y,  Z)  is  controllable  and 


and  the  pair  ( X ,  Z)  is  observable.  Thus,  (Uq'Fq'Uo,  C/oU0)  is 
observable.  Since  t/0" 1  F0“ 1 U0  is  upper  triangular,  the  first  column  of 
C/0U0  is  nonzero.  Suppose  oo  *  0  is  the  /th  entry  of  the  first  column  of 
C/0U0.  Then  the  same  entry  a,  of  C/,  U,  is  nonzero  for  sufficiently  small 
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<.  From  (17),  it  follows  that  j>,  has  fth  entry  jv,«{/)  =  (i,/(2||A/,||»a,  exp 
((l/,)#»i ,/)•  Thus.  Re  jJ,  has  the  Ah  entry 

Rc _MO  =  (V(2||AMI))l<*.l  c*P  (;  Rc  Mi ,0  cos  (;  lm  pui  i-urg  a,). 

Since  n,0  G  S,  Re  pu  >  0  for  small «.  From  elementary  analysis,  there 
exists  <o  >  0  such  that  0  <  «  <  e0  implies  the  existence  of  a  set  0,  with 
«fl,  <  A}  anti  || Rc  MOW  >  *  for  all  t  G  |0,  r|  -  0, . 

We  note  that  the  initial  condition  w,  may  be  complex.  In  general,  the 
natural  response  of  (16)  is  of  the  form  y(t)  =  r,(/)w  where  V,(t)  isareal- 
valued  matrix.  Hence,  Re  y(t)  -  r,(/)  Re  w,  and  if  we  set  [Aro,Zo,$o<  fo.1 T 
~  Re  wf,  we  obtain  an  output  y,  with  rth  entry  y,,(t)  =  Re  9„(.t). 
Therefore,  ||y,(/)||  2  |jl„(/)|  >  R  for  all  t  G  (0,  r)  -  0,.  Finally,  we 
note  that  || atq,  || .  ||?o,ll.  llio.ll.  Ilib.ll  S  ||Re  w,||  <  ||w,||  < 

2)  Our  approach  is  to  construct  an  input  function  u,  which  steers  the 
system  (20)  from  the  origin  to  some  state  w,  satisfying  the  conditions  of 
part  1),  the  transfer  occurring  on  an  arbitrarily  small  t  interval;  then  the 
system  will  be  allowed  to  evolve  front  w,  with  zero  input.  We  first 
consider  the  pair  (F*1,  F~ 'Bfl).  From  Lemma  1,  (F0'‘,  F0"'fl/0)  =  (A", 
XY).  This  pair  is  controllable  since  X  is  nonsinguiar;  the  corresponding 
controllability  matrix  is 

l XY  XlY  Xl,'tY\=*X\Y  XY  •••  X^^'Y) 

and  (X,  Y )  is  controllable.  Hence,  (F'1,  F~ 1 B/, )  is  controllable  for 
sufficiently  small  e.  Let 

MO- B/,F,'r  exp  (-IF-T)W,  (r)-'  exp  (-rF;r)  (18) 

where  the  Gramian  W,(r)  is  given  by  =  jj  exp  (-i/F*1) 

F-lB/,Bj F~T  exp  (-r)F~T)  di).  W,  is  nonsingular  for  small  e  since 
(F~\  F~'Bf,)  is  controllable  (see  [It,  p.  184]).  All  matrices  in  (18) 
converge  and  exp  ( -  tF~ T)  converges  uniformly  on  (0,  r]  as  t  -*  0*; 
hence,  &  converges  uniformly  to  4>o-  Thus,  there  exists  a  number  M,  < 
»  such  that  HiMOII  <  Af|  for  all  t  €  [0,  r]  and  t  sufficiently  small. 

Choose  Mi  <  oo  such  that  ||  C„  exp  (M„)||  <  A/2  for  small  t  and  all  t 
€  [0,  r|  where  C„  and  A„  are  given  by  (7).  Since  N‘ '  is  continuous,  we 
know  from  part  I)  that  for  sufficiently  small  t.  there  exist  real  vectors  Xn. , 
Zo. .  $o«  •  and  fa,  with  ||jto,  || ,  l|Zo>  II >  II ko, I! .  II  fo.  II  <  «i/(2Af,||AM||)  and  a 
set  ft,  with  mft,  <  5j/2  such  that  the  corresponding  output  j>,  of  (16) 
satisfies  ||J>,(/)1|  >  R  +  ( Mi/M\)6 ,  for  every  t  G  [0,  r]  -  ft,.  Let 


Then  the  output^,  may  be  written  =  yu  +  y/t  where  >,,(<)  =  C,(  exp 
Xq,'  and y/((/)  =  C/t  exp  ((f/e)  F~l)  From  (19),  ||*o,,||  <  6,/ 
Mi;  tlierefore,  |ji',<(/)||  <  (A/j/A/|)5(  for  every  t  €  [0,  r).  It  follows  that 
ILy/,(OII  >  WOll  -  lb., (Oil  >  R  for  each  t  G  |0.  r]  -  ft,. 

Next,  define  0,(0  =  y,(0*o/..  Then  ||jr0/,||  <  6\/Mt  guarantees  that 
110,(011  <  4|.  and  u,  steers  the  systems  =  F~xx  +  F~'Bf,u  from  the 
origin  as  /  »  0  to  Xo/t  at  /  *  r.  (See,  e.g.,  [II,  p.  556],)  Let 


The  divergence  of  the  output  of  the  closed-loop  system  described  in 
Theorem  3  is  referred  to  in  analysis  texts  as  "almost  uniform  convergence 
to  infinity."  In  view  of  the  arbitrarily  tight  bounds  that  muy  be  placed  on 
an  input  or  initial  condition  which  generate  this  divergent  behavior,  we 
conclude  that,  if  the  assumptions  of  the  theorem  are  met,  unbounded 
instability  at  the  output  of  a  closed-loop  configuration  can  result  from 
arbitrarily  small  noise  impinging  on  the  system. 

So  tar  we  have  demonstrated  that  the  existence  of  destabilizing 
perturbations  of  the  plant  and  compensator  is  guaranteed  if  a  certain  linear 
algebra  problem  admits  a  solution.  Indeed,  if  any  /I221  and  C%  arc 
chosen,  (3)  and  (4)  may  be  satisfied  by  simply  selecting  A 12  and  An 
arbitrarily  and  solving  for  Au,  Bu  and  Cj.  A  similar  remark  applies  to 
(13)  and  (14).  It  is  sufficient,  therefore,  to  find  A  22,  Bit  Q ,  F22,  G 2,  and 
Hi  such  that  I)  An  and  Fu  are  strictly  stable,  2)  (X,  Y,  Z)  is  controllable 
and  observable,  3)  X  is  nonsingular  with  an  eigenvalue  in  S,  and  4)  (4) 
and  (14)  are  satisfied.  Theorems  2  and  3  further  indicate  that,  if  l)-4)  are 
met,  the  resulting  instability  in  (16)  becomes  progressively  worse  as «  -♦ 
0  since  R  may  be  chosen  arbitrarily  large.  Thus,  arbitrarily  small 
uncertainty  can  lead  to  arbitrarily  large  instability. 

We  now  address  the  linear  algebra  problem  l)-4).  We  really  need  to 
find  only  one  solution  in  order  to  demonstrate  the  existence  of 
destabilizing  perturbations;  however,  it  is  possible  to  do  better.  To  obtain 
an  understanding  of  just  how  many  destabilizing  perturbations  actually 
exist,  let  (I),  (2),  (I  I),  and  (12)  have  orders  n,  n  +  it,  k,  and  k  +  k, 
respectively;  define  q  =  (n  +  rt)  (n  +  A  +  m  +  p)  +  (k  +  k)  (k  +  tc 
+  p  +  m).  Also,  consider  the  variety  in  consisting  of  all  (/In,  •  •  •, 
Ci,  Fu,  •  •  -,  Hi)  such  that  (3)-(6)  and  (13)  and  (14)  are  satisfied,  and  let 
V  C  •il<7  denote  the  intersection  of  that  variety  with  the  subset  in  which 
An  and  F!2  are  strictly  stable.  V  may  be  interpreted  as  the  set  of  all 
possible  state  augmentations  of  (1)  and  (11)  of  order  A  and  k, 
respectively.  Finally,  let  P  C  Ji’  be  the  set  of  all  points  for  which  (X,  Y, 
Z)  is  controllable  and  observable  and  X  is  nonsingular  with  an  eigenvalue 
in  S.  We  are  interested  in  properties  of  the  set  KHP. 

Theorem  4: 

1)  V  D  P  is  relatively  open  in  V. 

2)  V  c  T  is  nonempty  if  k  2  2  and  cither  a)  D  =  0  and  A  2  2  or  b)  D 
*  0  and  A  >  rank  D. 

Proof: 

1 )  This  is  obvious  since  T  is  open  in  i H". 

2)  Suppose  0  =  0  and  consider 


'  2V(s+l)'  0  •••  O' 

*2*5/(j+  1)*  0  •••  o' 

T(s)  = 

0  . 

,  l/(s)  = 

0  : 

6  •  •  •  •  0 

6  •  •  •  •  0 

Let  (A  j2,  Bi,  C2)  and  (F22,  G 1,  Hi)  be  controllable  and  observable 
realizations  of  T and  U,  respectively.  Then  An  and  F22  are  strictly  stable, 
-  CiA^Bi  =  T(C)  =  D,  and  -HiF^'Gi  =  1/(0)  =  0.  Note  that  T 
and  U  have  degrees  A  and  k.  Since  (X,  Y,  Z)  has  transfer  function 


Then  t|w,(/)||  <  St  and  u,  steers  the  second  subsystem  in  (7)  from  the 
origin  at  f  =  0  to  xe;,  at  /  =  rr.  w,  also  steers  the  first  subsystem  in  (7t 
from  the  origin  to  some  state  Xf>„  at  t  =  tr.  Since  j?0j,  is  given  by  the 
convolution  integral  i0».  *  So'  e*P  OA„)B„u,{t)  dt,  the  construction  of 
u,  and  uniform  convergence  of  exp  (tA„)  guarantee  that  x0„  —  0  as «  -* 
0  * .  Hence,  j?„  -»  0  uniformly  on  (0,  r]  as  «  -*  0  *  where  >„(/)  =  C„ 
exp  (lAt,)Jco„.  Applying  the  input  u,  steers  the  system  (20)  to  w,  = 
Af,[^,„j£o/,]rat  t  *  tr.  For  t  G  ]<r,  r],  the  corresponding  output  is y,(t) 
=  9,,{t  -  tr)  +  yf,(t  -  tr),  so  ||y,(0N  >R~  IIA.(f  -  <r)||  >  R  for 
small  <  and  all  t  G  (<r,  r]  -  («r  +  ft,).  Thus,  if  we  choose  e0  sufficiently 
smell  with  <0  <  6il2,  and  l),  =  [0,  tr)  U  (<r  +  ft,),  we  obtain  mfi,  <  S} 
and  i>«(0|  >  R  for  all  t  G  (0,  r]  -  ft,  whenever  0  <  <  <  t0.  □ 


K(s)  =  (/-r(s)(/(s))-T(s) 

~2ts(s+  l),'/((s+  1)''  *  (-2"*(s)  0  0* 

0  . 

6  •  ■  •  0 

and  V  has  characteristic  polynomial  A(s)  =  (s  +  l)***  -  2’l*(s,  it 
follows  that  {X,  Y,Z)  is  controllable  and  observable  and  X  is  nonsingular 
with  a  unit  eigenvalue. 

Now  suppose  D  *  0.  There  exist  nonsingular  matrices  M  and  N  such 
that 

[!>' !] 


IEEE  TRANSACTIONS  ON  AUTOMATIC  CONTROL.  VOL  33.  NO  ft.  JUNE  1988 


563 


where  r  =  rank  D.  Lei 

*l/(s+ty''  ’  1 

“ 

I/(s+  1) 

0 

7(s)  =  M  1 

1/(4+  1) 

0 

-r 

N 


U(s)  =  N 


2*'*  '• 'j/(jr+l)‘  0  •••  O' 
0 


M. 


Then 
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has  characteristic  polynomial  A(s)  =  (s  +  1)'  ‘((s  -t-  I)"'*  ',l  - 
2* '$).  Reasoning  similarly  as  lor  pan  a),  we  conclude  that  A22  and 
Fy  are  strictly  stable,  (4)  and  (14)  hold,  ( X ,  Y,  Z)  is  controllable  and 
observable,  and  X  is  nonsingular  with  a  unit  eigenvalue. 

To  complete  the  proof,  we  need  only  choose  A (2,  /t2(,  F|2,  and  F}\ 
arbitrarily  and  solve  for  the  remaining  matrices  from  (3).  (4)  and  (13), 
(14).  □ 

Part  I)  of  Theorem  4  demonstrates  that,  in  a  certain  sense,  the  high- 
frequency  effects  which  bring  about  closed-loop  instability  do  not 
correspond  to  the  complement  of  a  generic  set,  and  hence  cannot  be 
dismissed  as  merely  a  pathological  case. 

V.  Conclusions 

We  have  shown  that  input-output  information  alone  is  insufficient  for 
designing  robust  linear  compensators.  This  conclusion  leads  one  immedi¬ 
ate!  i  to  ask  what  further  information  is  actually  required  to  allow  a  robust 
design.  Although  we  cannot  give  a  clear  answer  yet,  we  can  offer  some 
insight.  The  development  of  our  results  indicates  the  high-frequency 
behavior  in  (2)  and  (12)  plays  a  role  in  destabilization.  Such  behavior  is 
closely  related  to  the  infinite-frequency  structure  of  (2)  and  ( 12)  with  t  = 
0  (see,  e.g.,  (14)).  One  might  therefore  suspect  that  some  knowledge  of 
the  poles  and  zeros  at  infinity  in  cither  the  plant  or  compensator  is 
essential.  The  exact  form  of  such  information  and  whether  it  can  be  easily 
measured  are  important  topics  for  further  research. 
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D.  W.  LUSE 

Abstract— A  multivariable  Nyquist  criterion  is  given  which  applies  to  a 
wide  class  of  open-loop  transfer  matrices,  including  tome  which  are 
improper  or  have  Infinitely  many  poles  in  Ihe  right-half  plane.  Applica¬ 
tion  of  the  test  to  specific  examples  requires  only  elementary  .knowledge 
of  complex  variable  theory. 

I.  INTRODUCTION 

The  Nyquist  test  (6]  for  closed-loop  stability  of  systems  in  terms  of 
open-loop  properties  has  been  generalized  for  applicability  to  many  types 
of  feedback  situations.  For  linear  systems,  the  main  three  directions  of 
research  have  been  toward  multivariable  systems  as  originated  by 
MacFarlanc  |l|,  distributed  systems  (e.g.,  |2J-|51).  and  nonscalar  gain 
variations  (e.g.,  |7]).  The  more  sophisticated  stability  tests  often  Involve 
technical  concepts  such  as  the  Fredholm  index  and  almost  periodic 
functions.  The  test  given  in  this  note  applies  to  a  very  wide  class  of 
transfer  matrices,  but  the  stability  conditions  themselves  involve  only  the 
winding  numbers  of  simple,  dosed,  piecewise-sniooth  curves.  This 
allows  a  number  of  theoretical  difficulties  to  be  ignored  in  the  theorem 
statement,  and  to  be  evaded  in  practice  for  many  particular  examples. 
These  advantages  are  gained  through  (wo  sacrifices  of  hypothesis.  When 
stability  of  distributed  systems  is  considered,  the  type  (or  types)  of  L„ 
stability  is  a  concern  |3).  [I0|,  For  this  note,  stability  is  restricted  to  be 
finite-gain  Li  stability.  The  second  sacrifice  is  that  the  open-loop  transfer 
matrix  must  be  evaluated  at  some  points  in  the  open  right-half  plane. 

Section  II  of  this  note  states  and  briefly  proves  a  version  of  the  Nyquist 
criterion.  Also  included  are  four  simple  examples.  Three  of  these  are 
pathological  examples  to  illustrate  the  strengths  of  the  theory,  while  the 
fourth  one  shows  why  (idealized  models  of)  feedback  loops  containing 
wave  propagation  behavior  are  not  robust  to  small  time  delays. 

If.  nyquist-Type  Stability  Test 

A  number  of  definitions  from  distribution  theory  are  needed  at  this 
point.  These  allow  for  a  brief  and  general  theorem  statement.  The 
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.jlly  sta&la  se'it grcuc.  Special  attention  will 
3e  given  to  relations  Between  the  hypotheses  used 
for  the  convergence  results  and  convergence  pro¬ 
perties  of  the  transfer  functions  for  the  approx- 
!^ating  systems. 
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Approximation  techniques  for  parabolic  control 
systems:  a  variational  approach 

In  this  paper  we  consider  the  linear  quad¬ 
ratic  regulator  problem  for  a  class  of  boundary 
control  problems  for  parabolic  systems.  The 
problem  is  formulated  using  a  variational 
approach  and  an  approximation  theory  is  devel¬ 
oped  for  solutions  of  the  associated  operator 
Riccati  equation.  Our  study  includes  strongly- 
damped  elastic  systems. 
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Uniform  stabilization/ exact  controllability  from 
the  boundary  of  thermoelastic  plates 

We  consider  the  small  vibrations  of  a- thin 
homogeneous,  thermally  and  elastically  isotropic 
plate  of  uniform  thickness.  It  is  known  that 
thermal  dissipation  alone  is  sufficient  to  strong 
ly,  but  not  uniformly,  stabilize  the  elastic  com¬ 
ponents  of  the  motion.  It  will  be  shown  that 
theintroductlon  of  additional  dissipation  through 
the  action  of  bending  and  twisting  moments,  and 
shear  force,  at  the  boundary  leads  to  uniform 
asymptotic  stability  and  to  explicit  asymptotic 
energy  estimates.  The  more  difficult  question  of 
exact  controllability  of  the  elastic  dynamics  by 
means  of  such  boundary  forces  and  moments  will 
also  be  considered.  It  may  be  shown  that  if  the 
thermal  diffusi vi ty  is  sufficiently  small,  the 
reachable  set  of  the  elastic  components  of  the 
motion  is  essentially  the  same  as  for  purely 
elastic  plate  motion. 
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On  the  Topological  Structure  of  the  Orbit  Space 
of  Controllable  Generalized  Linear  Systems 

Aba  tract*  We  study  the  topology  of  the  orbit 
space  of  controllable  descriptor  systems  modulo 
restricted  system  equivalence.  Using  a  scaling 
action,  we  prove  that  this  space  is  an  analytic 
manifold.  Using  the  Welerstrass  decomposition, 
we  obtain  an  analytic  stratification  of  this 
manifold.  By  decomposing  the  strata  into 
generalized  Hermits  cells,  and  using  tools  from 
Borel-Moore  homology,  we  compute  the  singular 
homology  groups  for  this  space  in  the  complex 
cases.  Consequently,  the  orbit  space  of 
controllable  descriptor  systems  is  a  smooth 
compactlf icatlon  of  the  orbit  space  of  controll¬ 
able  state  space  systems  modulo  change  of  basis 
in  the  state  space. 
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A  Oeoaetrlc  Decomposition  Theorem  for  Singularly 
Perturbed  Systems  with  Applications  to  Convergence 
of  Solutions 

Let  E  and  A  be  nxn  natrix-valued  analytic 
functions  of  a  nonnegative  real  parameter  t  with 
sEU)-A(e)  a  regular  matrix  pencil  for  all  e.  The 
limiting  behavior  of  the  solutions  of  the  faaily 
of  differential  equations 

E{«)x»A(e)x  (1) 

at  «— *0*  is  at  best  only  partially  understood. 
Soae  sufficient  conditions  are  known  for  certain 
types  of  convergence,  but  no  necessary  and 
sufficient  conditions  are  known  for  nontrivial 
topologies-  We  first  present  a  new  geometric 
decomposition  theorem  for  ( 1 ) .  We  then  show  that 
this  result  enables  us  to  prove  a  necessary  and 
sufficient  condition  under  which,  for  any  initial 
vector,  the  corresponding  solutions  of  (1) 
converge  uniformly  on  compact  subsets  of  (0.») 
with  bounded  peaking  in  [0,1]. 
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